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EXECUTIVE  SUMMARY 


Orthogonal  frequency-division  multiplexing  (OFDM)  is  a  digital  communication 
technique  that  employs  the  simultaneous  transmission  of  multiple  sub-carriers  at  low  date 
rate  to  achieve  an  overall  high  data  rate  in  a  spectrally  efficient  manner.  While  using  low 
data  rate  sub-carriers  is  effective  in  multi-path  environments  and  in  combating  frequency 
selective  fading,  OFDM  systems  are  sensitive  to  loss  of  synchronization.  Lack  of 
synchronization  in  frequency  or  time  results  in  a  noticeable  degradation  in  performance 
of  the  OFDM  communication  system.  Lack  of  frequency  synchronization  occurs  due  to 
oscillator  mismatch,  phase  noise,  and  Doppler  shift,  resulting  in  loss  of  orthogonality 
between  the  sub-carriers.  Lack  of  timing  synchronization  results  in  missed  signals  due  to 
failure  to  detect  and  then  properly  demodulate  the  arriving  packet.  To  maintain 
synchronization  with  the  transmitter,  the  OFDM  receiver  must  therefore  detect  the  packet 
(packet  or  frame  synchronization),  estimate  and  correct  any  carrier  frequency  error, 
perform  phase  synchronization  and  achieve  accurate  symbol  timing. 

In  this  thesis,  the  sensitivity  of  OFDM  to  lack  of  synchronization  is  examined, 
and  several  synchronization  algorithms  are  evaluated  for  a  commercial,  widespread 
application  of  OFDM  -IEEE  802. 1 1  g  Wireless  Local  Area  Networks  (WLANs).  The 
governing  802.11  IEEE  standards  do  not  promulgate  the  type  of  algorithms  that 
commercial  vendors  may  use  to  obtain  and  maintain  synchronization  with  the  transmitted 
signal.  Instead,  training  symbols  and  pilots  are  mandated  by  the  IEEE  standard  as  part  of 
the  802. 1  lg  packet,  and  each  vendor  may  implement  the  most  efficient  algorithm  to 
achieve  the  best  performance,  in  terms  of  reception  range  and  data  rate.  Vendors  who 
leverage  these  mandated  symbols  and  pilots  employ  data-aided  synchronization,  and 
because  the  training  symbols  are  located  at  the  start  of  the  802.1  lg  packet,  this  allows  for 
quick  synchronization  in  a  WLAN  environment.  The  underlying  premise  of  this  thesis  is 
that  the  use  of  improved  synchronization  algorithms  will  result  in  greater  performance  - 
including  greater  detection  range  of  the  OFDM  signal.  In  the  following  pages,  no  direct 
comparison  is  made  to  any  commercial  synchronization  algorithms  due  to  their 
proprietary  nature  -  instead  promising  algorithms  that  result  in  detection  of  the  802. 1 1  g 


packet  at  low  signal-to-noise  ratio  (SNR)  to  simulate  extended  range  are  examined. 
Additionally,  promising  frequency  synchronization  algorithms  are  evaluated  as  a 
necessary  step  to  successfully  demodulate  the  OFDM  signal  in  the  low  SNR 
environment.  Successful  detection  and  frequency  synchronization  of  the  OFDM  signal  in 
a  low  SNR  environment  points  the  way  to  follow  on  research  in  which  a  direct 
comparison  between  the  promising  algorithms  and  current  commercial  synchronization 
algorithms  can  be  made  using  commercial-off-the-shelf  (COTS)  equipment. 

The  synchronization  algorithms  examined  in  this  thesis  are  designed  in  MATLAB 
and  then  subject  to  a  performance  evaluation  using  multiple  test  signals  at  low  signal-to- 
noise  ratio  to  simulate  detection  and  carrier  frequency  estimation  and  correction  at 
extended  range.  Due  to  the  necessity  to  synchronize  quickly  in  a  WLAN  environment 
and  because  the  training  symbols  must  be  transmitted  in  accordance  with  the  IEEE 
standard,  the  evaluated  algorithms  are  data-aided  and  focus  on  the  training  symbols  found 
at  the  beginning  of  the  802. 1  lg  packet. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Wireless  Local  Area  Networks  (WLANs)  are  rapidly  proliferating  throughout 
commercial,  governmental,  and  educational  organizations.  To  satisfy  increased 
consumer  expectations,  including  multimedia  content,  engineers  must  strive  to  deliver  the 
greatest  data  rates  possible,  using  robust  methodologies,  while  minimizing  costs  and 
bandwidth.  Unfortunately,  as  data  rates  increase,  wireless  transmissions  encounter 
increasing  problems  with  fading  [1].  To  combat  fading  problems  at  high  data  rate,  multi¬ 
carrier  modulation  is  used.  Multi-carrier  modulation  is  based  on  the  parallel  transmission 
of  data  using  sub-carriers  and  can  be  done  in  a  spectrally  efficient  manner  using 
orthogonal  frequency-division  multiplexing  (OFDM)  [2],  Consequently,  OFDM  has 
become  the  modulation  of  choice  for  many  wireless  networks. 

The  relevant  governing  standard  is  the  IEEE  802.11  standard  [3],  including  the 
IEEE  802.11a  Medium  Access  Control  (MAC)  and  Physical  Layer  (PHY)  specification 
[4]  and  the  IEEE  802.1  lg  amendment  governing  the  further  higher  data  rate  extension  in 
the  2.4  GHz  unlicensed  Industrial,  Scientific  and  Medical  (ISM)  band  [5],  The  IEEE 
802.1  lg  standard  allows  data  rates  of  up  to  54  Mbps  (million  bits  per  second)  at  strong 
signal  to  noise  ratio  (SNR)  [5]  but  the  range  of  these  WLANs  is  limited  due  to  imposed 
power  constraints  [5],  Performance  of  the  OFDM  system  often  is  also  impeded  by 
sensitivity  to  loss  of  synchronization  between  the  transmitter  and  receiver  [6]. 

B.  THESIS  OBJECTIVE 

The  goal  of  this  thesis  is  to  analyze  promising  synchronization  techniques  that  can 
be  used  to  extend  the  range  and  performance  of  the  802. 1  lg  WLAN.  These  promising 
synchronization  techniques  are  evaluated  in  a  low  SNR  enviromnent  to  simulate 
reception  of  the  OFDM  signal  at  extended  range.  If  successful  detection  of  the  signal  is 
accomplished  then  the  packet  detection  synchronization  algorithm  is  judged  to  be 
promising  and  may  be  reserved  for  future  analysis  in  direct  comparison  with  existing 
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commercial  synchronization  algrithms.  The  same  standard  is  applied  to  promising 
frequency  synchronization  algorithms,  a  necessary  step  in  the  successful  demodulation  of 
the  OFDM  signal. 

C.  RELATED  WORK 

The  OFDM  system  is  vulnerable  to  synchronization  errors  [7].  At  long  range,  and 
consequently  low  signal  to  noise  ratio  (SNR),  the  OFDM  receiver  must  reliably  detect  the 
incoming  packet,  and  synchronization  is  the  first  task  perfonned  [8].  Moreover,  lack  of 
frequency  synchronization  between  transmitter  and  receiver  can  contribute  to  lower 
signal  strength  and  destructive  interference  from  adjacent  carriers  at  the  receiver,  further 
reducing  performance  [9]. 

The  use  of  multi-carriers  allows  the  OFDM  signal  to  achieve  high  data  rates  in  the 
presence  of  multi-path  propagation  of  the  transmitted  signal,  but  also  exposes  the  OFDM 
signal  to  vulnerabilities  with  respect  to  frequency  offset  and  phase  noise.  This  loss  of 
frequency  synchronization  leads  to  a  phenomenon  known  as  inter-carrier  interference 
(ICI)  [8].  Additionally,  the  multi-carrier  system  must  address  packet  detection  and 
symbol  timing  in  the  context  of  low  SNR  at  extended  range.  Failure  to  properly 
synchronize  in  time  with  the  OFDM  symbol  will  result  in  inter-symbol  interference  (ISI) 
[8].  Therefore,  to  achieve  optimum  performance,  the  OFDM  receiver  must  do  both  - 
time  and  frequency  synchronization.  The  IEEE  standard  does  not  specify  the  specific 
synchronization  techniques  to  be  used,  so  a  variety  of  time  and  frequency 
synchronization  techniques  will  be  analyzed  at  low  SNR  using  MATLAB. 

D.  THESIS  ORGANIZATION 

Chapter  II  reviews  the  problems  of  fading  in  the  presence  of  multi-path 
propagation,  OFDM  signal  generation  and  various  critical  parts  of  the  OFDM  signal 
structure.  Chapter  III  addresses  in  detail  the  difficulties  of  receiving  the  OFDM  signal  in 
the  presence  of  ISI  and  ICI.  Chapter  IV  comparatively  analyzes  various  time  and 
frequency  synchronization  techniques.  Finally,  Chapter  V  reviews  analyzed  data  and 
provides  a  summary  of  results. 
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II.  BASICS  OF  ORTHOGONAL  FREQUENCY-DIVISION 
MULTIPLEXING  (OFDM) 

A.  THE  PROBLEM  OF  FADING  AND  HIGH  DATA  RATES 

Fading  of  the  transmitted  signal  occurs  at  the  receiver  when  there  is  movement 
between  transmitter  and  receiver  (Doppler  shift)  and  multi-path  propagation.  Multi-path 
effects  are  a  natural  phenomenon  caused  by  reflectors  and  scatterers  within  the  wireless 
channel.  Multiple  “copies”  of  the  original  transmitted  symbol  arrive  at  the  wireless 
receiver.  These  “copies”  or  echoes  then  add  destructively  or  constructively  with  the 
original  transmitted  symbol.  [1] 

The  delay  between  the  first  reception  of  the  transmitted  symbol  and  the  last  echo 
is  known  as  the  delay  spread  of  a  particular  channel  and  simply  represents  the  time 
difference  between  the  shortest  (usually  direct  line  of  sight)  and  longest  transmission  path 
[10].  Equation  (2.1)  shows  the  relationship  between  delay  spread  (cr)  and  difference  in 
transmission  distance  (Ad) : 

A  d  =  ac  (2.1) 

where  c  is  the  speed  of  light.  Due  to  longer  transmission  distances,  typical  RMS  values 
of  delay  spread  in  an  outdoor  environment  are  on  the  order  of  microseconds  while  indoor 
environments  have  average  values  of  100  nanoseconds  but  could  exceed  1500 
nanoseconds  in  a  worst  case  situation  [1]. 

A  smaller  delay  spread  is  advantageous.  When  the  delay  spread  is  a  significant 
fraction  of  the  transmitted  symbol  duration,  the  multi-path  effects  can  substantially 
interfere  with  the  follow-on  symbol.  In  this  case  of  “inter-symbol  interference”,  the 
received  pulse  is  unevenly  distorted  across  its  bandwidth,  a  situation  known  as  frequency 
selective  fading.  [1] 

Inter-symbol  interference  (ISI)  is  exacerbated  by  high  data  rates  required  by 
multimedia  content.  This  is  because  increased  data  rates  are  often  implemented  by  using 
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shorter  duration  symbols,  resulting  in  greater  symbol  bandwidth  [1].  If,  for  example, 
binary  phase  shift  keying  (BPSK)  is  the  preferred  modulation  in  an  indoor  wireless 
transmission,  then: 


where  W  is  the  signal  equivalent  noise  bandwidth  in  Hz,  R  is  the  data  rate  in  bits  per 
second,  and  T  is  the  bit  duration  in  seconds.  If  20  Mbps  data  rate  is  desired,  20  MHz 
bandwidth  is  required  and  bit  duration  is  50  nanoseconds.  In  this  case  bit  duration  is  less 
than  average  delay  spread  values  indoors  therefore  such  a  communications  link  would 
experience  substantial  ISI.  To  reduce  ISI,  the  solution  is  either  to  reduce  the  data  rate 
(increasing  bit  duration)  [1],  to  provide  channel  equalizers  to  mitigate  the  ISI,  or  switch 
to  an  M  -  ary  type  of  modulation  [11].  Channel  equalizers  result  in  increased  receiver 
complexity,  M  -  ary  modulation  requires  higher  SNR  per  bit  [11],  and  increasing  bit 
duration  (thereby  decreasing  bandwidth  and  decreasing  data  rate)  sufficiently  to  eliminate 
ISI,  would  require  bit  durations  much  larger  than  the  RMS  values  of  delay  spread.  If,  in  a 
particular  indoor  channel,  cr^  =  1 00  x  ]  0  9 ,  then  to  effectively  mitigate  ISI  for  a  binary 
modulation  scheme,  bit  duration  must  adhere  to  the  following  rule  of  thumb:  [1] 

r>10^s  (2.3) 

In  such  a  case,  bit  duration  can  be  no  less  than  1000  nanoseconds  implying  a  maximum 
data  rate  of  1  Mbps  vice  the  desired  20  Mbps.  Thus,  ISI  imposes  an  upper  limit  on 
symbol  rate  and  consequently  on  data  rate  for  BPSK  and  various  other  modulations. 

B.  FREQUENCY-DIVISION  MULTIPLEXING  AND  DATA  RATES 

In  our  previous  example,  if  we  instead  use  N  sub-carriers  vice  a  single  carrier, 
and  transmit  data  in  parallel  on  all  of  the  sub-carriers,  the  sub-carrier  bit  duration  is  now: 

Th=NT  (2.4) 

The  resulting  increase  in  sub-carrier  bit  duration  reduces  individual  sub-carrier 
data  rates  by  a  factor  of  A  but  the  parallel  transmission  of  all  sub-carriers’  signals  yields 
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the  original  desired  data  rate.  Thus  each  sub-carrier’s  signal  experiences  negligible  ISI  in 
a  multi-path  environment  [12].  IEEE  802.1  lg  and  several  other  wireless  standards  utilize 
this  multi-carrier  technique  to  solve  the  ISI/frequency  selective  fading  problem. 

C.  SPECTRALLY  EFFICIENT  SPACING  FOR  SUB-CARRIERS 

Multi-carrier  modulation  (based  on  multiple  carriers,  called  sub-carriers)  is 
essentially  frequency  division  multiplexing  (FDM)  [13].  In  order  for  this  scheme  to  work 
effectively,  so  that  there  is  no  interference  between  sub-carriers,  the  sub-carriers  must  be 
spaced  apart  in  frequency.  Figure  1  depicts  a  scheme  that  has  very  low  interference 
between  sub-carriers,  but  uses  too  much  bandwidth,  i.e.  it  is  not  spectrally  efficient. 


Frequency  (arbitrary  units) 


Figure  1  Frequency  Division  Multiplexing 


With  respect  to  spectral  efficiency,  there  is  a  better  way  to  manage  the  allotted 
bandwidth.  If  we  use  Orthogonal  Frequency  Division  Multiplexing  (OFDM),  there 
would  be  significant  savings  in  bandwidth  while  still  maintaining  our  lower  bit  rate  per 
sub-carrier.  The  resulting  power  spectral  density  (PSD)  for  the  same  eight  sub-carriers 
depicted  earlier  would  look  like  Figure  2.  The  increased  spectral  efficiency  compared  to 
simple  frequency  division  multiplexing  (previous  figure)  is  obvious. 
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Frequency  (arbitrary  units) 

Figure  2  Orthogonal  Frequency  Division  Multiplexing 


Note  that  in  Figure  2,  each  sub-carrier  peak  coincides  with  nulls  in  the  power 
spectral  density  of  all  other  sub-carriers.  This  is  a  practical  result  of  orthogonality.  In 
OFDM,  orthogonality  is  maintained  between  the  sub-carrier  channels  by  properly  spacing 
the  sub-carrier  center  frequencies.  Each  sub-carrier  channel’s  signal  is  a  quadrature 
amplitude  modulated  (QAM)  signal  given  as  [2] 

sf(t)  =  AdL  (t)cos (coj)  -  Ad0  (t)sin(0.t)  (2.5) 

where  A  is  a  constant,  proportional  to  signal  amplitude,  coi  is  the  center  frequency  for 
sub-channel  i ,  and 

=+oo 

d,_  (0  =  X  dii,nPTs  0  -  nTa).  (2.6) 

n=— oo 

The  value  of  dL  can  take  on  the  values +1, ±3, ±5,..., ±>/M-l;  Mis  the  number  of 
symbols  for  a  given  modulation  (M-QAM);  Ts  is  the  symbol  duration; 


J 1  if  0  <  t  <  7) 
[0  otherwise; 


(2.7) 


and  the  definition  of  dQ  (t)  is  analogous  to  that  of  d,  (Q  as  the  quadrature  component. 
Using  the  trigonometric  identity  [14]: 
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cos(a  +  /?)  =  cos(a)  cos  /3  -  sin(a)  sin  /? 


st(t)  may  be  re-written  as: 


st{t)  =  a,.(t)cos[(a,.t  +  9:(t)\ 


Where  ai(t)  =  AJdIi  (t)  +  d0i  (t)  and  9,(0  =  tan 


2(?)and  #(t)  =  tan  1  .  To  maintain 


(0 


orthogonality  between  sub-carriers,  the  correlation  between  the  sub-carriers  must  be  zero. 
The  correlation  between  signals  on  sub-carriers  p  and  q  is 

+oo  +oo 

cP,q  =  J  sp(t)sg(t)dt=  J  ap(t)cos[copt  +  0p(t)]aq(t)cos[coqt  +  9q(t)]dt 


H=+0o("+1)rs 

=  E  f  ap(t)cos[copt  +  9p(t)]aq(t)cos[(oqt  +  9q(t)^it 


n=+ oo 

-s  1  aP,n  cos [copt  +  epn\aqn  cos [coqt  +  dq  n(t)\lt 


(2.10) 


where  apn  =  A^jd,^  +  dQ„  ,  0p  n  =  tan“ 


and  dj  and  dQ  are  the  respective 


in-phase  and  quadrature  OFDM  symbol  values  of  the  p‘h  sub-carrier  over  the  interval 
nTs  <t  <(n  +  l)Ts.  Therefore,  using  the  following  trigonometric  identity  [14], 

cosacos/?  =  ^-[cos(a  +  /?)  +  cos(a  —  /?)], the  correlation  may  be  re-written  as 


-j  «=+oo  (w+l)^ 

CP.«  aP,naq,n  J  cos[(®p  +  (oq)t  +  (9pn  +  9q  „)]  +  cos[(<yp  -<oq)t  +  (9pn  -  9qj]dt  (2.1 1) 


if  (<u(.  +  ojj)  »  — ,  as  is  typically  the  case,  then 


n-+  oo 


j  cos[(fi>,-fi>9)f  +  (0^-09>II)]A  (2.12) 

^  «=-00 
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=  0  if  [cop  -  a>q )  Ts  =  Ink,  where  k  e  set  of  integers. 

Considering  adjacent  sub-carriers,  letting  k  =  1  for  closest  spacing,  and  thereby 
maximizing  spectral  efficiency, 

f  _f  _*Vi  cop-1-R 
Jp+1  Jp  In  2  n  Ts 

The  above  result  shows  that  to  achieve  orthogonality  between  the  sub-carriers, 
while  maximizing  spectral  efficiency,  the  sub-carrier  center  frequency  spacing  must  be 
equal  to  the  OFDM  symbol  rate. 

D.  INVERSE  DISCRETE  FOURIER  TRANSFORM  KERNEL  FOR  OFDM 

To  properly  transmit  the  OFDM  symbols  in  a  spectrally  efficient  manner,  as 
discussed  in  the  previous  section,  it  is  essential  to  create  and  maintain  orthogonality 
between  the  multiple  sub-carriers.  The  previous  section  demonstrated  that  the  greatest 
spectral  efficiency  that  could  be  achieved  while  maintaining  orthogonality  resulted  when 
sub-carrier  spacing  in  frequency  was  equal  to  the  sub-carrier  symbol  rate,  Rs.  The  most 
obvious  solution  might  be  to  have  a  transmitter  and  receiver  each  having  N  local 
oscillators.  Each  oscillator  would  tune  an  assigned  individual  sub-carrier  frequency  and 
the  oscillators  would  be  separated  in  frequency  by  an  amount  equal  to  the  sub-carrier 
symbol  rate.  Such  a  scheme,  however,  would  be  unreasonably  expensive  and  complex 
[2].  The  inevitable  relative  drift  in  oscillator  frequencies  would  result  in  non-zero 
correlation  between  the  sub-carriers,  destroying  orthogonality,  and  resulting  in  ICI. 

There  is  a  superior  method  that  ensures  that  all  oscillator  frequencies  are  precisely 
Rs  apart,  thereby  maintaining  orthogonality  between  the  N  sub-carriers  while 
maintaining  spectral  efficiency.  This  method  involves  using  the  Fourier  transform  [15]. 
Consider  the  resulting  signal  on  the  ith  sub-carrier, 

st{t)  =  Adj  ((t)cos (a^)  -  AdQ  .(t)sin (®,t)  (2.13) 
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where  .4  is  a  constant,  proportional  to  signal  amplitude,  coi  is  the  i"'  sub-carrier’s  center 
frequency,  and  dQ  ;(t)  and  d,  t(t)  are  the  respective  quadrature  and  in-phase  pulse 

amplitude  modulated  signals  of  the  i"1  sub-carrier  for  a  given  M  -  ary  modulation. 
Alternatively,  Equation  (2.13)  can  be  more  conveniently  written  as 

5;.(0  =  ARe{[j/,(0  +  ^e,(0]^}-  (2-14) 

We  can  further  simply  the  representation  of  st  ( t )  by  recognizing  from  Equation  (2.6)  that 


du(t)  +  jdQi(t)  =  dt(n)PT'  (t  -  nTs ),  (2.15) 


where/),  (t)  is  the  rectangular  pulse  function  defined  in  Equation  (2.7)  and 
di(n)  =  dI  +  jdQ  .  The  new  representation  of  st(t)  is 


st(t)  =  +Re 


+oo 

YJd,(n)Pr(t-nT,)elv 


(2.16) 


If  we  only  consider  one  symbol  duration  where  mTs  <t<  ( m  +  1)7) ,  and  recognizing  that 


coi  =  2^/',  the  result  is 

st(t)  =  AKQ{d,(m)ej27rft}.  (2.17) 

Considering  all  N  sub-carriers, 

s(t) = = (2-18) 

i=0  i=0 

Additionally,  if  the  ith  sub-carrier  frequency  is  defined  as  ft  =  f0  +  iAf  then  s(t)  can  be 
written  as 


s(t )  =  +Re 


N-l 


Yd^e^e^ 


(2.19) 
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As  previously  discussed,  orthogonality  is  achieved  when  Af  =  jr-  If  we  sample  s(t)at 


samples  t  =  m—,  where  m  =  0,1,...,  N  - 1,  the  result  is 
N 


(  T  ) 

m  — 

l  N) 

=  ATRe-j 

N- 1 


j2m—  j2rrf0m-f- 
7  N  a  N 


(2.20) 


The  Inverse  Discrete  Fourier  Transfonn  (IDFT)  is  defined 

1  '2jii m 

as  IDFT^d^m)]  =  — '^dj{m)e  N . 

N  i= o 

Consequently,  Equation  (2.20)  may  be  re-written  as 


(  T  \ 

m— 

v  Nj 


=  NAKe{lDFT[di(m)\e,2*f'f 


(2.21) 


or,  using  the  complex  envelope, 


f  t\ 
V  Nj 


=  NA  ■  IDFT  [dt  (m)\ . 


(2.22) 


Equation  (2.22)  shows  that  the  IDFT  algorithm  forms  the  kernel  for  transmitting 
the  serial  data  stream  into  N  Orthogonal  Frequency  Division  Multiplexed  data  streams. 
The  IDFT  algorithm  ensures  that  all  oscillator  frequencies  are  precisely  R,  apart,  thereby 
maintaining  orthogonality  between  the  N  signals.  The  IDFT  algorithm  is  practical  today 
due  to  faster  signal  processing  and  computational  efficiency  via  application  specific  IFFT 
integrated  circuits  [16]  and  reconfigurable  hardware  such  as  field  programmable  gate 
arrays  [8],  This  makes  the  IDFT  method  less  expensive  than  a  conventional  multi¬ 
oscillator  technique.  Additionally,  the  receiver  simply  implements  an  FFT  to  recover  the 
data  [17].  Another  advantage  is  that  the  receiver  can  separate  carriers  by  use  of  baseband 
signal  processing  techniques  vice  using  more  complex  bandpass  filtering  techniques  [18]. 
At  the  receiver,  the  sub-carrier  signal  can  then  be  integrated  over  the  symbol  period  to 
produce  the  desired  sub-carrier  decision  statistic.  All  other  sub-carrier  signals  are  filtered 
out  due  to  their  orthogonal  nature.  This  does  not  entirely  solve  the  synchronization 
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vulnerabilities  of  the  OFDM  signal.  For  example,  symbol  timing  becomes  critical  to 
ensure  that  only  the  appropriate  OFDM  symbols  are  included  in  the  receiver’s  FFT 
processing.  Failure  to  properly  time  the  FFT  processing  could  result  in  significant  ISI  by 
including  values  of  either  the  previous  or  follow-on  OFDM  symbol  as  we  shall  discuss  in 
section  III  [8]. 

In  summary,  this  chapter  addressed  the  necessity  to  use  multi-carrier  modulation 
to  achieve  high  data  rates  in  fading  channels.  OFDM  accomplishes  this  in  a  spectrally 
efficient  manner  using  the  IDFT  kernel  to  achieve  orthogonality  between  the  sub-carriers. 
In  the  next  chapter,  the  drawbacks  of  this  scheme,  with  respect  to  synchronization,  are 
addressed. 
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III.  SYNCHRONIZATION  VULNERABILITIES  OF  OFDM 


As  discussed  in  the  previous  section,  high  data  rates  may  be  achieved  with  per 
sub-carrier  flat  fading  and  minimal  ISI  by  multiplexing  and  transmitting  data  in  parallel 
via  sub-carriers.  The  most  efficient  way  to  do  this  is  through  OFDM,  using  an  IFFT-FFT 
processing  algorithm  [2],  This  ensured  orthogonality  can  be  maintained  between  the  sub¬ 
carriers  while  transmitting  the  sub-carriers  in  a  spectrally  efficient  manner. 

Unfortunately,  because  of  the  efficient  use  of  available  bandwidth  and  the 
simultaneous  transmission  of  N  sub-carriers,  OFDM  is  extremely  sensitive  to 
imperfections  in  frequency  synchronization.  Additionally,  the  multiplexing  of  symbols 
onto  multiple  sub-carriers,  in  a  multi-path  environment,  leads  to  a  system  that  is 
vulnerable  to  imperfections  in  timing  synchronization  as  wellEquation  Section  3.  [19] 

A.  INTER-CARRIER  INTERFERENCE 

Lack  of  frequency  synchronization  causes  inter-carrier  interference  (ICI).  OFDM 
is  vulnerable  to  ICI  because  of  the  simultaneous  and  spectrally  efficient  means  of 
transmitting  the  sub-carriers  as  seen  below. 

Frequency 


Figure  3  Depiction  of  ICI  -  from  Reference  [20] 
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In  Figure  3,  taken  from  [20],  local  oscillator  (LO)  drift  has  caused  the  receiver’s 
center  frequency  of  a  particular  sub-carrier  to  drift  from  the  maximum  value  (peak  of  the 
sine  function)  for  the  sub-carrier  of  interest.  As  the  diagram  shows,  orthogonality 
between  the  sub-carriers  has  been  destroyed.  In  the  diagram  above,  if  the  intention  is  to 
detect  the  middle  frequency  sub-carrier  there  will  be  interference  from  the  other  two  sub¬ 
carriers  (ICI).  Figure  4  is  provided  to  show  receiver  sampling  when  ICI  is  not  present. 
With  multiple  sub-carriers  present,  each  sub-carrier  is  sampled  by  the  receiver  at  its 
maximum  amplitude  (A fc  =0)  while  the  remaining  sub-carriers  have  zero  amplitude  at 

the  sampling  frequency,  preserving  orthogonality  and  maximizing  SNR.  This  is  much 
more  desirable  than  the  situation  depicted  in  Figure  3. 


Figure  4  Reception  without  ICI 


In  Figure  3,  a  single  sub-carrier  frequency  was  shown  to  drift.  The  IF  FT  /  FFT 
algorithm  previously  discussed  eliminates  the  need  for  N  individual  oscillators  for  N 
sub-carriers,  negating  the  problem  of  individual  sub-carrier  oscillator  drift.  This  scheme, 
however,  cannot  address  ICI  caused  by  the  possibility  of  a  frequency  offset  of  the 
receiver’s  local  oscillator,  phase  noise  inherent  in  the  transmission  system,  or  Doppler 
shift  between  receiver  and  transmitter. 
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1. 


Local  Oscillator  Offset 


Local  oscillator  (LO)  frequency  offset  is  a  frequency  difference  between  the  LO 
of  the  receiver  and  the  LO  of  the  transmitter.  This  will  result  in  a  frequency  offset  for 
every  sub-carrier  as  is  depicted  in  Figure  5. 


Figure  5  Sub-carrier  reception  with  ICI 


The  frequency  difference  between  the  receiver’s  and  transmitter’s  LO’s  is  the 
mathematical  equivalent  of  multiplying  the  post  IFFT  transmitter  complex  equivalent 
lowpass  signal  in  the  time  domain  by  a  complex  exponential  whose  frequency  is  equal  to 
the  frequency  offset  (difference).  This  is  the  same  as  shifting  the  received  signal  in 
frequency.  [21] 

Additionally,  close  examination  of  Figure  5  shows  that  there  is  a  reduction  in 
received  signal  strength  of  the  desired  sub-carrier  because  the  sub-carrier  is  no  longer 
sampled  at  the  peak  of  the  sine  squared  function.  This  reduction  in  signal  strength, 
combined  with  the  resulting  ICI,  results  in  lower  detection  rates  [8].  The  degradation  in 
SNR,  for  relatively  small  frequency  errors,  can  be  approximated  as  [13] 


D  =— l^— 

freq  3  In  10 


de¬ 


an 
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In  Equation  (3.1),  Dfreq  is  the  degradation  in  SNR,  in  dB,  due  to  frequency  offset, 

J7 

where  A f  is  the  frequency  offset  in  Hz,  — —  is  the  symbol  energy  to  noise  density  ratio  in 

Nq  - 

Watts/Watt,  and  Ts  is  the  symbol  period  in  seconds. 


Examination  of  Equation  (3.1)  shows  that  the  degradation  in  SNR  is  proportional 
to  the  square  of  the  frequency  offset.  Additionally,  recognizing  that  the  degradation, 

E 

Dp.  ,  is  proportional  to  — A  which  will  be  largest  when  reception  is  strong  and  higher 

N o 

order  M  -  ary  modulations  are  used.  Therefore  we  expect  degradation  in  SNR  due  to  the 
offset  to  be  most  pronounced  for  64  -  QAM  and  least  pronounced  for  BPSK. 


Figure  6  is  an  illustration  of  Equation  (3.1).  It  shows  the  sensitivity  of  various 
modulation  schemes  to  frequency  offset  at  a  given  bit  error  rate  (BER)  of  1CT6.  The 


— -required  for  BPSK/QPSK,  16-QAM,  and  64-QAM  modulation  schemes  to  achieve  a 


BER  of  1 0  6 ,  uncoded,  in  an  additive  white  Gaussian  noise  (AWGN)  channel  with  no 
fading  and  perfect  time,  phase  and  frequency  synchronization,  are  10.5  dB,  14.7  dB,  and 
18.8  dB  respectively  [22].  The  results  are  consistent  with  other  sources,  such  as  [8]  and 
[13]. 
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Figure  6  SNR  degradation  due  to  relative  frequency  error,  A/T 


Figure  6  indicates  that  for  a  degradation  of  no  more  than  0.5  dB,  the  frequency 
offset  cannot  exceed  0.85%  of  sub-carrier  spacing  for  64-QAM,  or  2.6565  kHz.  For  16- 
QAM,  frequency  offset  cannot  exceed  1.76%  of  sub-carrier  spacing,  or  5.500  kHz. 
Finally,  for  BPSK/QPSK,  frequency  offset  cannot  exceed  5.58%  of  sub-carrier  spacing, 
or  17.437  kHz.  Hence,  OFDM’s  sensitivity  to  frequency  offset,  as  depicted  in  Figure  6, 
is  apparent  and  confirmed  by  a  variety  of  sources  such  as  [8]  and  [9].  The  IEEE  802. 1  lg 
standard  stipulates  maximum  center  frequency  tolerance  of  ±25PPM  [5].  Under  the 
IEEE  standard,  the  maximum  possible  error  between  receiver  and  transmitter  occurs  if 
each  oscillator  experiences  the  maximum  deviation  of  25  PPM  with  opposite  sign, 
resulting  in  a50  PPM  worst  case  frequency  offset.  At  2.4  GHz,  this  translates  to  a 
maximum  allowable  frequency  error  of 
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A/;.  =  50  •  1 0  6  •  2.4  - 1 09  Hz  =  120  kHz. 


(3.2) 


From  Equation  (3.2)  and  Figure  6,  we  can  infer  that  despite  the  fact  that 
commercial  receivers’  and  transmitters’  respective  oscillators  will  meet  the  IEEE 
standard,  a  frequency  offset  of  60  kHz,  or  19.2  %  of  sub-carrier  spacing,  will  cause 
massive  degradation  of  over  200  dB  in  receiver  performance.  The  worst  case  scenario  of 
120  kHz,  or  38.4%  of  sub-carrier  spacing,  will  cause  a  degradation  of  over  1000  dB  in 
receiver  performance.  In  either  case,  802.1  lg  systems  must  correct  for  frequency  offset 
of  the  respective  local  oscillator  with  frequency  synchronization  or  performance  will  be 
severely  affected. 

2.  Phase  Noise 

The  additional  problem  of  phase  noise  is  related  to  the  quality  and  cost  of  the 
local  oscillators  (LO)  used.  The  ideal  oscillator  produces  only  the  required  frequency  and 
its  power  spectral  density  is  a  Dirac  delta  function  at  the  LO  frequency.  The  output  of  an 
actual  oscillator,  not  only  produces  the  desired  LO  frequency,  but  also  outputs  power 
concentrated  in  a  band  about  that  frequency.  These  additional  unwanted  frequency 
components  constitute  phase  noise  and  are  present  even  after  LO  frequency  offset 
correction  is  performed  by  the  receiver  [23].  Phase  noise  has  two  undesirable  effects  on 
receiver  perfonnance  -  ICI  and  constellation  rotation  [24], 

Recalling  from  Equation  (2.22),  if  the  complex  envelope  of  the  transmitted 

T 

OFDM  signal  for  a  given  symbol  is  sampled  at  t  =  m  —  ,  where  m  =  0,1,...,  A- 1,  then  the 
complex  envelope  of  the  transmitted  signal  becomes 

s(m)  =  NA^J^di(m)exv  j2ni^  j.  (3.3) 

In  Equation  (3.3),  the  exp (•)  term  is  used  for  clarity  vice  e']  when  the  structure 
of  the  exponents  would  cause  very  small  fonts  to  be  used.  If  we  sample  instead  at 
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t  =  mTs  vice  t  =  m  —  for  simplicity,  and  disregarding  constants,  Equation  (3.3)  can  be 


N 


re-written  as 


N- 1 


r(m)  =  XJ' exP 


/=0 


.  2  n  . 

l  — im 

v  N  j 


(3.4) 


The  received  complex  envelope  of  the  signal,  neglecting  additive  noise,  oscillator 
mismatch,  and  channel  fading,  is  only  affected  by  phase  noise  in  the  receiver  as  [25] 


r(m)  =  s(m)eJ 


(3.5) 


At  the  receiver,  the  signal  is  Orthogonal  Frequency  Division  de-multiplexed  including  the 
FFT  as 


y‘  =  N^{p)eX P[-JT¥ 


(3.6) 


which  can  be  expanded  to 


yk  =  77  S  eMP)‘ ?  ( P )  exP  [  -J  ~77  Pk 


N^o 


^  2n  k 
V  N  j 


(3.7) 


Using  Equation  (3.4),  Equation  (3.7)  can  be  re-written,  using  different  indices,  as 


N- 1  JV— 1 


^=77Ee  Z^exp 

p=0  r=0 


2jT 


r 

J~rrrP 

V  N 
( 


\ 

'  .2  n  ^ 

exp 

-J  —  Pk 

J 

l  N  J 

1  N- 1  JV— 1  /  O-rr 

=  77  £  S  dr  ^v[j—p{r-k) 


N^0 


(3.8) 


where  Jris  the  sub-carrier  symbol  defined  with  Equation  (2.15).  If  the  phase  offset, 
</>(p),  is  small,  it  may  be  approximated  as  eJ^m)  ~  1  +  j(f>{ p)  [25],  In  that  case,  yk becomes 


1  N~l  N~l  (  Itt  \  i  N~l  N~l  (  Itt  \ 

77Z<ZexP  j^p{r~k)  +  ^2H2>(^)ex:  p  j  —  p{r-k)\  (3-9) 

iv  r=Q  p=0  \  IV  J  !V  r= 0  p=0  V  tv  ) 


r= 0  /)=0 

The  first  tenn  on  the  right  hand  side  of  Equation  (3.9)  can  be  simplified  to 
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where  sk  is  the  total  error  caused  by  phase  noise  and  is  added  to  every  sub-carrier 
symbol.  It  represents  the  second  term  of  Equation  (3.10)  and  can  be  described  as  a 
complex  exponential  summed  up  for  every  sub-carrier.  As  a  result,  sk  can  be  re-written 
as 


7  N~l  N ~l  (  1  IT  \ 

^=Tf2XZ^)exP  J^p(r~k)  •  (3-12) 

tV  r=o  p= o  V  tV  y 

Hence,  each  sub-carrier  value  has  an  error  term  in  it  due  to  the  effects  of  phase  noise, 
which  in  turn  has  two  specific  effects.  The  first  effect  occurs  for  the  k,h  term  of^  r{r=k) , 

or 

7  k  N~X  f  2  7T  3 

£k,r(r=k)=  -^:Z^rZ^)exP  7  — =  jdk®.  (3.13) 

tv  r=k  p=0  v  tv  7 

where  the  average  phase  error, O,  is 
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(3-14) 


®  =  Yj0(p)- 


p= o 


Equation  (3.13)  shows  that  a  phase  error  of  jO  is  added  to  every  sub-carrier  value 
and  that  this  phase  error  is  proportional  to  the  received  sub-carrier  value  dk.  Each  sub¬ 
carrier,  regardless  of  the  value  of  dk  is  equally  affected  by  <f>.  Hence,  O  represents  the 

angular  magnitude  of  the  constellation’s  rotation  due  to  phase  noise  [25],  When  the 
rotation  of  the  constellation’s  data  points  exceeds  the  modulation  type’s  decision 
boundaries,  symbols  will  be  demodulated  incorrectly  [8],  Fortunately,  every  sub-carrier 
experiences  equal  constellation  rotation  and  this  can  be  leveraged  to  develop  effective 
techniques  to  correct  constellation  rotation  and  maintain  performance  [13]. 

The  second  effect  caused  by  the  error  tenn  sk  r  is  ICI.  The  ICI  is  caused  by  all  of 
the  remaining  sub-carriers  interfering  with  the  sub-carrier  being  demodulated.  In  this 
case,  a,  ...  becomes 

’  k,r[r^k) 


N- 1  N- 1 


.  2k 


£k,r(r*k)  =^2H2>0h>eXP|  j-pir-k ) 


r= 0 
r*k 


P= o 


N 


(3.15) 


Equation  (3.15)  shows  that  the  ICI-error  is  a  frequency  shifted  version  of  </>(p), 
the  frequency  shift  being  proportional  to  (r  -&) .  This  complex  number  is  added  to  each 
sub-carrier’s  useful  value,  reducing  its  amplitude  and  destroying  orthogonality  [25], 

To  minimize  degradations  caused  by  phase  noise,  wireless  engineers  will  model 
and  simulate  phase  noise  for  various  commercial  oscillators  to  triangulate  price, 
accuracy,  and  complexity  of  the  receiver’s  synchronization  technique  to  meet 
performance  and  cost  requirements.  To  properly  model  the  phase  noise  in  an  oscillator, 
the  phase  noise  is  treated  as  a  random  process  which  causes  the  oscillator  frequency  to 
deviate  from  the  desired  frequency.  The  strength  of  the  undesired  frequency  deviations 
can  be  shown  in  a  depiction  of  the  power  spectral  density  of  the  oscillator  (including 
phase  noise).  Ideally,  the  oscillator  would  only  generate  the  desired  LO  frequency,  but 
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actually  generates  random  signals  about  the  desired  frequency.  This  is  shown  in  the  PSD 
for  a  phase  locked  oscillator  depicted  in  Figure  7  below. 


Figure  7  Single-sided  LO  PSD  with  Phase  Noise  (/?  =  1  Hz) 


The  half  power  point,  or  Beta,  as  depicted  in  Figure  7,  can  be  used  to  estimate  the 
SNR  degradation  [9].  The  degradation  in  SNR,  caused  by  a  failure  to  correct  for  phase 
noise  of  the  LO,  is  approximated  as  [13] 


Dphase  6  In  1 0  '4^  ^  Nn  dB’ 


(3.16) 


where  ft  has  been  defined  above,  Ts  is  the  symbol  period  in  seconds,  and  is  the 
symbol  energy  to  noise  density  ratio  in  Watts/Watt.  Figure  8  depicts  the  degradation  in 
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SNR  due  to  phase  noise,  using  the  same  modulations  types  and  values  of  BER 
£ 

and  — -  that  we  used  in  the  previous  analysis  for  frequency  offset. 


Figure  8  SNR  degradation  caused  by  phase  noise 


The  horizontal  axis  of  Figure  8  is  normalized  to  the  802.1  lg  sub-carrier  spacing 
of  312.5  kHz  for  clarity.  If  an  arbitrary  SNR  degradation  of  .5  dB  is  acceptable,  then 
from  Figure  8,  the  one-sided  3  dB  oscillator  half-power  point  or  linewidth  (/?)to  sub¬ 
carrier  spacing  ratio  must  be  less  than  .0019  for  64-QAM.  This  translates  to  a  maximum 
acceptable  linewidth  of  594  Hz.  For  16-QAM  and  BPSK,  the  maximum  acceptable 
linewidth  is  relaxed  to  750  Hz  and  1.50  kHz  respectively. 

The  narrow  linewidth  limitations  required  to  avoid  noticeable  degradations  in 
SNR,  and  consequently  receiver  performance,  suggest  that  an  obvious  solution  would  be 
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to  purchase  the  highest  quality  LO  for  the  802. 1 1  g  receiver  and  transmitter.  The 
economic  feasibility  of  any  broadband  fixed  wireless  access  system,  however,  depends 
heavily  on  the  cost  of  individual  subscriber  units  and  their  components.  Unfortunately, 
the  cost  of  components,  including  oscillators,  is  directly  related  to  their  stability  and 
accuracy.  Hence,  any  commercial  OFDM  wireless  access  system  shall  have  to  meet 
acceptable  performance  standards  with  moderately  priced  components,  necessitating 
receiver  synchronization  techniques  to  correct  for  phase  noise  inherent  in  the  receiver’s 
LO.  [27] 


3.  Doppler  Shift 

The  final  phenomena  causing  ICI  is  Doppler  shift.  Doppler  shift  will  occur  when 
there  is  relative  motion  between  receiver  (Rx)  and  transmitter  (Tx).  The  Doppler  shift  is 
the  difference  between  the  received  signal  center  frequency  and  the  transmitted  signal 
center  frequency.  Considering  the  802. 1  lg  system,  the  Doppler  shift  for  a  particular  sub¬ 
carrier  is  given  as 

f«x  ~  fr,  =A fD=~  fs  cos «•  (3 • 1 7) 

c 

In  Equation  (3.17),  AfD  represents  the  deviation  from  sub-carrier  center  frequency 
due  to  Doppler  shift,  fs  is  the  transmitted  sub-carrier  center  frequency,  c  is  the  speed  of 

light,  and  v  cos «  is  the  relative  velocity  between  transmitter  and  receiver.  Equation 
(3.17)  only  considers  either  the  Rx  or  Tx  in  motion.  Most  scenarios  in  a  WLAN 
environment  would  involve  a  mobile  unit  (such  as  a  laptop  or  hand  held  device)  and  a 
station  unit  (such  as  a  wireless  access  point  or  router).  Hence  a  represents  the  angular 
difference  between  the  Rx’s  or  Tx’s  velocity  vector  and  a  straight  line  connecting  the 
two.  Figure  9  below  highlights  four  instances  of  Doppler  shift  in  which  the  Tx  is 
stationary  while  the  Rx  moves. 
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0°  <  ec  <  90° 


Case  1 


Case  4 


Figure  9  Doppler  shift  due  to  relative  motion  between  RX  and  TX 


If  v  ,  or  the  speed  of  the  Rx,  is  the  same  for  all  four  cases  of  Figure  9,  the 
maximum  Doppler  shift  occurs  when  a  =  0°,  as  depicted  in  Case  1  of  Figure  9.  The 
minimum  Doppler  shift  occurs  when  a  =  180°,  as  depicted  in  Case  4.  Cases  2  and  3 
highlight  how  a  can  vary  between  [(), tt]  and  thereby  change  the  value  of  A fD. 
Perpendicular  axes  are  shown  in  Figure  9  in  each  case  for  clarity. 

From  Equation  (3.17),  Doppler  shift  has  a  greater  effect  at  high  frequencies. 
Consequently,  when  considering  multiple  sub-carriers,  higher  frequency  sub-carriers  will 
experience  a  greater  frequency  shift  due  to  Doppler.  Equation  (3.17)  may  be  re-written 
for  multiple  sub-carriers  as 
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(3.18) 


¥dj  =-fcosa, 
c 

where  i  =  \...N  sub-carriers.  The  percentage  change  in  frequency  due  to  Doppler  shift, 
however,  will  be  the  same  for  all  sub-carriers,  when  each  sub-carrier’s  Doppler  shifted 
frequency  is  normalized  to  its  respective  transmitted  sub-carrier  center  frequency  [28].  In 
this  case  (for  simplicity,  let  a  =  0), 


A fpj  v 

fi  c 


(3.19) 


where  <^is  the  percentage  change  in  frequency  due  to  Doppler  shift.  From  Equation 
(3.19),  the  percentage  change  in  frequency  for  all  sub-carriers  is  only  affected  by  the 
relative  motion  between  the  Tx  and  Rx.  Further,  if  a  particular  sub-carrier  is  transmitted 
at  frequency  fs,  then  the  received  frequency  will  be  /'  {  \+  g)  due  to  Doppler.  For 
multiple  sub-carriers,  the  received  Doppler  shifted  frequencies  becomes 

(1  +  f  )(fc  +  htf. )  =  (1  +  f  )fc  +  (1  +  4 )  hAf„  (3.20) 


where  fc  is  the  carrier  frequency,  A fs  is  the  sub-carrier  frequency  spacing  and 
f-N 


h  = 


Unlike  Equation  (3.18),  Equation  (3.20)  describes  Doppler  shifted 


frequencies  in  which  a  carrier  frequency  is  bracketed  above  and  below  by  sul> 


carriers.  From  Equation  (3.20),  it  can  be  seen  that  the  carrier  frequency,  sub-carrier 
frequencies,  and  sub-carrier  spacing  are  affected  by  the  same  percentage  For  an 
OFDM  signal,  recall  the  following  time-domain,  transmitted  OFDM  signal  from 
Equation  (2.9): 


N- 1 

s(t )  =  ^a;.(t)cos[(<yf  +  icoA)t  +  0i(t)\ 

i= 0 


(3.21) 


where  coc  is  the  carrier  frequency,  coA  is  the  sub-carrier  spacing,  a:  is  dependent  upon 
signal  amplitude,  6j  varies  with  the  in-phase  and  quadrature  modulation  values  and  for 
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clarity,  sub-carriers  i  =  0 . . .  N  - 1  are  considered.  When  a  Doppler  shift  is  included,  the 
received  signal  becomes 

N- 1 

r(t)  =  Y  a,  (0  cos  [(1  +  4)  (coc  +  icoA  )t  +  0i  (f)],  (3 .22) 

(=0 

where  the  effect  of  the  Doppler  shift  on  0j  (t)has  been  ignored.  Substituting  Equation 
(3.20)  and  recalling  Equation  (2.8),  Equation  (3.22)  may  be  re-written  as 

N- 1 

r(t)  =  Y{al  (t)  cos  [(1  +  %)icoAt  +  0j  (t)]  cos  [( 1  +  £) 

i=o  (3.23) 

-  at  ( t )  sin  [(1  +  %)icoAt  +  6*  (t)]  sin  [(l  +  £)  coci^  | 

If  the  factor  cos[(l  +  £,  )icoAt  +  m]  can  be  considered  as  the  envelope  of  the 
carrier  cos[(  1  +  g)cort^  and  similarly  for  the  second  term  in  Equation  (3.23),  the 

expression  shows  that  both  carrier  frequency  and  envelope  frequency  are  affected  by  the 
same  Doppler  percentage  [28] 

Finally,  clock  timing  is  derived  from  the  OFDM  symbol  rate  f?s.  Recall  from 
Equation  (2.12),  that  to  maintain  orthogonality,  sub-carrier  spacing  must  be  an  integer 
multiple  of  Rs.  When  the  sub-carrier  spacing  is  shifted  in  frequency  due  to  some 
Doppler  percentage  orthogonality  between  the  sub-carriers  is  destroyed,  creating  ICI. 
[28] 

Consequently,  Doppler  shift  affects  the  frequencies  of  the  carrier,  sub-carriers, 
envelope,  and  symbol  timing  by  the  same  Doppler  percentage  thereby  creating  ICI. 
The  degradation  in  performance  caused  by  Doppler  shift  has  been  estimated  by  many 
sources.  It  has  been  shown  that  the  fractional  degradation  in  received  power  due  to  ICI  in 
an  OFDM  system  with  flat  fading  in  any  Doppler  spectra  is  universally  upper  bounded  by 
[29] 

(3-24) 
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where  PICI,  as  a  function  of  A/',  and  Ts,  is  the  unitless  degradation  of  received  signal 

power  caused  by  Doppler  shift  induced  inter-carrier  interference.  In  Equation  (3.24), 
received  signal  power  is  normalized  to  1  watt  for  clarity.  Hence,  a  computed  magnitude 
of  .1  for  PICI  translates  to  a  loss  of  10%  in  received  signal  power.  Although  Equation 
(3.24)  reveals  that  the  PICI  is  a  function  of  A fD  and  Ts,  symbol  period  in  802.1  lg 
systems  is  fixed  at  3.2  //sec  [4].  Additionally,  commercial  WLANs  are  designed  to  be 
used  inside  buildings,  or  small  outside  areas,  with  fixed  access  points.  Mobile 
subscribers  are  usually  stationary  or  moving  at  pedestrian  speeds.  A  brisk  walking  speed 
of  3  m/s  yields  a  worst  case  Doppler  shift  of 

_  3m/ sec — 2.4-109  =  24  Hz.  (3.25) 

D  3 -108  m/sec 

The  resulting  maximum  fractional  degradation  in  received  power  due  to  ICI, 
using  Equation  (3.24),  turns  out  to  be 

PICI  <^(2/r-24Hz-3.2  //sec)2  =  .0000000194,  (3.26) 

or  .00000194  %  loss  in  received  signal  power.  This  negligible  effect  is  further  confirmed 
by  visual  inspection  of  Figure  6  as  24  Hz  represents  only  .007  %  of  sub-carrier  spacing. 
At  that  amount  of  frequency  offset,  there  is  little  effect  on  overall  system  performance. 
Therefore,  ICI  caused  by  Doppler  shift  will  be  much  less  than  ICI  caused  by  phase  noise 
or  LO  offset  in  commercial  systems  [29]. 

4.  Summary 

In  summary,  inter-carrier  interference  (ICI)  is  caused  by  three  phenomena  -  local 
oscillator  offset,  phase  noise  and  Doppler  shift.  While  Doppler  shift,  under  normal 
conditions,  was  shown  to  have  negligible  effect  on  system  performance,  LO  offset  and 
phase  noise  exposed  the  vulnerabilities  of  OFDM  systems.  This  is  due  to  the  spectrally 
efficient  manner  that  the  sub-carriers  are  transmitted  over  the  wireless  channel.  Offset 
sampling  by  the  OFDM  receiver  produced  a  reduction  in  signal  strength  and  ICI  from  the 
remaining  sub-carriers.  In  the  case  of  LO  offset,  it  was  shown  that  even  small  deviations 
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in  frequency  offset  caused  noticeable  degradations  in  perfonnance,  particularly  for  higher 
order  modulations  (e.g.  64-QAM).  It  was  also  shown  that  this  degradation  in 
performance  could  occur  despite  having  met  the  IEEE  802. 1 1  g  oscillator  standard. 
Consequently,  frequency  synchronization  is  an  essential  task  for  the  OFDM  receiver. 
Phase  noise  was  shown  to  produce  noticeable  degradation  in  system  performance  as  well. 
The  degradation  from  phase  noise  was  directly  related  to  the  3  dB  linewidth  of  the 
receiver’s  oscillator.  The  obvious  solution  of  purchasing  high  quality  oscillators  was 
discarded  due  to  their  high  cost  necessitating  phase  tracking  by  the  OFDM  receiver. 
Frequency  [8]  and  phase  synchronization  are  essential  tasks  for  commercial  OFDM 
receivers  today  and  various  techniques  have  been  proposed  [24].  These  techniques  shall 
be  explored  in  Section  IV. 

B.  INTER-SYMBOL  INTERFERENCE 

The  previous  section  highlighted  the  vulnerabilities  of  OFDM  to  ICI.  From  the 
discussion  in  chapter  II,  and  from  Equations  (2.3)  and  (2.4),  it  is  shown  that  OFDM  is 
particularly  resistant  to  the  effects  of  inter-symbol  interference  (ISI).  Individual  sub¬ 
carrier  symbol  duration  is  longer  than  if  a  single  high  data  rate  carrier  was  used, 
mitigating  multi-path  effects. 

Regrettably,  ISI  is  not  entirely  eliminated.  Channel  impulse  response  (delay 
spread  of  the  channel)  could  affect  follow-on  OFDM  symbols  [8],  In  addition,  symbol 
timing  offset  can  introduce  ISI  [30], 

1.  Long  Delay  Spreads  and  ISI 

Long  delay  spreads  for  certain  channels  could  induce  ISI  for  the  802.1  lg  system. 

Per  the  IEEE  802.1  lg  standard,  the  sub-carrier  symbol  duration  is  3.2  //sec. 
lengthy  symbol  duration  will  generally  ensure  flat  fading  within  a  particular  sub-carrier, 
precluding  ISI.  Recall  that  Equation  (2.3)  stipulates  a  rule  of  thumb  for  flat  fading. 
Inspection  of  Table  1  shows  that  for  all  but  the  longest  channel  delay  spreads,  the  flat 
fading  rule  of  thumb  is  adhered  to  for  sub-carrier  duration.  This  is  further  confirmed  by 
the  fact  that,  according  to  Table  1,  a  typical  office  building  has  a  median  RMS  delay 
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spread  of  25  //sec.  This  corresponds  to  only  0.7%  greater  path  length  than  line  of  sight 
for  802.11g’s  sub-carrier  symbol  duration.  However,  for  channels  with  extremely  long 
RMS  delay  spreads,  ISI  could  manifest  itself  when  OFDM  symbols  overlap  in  time.  [31] 

The  following  table  lists  delay  spreads  for  various  channels  between  900  MHz 
and  1.9  GHz  but  can  be  extrapolated  to  higher  frequencies  as  an  upper  bound. 


Table  1  Typical  Values  of  RMS  Delay  Spread  -  from  Reference  [1] 


Environment 

RMS  Delay  spread  (cr^) 

Notes 

Urban 

1300  //sec  average 

600  //sec  std.  dev. 

3500  //sec  maximum 

New  York  City 

Urban 

10-25  //sec 

Worst  case,  San  Francisco,  CA 

Suburban 

200-310  //sec 

Averaged  typical  case 

Indoor 

10-50  //sec 

25  //sec  median 

Office  building 

Indoor 

270  //sec  maximum 

Office  building 

Indoor 

70-94  //sec  average 

1470  //sec  maximum 

Three  San  Francisco  buildings 

From  Table  1,  we  conclude  that  802.1  lg  systems  will  be  resistant  to  ISI  except  in 
most  channels  with  particularly  long  delay  spreads.  To  ensure  ISI  is  eliminated  between 
follow-on  OFDM  symbols,  a  guard  period  is  inserted  between  the  symbols.  The  IEEE 
standard  specifies  the  duration  of  the  guard  interval  -.8  //sec  (7f).  Figure  10  shows 
how  the  guard  interval  provides  ISI  protection  for  follow-on  OFDM  symbols. 
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Figure  10  Guard  interval  and  channel  response  -  from  Reference  [20] 


ISI  from  each  OFDM  symbol  is  depicted  above  by  a  generic  channel  impulse 
response  (CIR),  decreasing  in  intensity  with  time.  Note  that  once  a  properly  chosen 
guard  interval  has  been  placed  before  each  follow-on  OFDM  symbol,  the  effects  of  ISI 
are  essentially  eliminated.  As  Table  1  depicts,  ISI  could  appear  in  the  symbol  estimate  if 
a  channel  has  an  extremely  long  delay  spread,  but  a  properly  designed  system  will  have 
minimum  ISI  due  to  the  weak  values  of  the  channel  response  far  into  the  guard  interval. 
Critically,  the  length  of  the  guard  interval  must  be  greater  than  the  maximum  acceptable 
delay  spread  value  [13].  Total  symbol  duration  (Ttotal)  will  now  be: 

Tto,al=Tgl+Ts  (3-27) 

where  Ts  is  the  original  OFDM  symbol  duration  of  3.2  //sec  and  T .  is  the  guard  interval 

duration.  In  theory,  the  guard  interval  duration  could  be  made  to  vary  with  the  channel 
impulse  response  (time  delay  spread  of  the  channel),  however,  per  the  802.1  la  standard, 
from  which  the  802.1  lg  standard  adopts,  guard  interval  duration  is  fixed  at  20%  of  Ttotal . 

[4]  The  drawback  in  accommodating  the  guard  interval  is  overhead  with  respect  to  data 
transmission  rate,  but  well  worth  the  price  to  mitigate  ISI. 

2.  Symbol  Timing  Offset  and  ISI 

Symbol  timing  offset  can  also  introduce  ISI  within  an  802.1  lg  system.  As  the 

receiver  begins  its  demodulation  process,  the  FFT  window  length  shall  equal  the  length 
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of  the  OFDM  symbol.  Ideally,  the  FFT  window  could  start  with  the  first  data  sample  of  a 
particular  symbol.  The  window  would  end  with  the  last  data  sample  of  the  same  symbol. 
Figure  1 1  depicts  this. 


Ideal  FFT  window 

GI3 

OFDM  symbol  2 

GI2 

OFDM  symbol  1 

<• -  time 

Figure  11  Ideal  FFT  window  timing 

Note  that  in  Figure  11,  when  demodulating  OFDM  symbol  2,  no  FFT  samples  are 
taken  from  Guard  Interval  2.  This  provides  the  maximum  ISI  protection,  as  samples  of 
the  guard  interval  may  be  affected  by  multi-path  effects  from  OFDM  symbol  1 . 

In  practice,  it  is  impossible  to  fix  the  symbol  timing  point  at  exactly  the  first 
sample  of  an  OFDM  symbol.  There  is,  invariably,  some  deviation  of  symbol  timing 
around  a  mean  value.  [8] 

In  Figure  11,  early  FFT  window  timing  means  including  samples  from  guard 
interval  2,  while  omitting  some  of  the  last  samples  from  OFDM  symbol  2.  Late  FFT 
window  sampling  would  mean  including  samples  from  guard  interval  3  and  omitting  the 
initial  samples  of  symbol  2. 

The  sampled  guard  interval  could  have  been  designed  to  contain  no  signal  data 
but  then,  in  a  multi-path  environment,  orthogonality  would  have  been  lost  between  the 
sub-carriers,  resulting  in  inter-carrier  interference  (ICI).  Instead,  the  IEEE  standard 
mandates  that  the  OFDM  symbol  is  cyclically  extended  into  the  guard  interval.  Flence 
the  guard  interval,  serves  two  functions  -  protect  the  OFDM  symbols  against  ISI  and 
mitigate  ICI  in  a  multi-path  environment  by  maintaining  orthogonality  between  sub¬ 
carriers.  [13] 
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From  Equation  (2.12),  orthogonality  between  sub-carriers,  is  achieved  when  the 
frequencies  of  the  sub-carriers  are  integer  multiples  of  each  other.  To  maximize  spectral 
efficiency,  the  first  sub-carrier  is  set  to  integer  1,  yielding  a  sub-carrier  spacing  of  1/7) . 
Neglecting  the  carrier  frequency,  additional  sub-carriers  would  be  at  frequency 
i/Ts,  where  i  =  1 ...  N  sub-carriers.  Figure  12  shows  an  OFDM  receiver’s  FFT  interval 

with  three  sub-carriers.  Although,  these  sub-carriers  would  be  summed  together  in  an 
OFDM  signal,  they  are  depicted  separately  for  clarity.  The  number  of  cycles  in  the  FFT 
interval  of  sub-carriers  2  and  3  are  integer  multiples  of  the  first  sub-carrier,  achieving 
orthogonality.  Furthermore,  a  guard  period  exists  prior  to  the  FFT  interval.  In  Figure  12, 
the  guard  interval  is  filled  with  zero-valued  samples. 


Figure  12  3  sub-carriers  in  FFT  interval  -  no  ICI  present 
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If  a  sub-carrier’s  arrival  at  the  receiver  is  delayed  due  to  multi-path,  there  will  no 
longer  be  integer  multiples  of  sub-carrier  cycles  in  the  FFT  interval,  destroying 
orthogonality  between  the  sub-carriers.  This  is  depicted  in  Figure  13 


I  Guard  I 
I  interval 


FFT  interval 


► 


Figure  13  3  sub-carriers  with  ICI  caused  by  sub-carrier  delay 


In  Figure  13,  sub-carrier  3  is  delayed  and  no  longer  has  an  integer  multiple 
number  of  cycles  within  the  FFT  interval  relative  to  the  other  sub-carriers.  Hence,  sub¬ 
carrier  3  will  cause  cross-talk,  or  create  ICI,  with  sub-carriers  1  and  2  and  vice  versa. 
The  solution  is  to  cyclically  extend,  or  copy  samples  from  sub-carrier  3  and  append  it  to 
the  beginning  of  sub-carrier  3.  The  result  is  shown  in  Figure  14  where  delayed  sub¬ 
carrier  3  now  has  an  integer  number  of  cycles  within  the  FFT  interval.  [13] 
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Figure  14  3  sub-carriers  with  CP  added  to  delayed  sub-carrier 


Because  the  guard  interval  actually  consists  of  samples  of  the  follow  on  OFDM 
symbol  it  is  more  correctly  known  as  the  cyclic  prefix  (CP).  In  Figure  1 1  above,  guard 
interval  2  actually  consists  of  a  repetition  of  a  portion  of  OFDM  symbol  2’s  sample 
values.  This  raises  the  potential  for  ISI  contamination  from  OFDM  symbol  1  but  as 
mentioned  previously,  channel  response  values  should  be  weak  far  into  the  guard 
interval.  Samples  from  the  guard  interval,  now  cyclic  prefix,  contain  the  appropriate 
values  for  correct  demodulation.  Consider  the  case  of  late  FFT  window  timing,  as 
depicted  in  Figure  15  (b).  The  FFT  window  includes  samples  from  CP  3  and  omits  the 
beginning  of  OFDM  symbol  2.  Cyclic  Prefix  3,  however,  consists  of  samples  from 
follow-on  OFDM  symbol  3  (not  depicted).  The  result  from  this  late  symbol  timing  is 
severe  ISI  from  CP  3,  which  contains  samples  from  the  follow  on  symbol.  [8] 
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OFDM  symbol  1 

(b)  Unacceptable  (late)  FFT  window  -  severe  ISI 

Figure  15  Early  vs.  Late  FFT  window  timing  and  ISI 


3.  Summary 

The  solution  to  avoiding  ISI  caused  by  late  symbol  timing  offset  is  to  ensure  that 
the  FFT  window  starts  early.  This  means  including  samples  from  the  appropriate  cyclic 
prefix.  Orthogonality  is  still  maintained  because  the  CP  contains  the  appropriate  samples 
values  for  correct  demodulation.  If,  as  Table  1  suggests,  the  length  of  the  CP  is  longer 
than  the  CIR,  ISI  is  significantly  reduced  allowing  the  OFDM  receiver  to  take  advantage 
of  the  circular  aspects  of  the  FFT.  With  ISI  mitigated,  the  linear  convolution  between  the 
symbol  and  channel  becomes  circular.  In  other  words,  the  FFT  window  should  begin 
well  inside  the  CP  (when  the  values  of  the  channel  response  due  to  the  prior  symbol  are 
weak)  and  end  prior  to  the  CP  of  the  follow  on  symbol  to  avoid  ISI. [8] 

In  summary,  this  chapter  showed  clearly  that  802.1  lg  receivers  must  implement 
some  form  of  symbol  synchronization.  Failure  to  properly  synchronize  in  time  and 
frequency  resulted  in  ISI  and  ICI.  The  next  chapter  will  address  various  synchronization 
algorithms  to  overcome  these  problems. 
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IV.  SYNCHRONIZATION  OF  THE  OFDM  SIGNAL 


A.  INTRODUCTION 

As  shown  in  Chapter  III,  failure  of  the  receiver  to  properly  synchronize  in  time 
and  frequency  can  seriously  impair  the  performance  of  an  OFDM  system.  To  ensure  ISI- 
free  detection,  precise  timing  synchronization  must  occur  to  detennine  the  symbol 
boundaries  within  an  OFDM  packet.  To  ensure  sub-carrier  orthogonality  and,  hence, 
ICI-free  detection  and  performance,  frequency  synchronization  must  occur  [32].  Prior  to 
these  tasks,  the  packet  must  be  detected  by  the  receiver.  This  is  the  first  synchronization 
task  performed  that  future  synchronization  tasks  are  dependent  upon  [8].  Consequently, 
to  properly  detect  and  demodulate  the  OFDM  signal,  the  receiver  and  transmitter  must  be 
properly  synchronized.  The  algorithms  that  make  up  the  OFDM  synchronization  process 
represent  a  major  design  problem  that  has  to  be  solved  to  field  a  successful  product  [8]. 

1.  OFDM  Synchronization  and  Single  Carrier  Synchronization 

Synchronization  in  the  OFDM  multi-carrier  system  can  be  very  different  from 
synchronization  in  a  single  carrier  system.  Specifically,  the  nature  of  the  OFDM  signal 
wavefonn  makes  the  majority  of  the  traditional  single  carrier  synchronization  algorithms 
unusable  [8].  While  the  parallel  transmission  of  N  sub-carriers  results  in  longer  symbol 
duration,  and  consequently  less  sensitivity  to  timing  offset  than  in  a  single  carrier  system 
[31],  OFDM  is  very  sensitive  to  frequency  offset,  as  shown  in  Chapter  III.  Single  carrier 
systems  may  experience  a  similar  degradation  in  SNR  from  frequency  offset,  but  do  not, 
in  general,  suffer  the  loss  of  orthogonality  and  the  resulting  interference  from 
neighboring  sub-carriers.  Also,  many  of  the  OFDM  synchronization  functions  may  be 
performed  in  the  time  or  frequency  domain,  a  level  of  flexibility  that  is  not  found  in 
single  carrier  systems  [8],  Ultimately,  the  variety  of  synchronization  schemes  available 
within  OFDM  presents  a  trade-off  between  performance  and  computational  complexity. 
[8] 
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2.  Synchronization  for  WLANs  and  Broadcast  Systems 

OFDM  is  used  for  both  broadcast  communication  systems  and  ad  hoc  wireless 
networks.  The  nature  of  the  packet  switched  WLAN  will  determine  the  appropriate 
synchronization  algorithms.  Broadcast  systems,  such  as  Digital  Audio  Broadcasting  or 
Digital  Terrestrial  Television  Broadcast,  can  employ  techniques  that  may  require  a 
relatively  long  time  to  acquire  and  then  track  the  signal.  This  luxury  is  not  permitted  for 
WLANs.  The  packet  switched  WLAN  can  expect  random,  short  bursts  of  data  [8].  In 
802.1  lg  the  maximum  permitted  length  of  the  Medium  Access  Control  Protocol  Data 
Unit  (MPDU)  is  4095  bytes  [3].  The  preamble  and  header  fields  comprise20  //sec  [4], 
resulting  in  a  627  //sec  signal  duration  at  the  maximum  data  rate  of  54  Mbps  [12].  The 
header  is  required  for  successful  demodulation  of  the  OFDM  data  symbols,  leaving  only 
the  preamble  duration  of  16  //sec  for  synchronization  [12].  Therefore,  synchronization 
must  occur  very  rapidly  with  minimal  overhead  to  maintain  high  throughput  [12]. 
Consequently,  the  IEEE  802.1  lg  preamble  has  been  designed  to  facilitate 
synchronization  [8].  This  fundamentally  differs  from  broadcast  systems,  in  which  the 
signal  is  always  being  transmitted. 

3.  Synchronization  Methodology 

No  explicit  synchronization  technique  is  mandated  by  the  IEEE  802.1  lg  standard. 
Many  different  types  of  synchronization  algorithms  have  been  proposed  to  maximize 
performance  [33]. 

For  the  OFDM  WLAN,  synchronization  is  usually  accomplished  in  two  phases: 
acquisition  and  tracking.  A  good  acquisition  technique  should  be  robust  with  as  wide  an 
acquisition  range  as  possible  and  should  obtain  a  feasible  timing  or  frequency  estimate 
quickly.  Furthermore,  the  better  the  estimate,  the  better  the  overall  bit  error  performance. 
[33] 

Synchronization  algorithms  can  be  divided  into  three  types: 

•  Non  data-aided  synchronization  that  depends  on  the  spectral 
characteristics  of  the  received  signal 
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•  Cyclic  prefix  based  synchronization  that  exploits  the  redundancy  present 
in  the  OFDM  packet. 

•  Data-aided  synchronization  that  exploits  training  data  or  pilot  symbols 
within  the  OFDM  packet.  [12] 

a.  Non  Data-Aided  Synchronization 

Non  data-aided  synchronization  does  not  need  training  symbols  or  pilot 
tones,  reducing  overhead  and  increasing  data  throughput.  Because  training  symbols  are 
not  required,  non  data-aided  synchronization  appears  well  suited  for  broadcast  type 
systems  because  the  receiver  does  not  have  to  wait  for  the  training  symbols  to  appear  at 
the  beginning  of  the  data  stream.  Unfortunately,  the  performance  of  non  data-aided 
synchronization  algorithms  is  substandard  in  the  presence  of  multi-path  interference  and 
their  accuracy  is  generally  not  sufficient  to  guarantee  orthogonality  during  OFDM  packet 
based  transmission  [12]. 

b.  Cyclic  Prefix  Based  Synchronization 

Cyclic  prefix  based  synchronization  exploits  the  redundancy  within  the 
OFDM  data  symbol.  As  described  in  chapter  III,  the  cyclic  prefix  has  been  inserted  to 
protect  follow  on  data  symbols  from  ISI  but  also  contains  copies  of  the  data  samples  to 
preserve  sub-carrier  orthogonality.  Since  the  cyclic  prefix  is  a  repetition  of  some  fraction 
of  the  OFDM  symbol,  analyzing  the  change  in  the  repeated  data  can  be  used  to  maintain 
synchronization.  Generally,  the  complexity  of  cyclic  prefix  based  synchronization  is  less 
than  data-aided  or  non  data-aided  synchronization,  allowing  for  lower  cost  hardware  [12]. 
There  is  no  additional  overhead  to  be  transmitted  because  the  cyclic  prefix  is  already 
mandated  by  the  IEEE  standard.  If  the  OFDM  receiver  were  to  simply  use  this  low 
complexity  technique,  one  immediate  problem  would  be  packet  detection.  This  initial 
synchronization  task  would  be  dependent  upon  a  cyclic  prefix  at  the  start  of  the  OFDM 
packet.  [4]  and  [5]  make  no  specific  provision  for  a  CP  at  the  beginning  of  the  preamble. 
This  is  confirmed  by  Figure  17  which  depicts  the  802. 1 1  g  preamble  with  no  CP  for  the 
start  of  the  OFDM  packet.  Furthermore,  unless  cyclic  prefix  synchronization  is 
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combined  with  other  synchronization  techniques,  it  will  take  several  guard  intervals 
before  synchronization  is  obtained,  something  that  is  not  practical  for  packet-based 
OFDM  [13].  Because  of  this,  cyclic  prefix  based  synchronization  is  best  suited  for 
broadcast  systems  [12]. 

c.  Data-Aided  Synchronization 

Data-aided  synchronization  algorithms  depend  on  special  synchronization 
symbols  that  must  be  added  to  the  data  packet  and  /  or  pilot  tones  embedded  within  the 
OFDM  sub-carriers.  The  use  of  data-aided  algorithms  to  achieve  synchronization 
therefore,  by  definition,  reduces  data  throughput  and  requires  increased  overhead.  On  the 
other  hand,  data-aided  synchronization  provides  the  widest  acquisition  range,  best 
frequency  tracking  and  can  meet  the  requirement  of  the  packet-based  802.1  lg  system  for 
synchronization  within  16 //sec  of  packet  arrival  [12].  Because  the  802. 1  lg  IEEE 
standard  mandates  the  use  of  synchronization  symbols  in  the  preamble  and  designates 
certain  sub-carriers  as  pilots,  the  overhead  associated  with  data-aided  synchronization  is 
already  borne  by  the  packet  in  order  to  meet  the  IEEE  standard.  A  general 
communications  engineering  design  principle  is  that  the  receiver  should  use  all  the 
available  information  to  its  advantage  [8].  As  a  result,  we  can  expect  many 
synchronization  algorithms  to  exploit  these  parts  of  the  OFDM  packet  structure  because 
of  the  requirement  to  achieve  fast  and  accurate  synchronization  while  minimizing 
additional  overhead.  Viable  data-aided  synchronization  techniques,  because  they  seem  to 
be  the  likely  candidates  to  improve  perfonnance  of  the  802.1  lg  system,  will  be  the  focus 
of  the  remainder  of  this  chapter  and  those  parts  of  the  OFDM  packet  that  data-aided 
synchronization  techniques  exploit  will  be  addressed  in  the  following  section. 

B.  IEEE  802.11G  PACKET  STRUCTURE  AND  SYNCHRONIZATION 

1.  Physical  Layer  Convergence  Protocol  Receive  Procedure 

The  IEEE  802. 1  lg  standard  promulgates  an  OFDM  packet  structure  that  allows 
for  synchronization  in  time  and  frequency.  The  OFDM  packet  is  sent  from  the 
transmitter  to  the  receiver  in  accordance  with  the  802. 1  lg  physical  layer  convergence 
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function  (PLCF),  which  is  supported  by  a  specified  physical  layer  convergence  protocol 
(PLCP).  The  PLCP  defines  the  method  of  mapping  802. 1 1  g  medium  access  control 
(MAC)  sub-layer  protocol  data  units  (MPDUs)  into  a  frame  format  suitable  for  sending 
data  between  transmitter  and  receiver  given  the  characteristics  of  a  particular  wireless 
channel.  The  characteristics  of  the  wireless  channel,  and  the  method  of  transmitting  data, 
form  the  physical  medium  dependent  (PMD)  portion  of  the  PLCF.  [4] 

The  OFDM  PLCP  receive  procedure  is  depicted  below. 
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Figure  16  PLCP  receive  procedure  -  from  Reference  [4] 


Figure  16  shows  that  the  PLCP  includes  a  12  symbol  preamble,  a  signal  field  that 

is  convolutionally  coded  and  modulated  using  BPSK,  and  data.  Information  about  the 

data  rate  and  length  of  the  data  field  is  included  in  the  signal  field,  which,  coupled  with 

the  service  field,  comprises  the  header  of  the  OFDM  packet.  The  coded  data,  consisting 

of  OFDM  symbols,  pad,  and  tail  bits,  comprise  the  PLCP  data  service  unit  (PDSU).  The 

OFDM  receiver  uses  the  preamble  to  generate  the  Received  Signal  Strength  Indication 

(RSSI).  In  a  wireless  network,  such  as  802. 1  lg,  receivers  will  not  know,  a  priori,  when  a 

packet  will  be  received.  A  strong  value  for  RSSI  indicates  the  successful  detection  of  an 
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OFDM  packet,  a  necessary  requirement  before  additional  synchronization  tasks  can  be 
attempted.  Figure  17  closely  depicts  the  preamble  and  header,  including  cyclic  prefixes, 
and  highlights  the  synchronization  tasks  that  must  be  perfonned  prior  to  successful 
demodulation  of  the  transmitted  data  symbols. 


Figure  17  OFDM  preamble  and  header  structure  -  from  Reference  [4] 


The  preamble  consists  of  10  repetitions  of  the  same  symbol,  labeled  in  the 
diagram  as  tx...tw, 2  additional  symbols  repeated,  labeled  as  7j  and  T2,  and  GI 2,  which 

serves  as  a  cyclic  prefix  for  symbol  7],  Underneath  the  preamble  symbols,  the  text  in 
Figure  17  delineates  the  various  synchronization  tasks  that  can  be  performed  to  train  the 
receiver.  Due  to  the  short  duration  of  0.8  //sec  each  for  symbols  tl...tw,  these  are 
known  as  short  training  symbols  while  7]  and  7],  each  with  symbol  duration  3.2  //sec, 
are  known  as  the  long  training  symbols.  [4] 

It  is  these  repeated  symbols  and  the  embedded  pilot  tones  within  the  sub-carriers 
that  many  of  the  data-aided  802. 1  lg  synchronization  algorithms  exploit  to  quickly 
synchronize  to  the  WLAN  signal. 

2.  Short  Training  Symbols 

The  10  short  training  symbols  mentioned  in  the  previous  section  form  the 
beginning  of  the  OFDM  packet.  Each  of  the  symbols  is  16  samples  long  and  is  used  for 
packet  detection,  automatic  gain  control  (AGO)  and  diversity  selection.  The  symbols  can 
also  be  exploited  for  coarse  frequency  synchronization  and  symbol  timing.  [12] 
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From  Figure  17,  each  symbol  has  duration  of  0.8  //sec,  resulting  in  a  total 
duration  of  8.0  //sec  for  the  short  training  symbols.  The  symbols  are  generated  by  the 
802.1  lg  transmitter  using  12  sub-carriers  and  modulated  by  the  sequence, 

S-26,26  =  713/6  {0,  0,  1+7,0,  0,  0,  -1-7,0,  0,  0,  1+7,0,  0,  0,  -1-7,0,  0,  0,  -1-7,0, 
0,  0,  1+7,0,  0,  0,  0,  0,  0,  0,  -1-7,0,  0,  0,  -1-7,0,  0,  0,  1+7,0,  0,  0,  1+7,0,  0,  0, 
1+ 7,  0,  0,  0,  1+ 7,  0,  0}.  The  reason  for  multiplying  the  coefficient  values  of  the 
sequence  by  ^13/6  is  to  nonnalize  the  power  of  the  resulting  OFDM  symbol  because 

only  12  of  52  sub-carriers  are  used.  The  resulting  symbol  is  mapped  to  the  appropriate 
input  of  the  64  point  IFFT  that  serves  as  the  kernel  for  the  transmitter  OFDM  processing. 
The  mapping  procedure  is  depicted  below.  [4] 


Tint  Domain  Outputs 


Figure  18  Coefficient  mapping  within  the  IDFT  -  from  Reference  [4] 


In  Figure  18,  the  values  for  coefficients  1...26  are  mapped  to  the  same  numbered 
inputs  while  the  coefficients  -26... -1  are  mapped  to  inputs  3 8... 63.  The  remaining 
inputs,  27... 37  and  0  (dc)  are  set  to  zero.  [4] 

The  in-phase  component  of  the  time  domain  output  of  the  IFFT  is  shown  in 
Figure  19. 
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Time  (in  seconds)  xIO  6 


Figure  19  Short  training  symbols  (time  domain,  in-phase) 


As  depicted  in  Figure  19  and  in  accordance  with  [4],  the  output  of  the  IFFT  has  been 
cyclically  extended  to  161  samples,  or  8  //sec,  and  the  following  windowing  function  has 
been  applied: 


PT{n)  =  < 


.5 

1 

.5 


n  =  0 

1  <  //  <  1 60  > 
n  =  160 


(4.1) 


In  Equation  (4.1),  //represents  the  sample  index  and  PT  is  the  windowing 

function  discussed  in  Equation  (2.6),  now  a  function  of  sample  index  vice  time.  From 
Figure  19,  the  short  training  symbols  are  repetitive  with  ten  repetitions,  each  16  samples 
in  duration.  Since  the  802.1  lg  receiver  knows  a  priori  the  sample  values  of  the  training 

symbols,  this  can  be  leveraged  to  help  achieve  synchronization  [13]. 
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3.  Long  Training  Symbols 

The  long  training  symbols  are  included  in  the  OFDM  packet  for  fine  frequency 
synchronization  and  are  concatenated  to  the  short  training  symbols  [4],  Note  that  the 
long  training  symbols  include  a  cyclic  prefix  designed  to  protect  the  long  training 
symbols  from  ISI  generated  by  the  short  training  symbols.  Figure  17  indicates  that  the 
duration  of  the  cyclic  prefix  is  1.6  //sec.  This  duration  compares  favorably  to  the  values 
of  delay  spread  in  Table  1.  The  two  long  training  symbols  are  identical  and  each  symbol 
is  3.2  //sec  long,  resulting  in  a  long  training  symbol  total  duration  (including  cyclic 
prefix)  of  8.0  //sec.  The  long  training  symbols  are  generated  by  the  802.1  lg  transmitter 
using  the  sequence  L_26  26  =  {1,  1,  -  1,  -  1,  1,  1,  -  1,  1,  -  1,  1,  1,  1,  1,  1,  1,  -  1,  -1,  1,  1, 


1,  1,  -  1,  1,  1,  1,  1,  0,  1,  -  1,  -  1,  1,  1,  -  1,  1,  -  1,  1,  -  1,  -  1,  -  1,  -  1,  -  1,  1,  1,  -  1, 


-  1,  1,  -  1,  1,  -  1,  1,  1,  1,  1}  [4],  After  mapping  to  the  appropriate  IFFT  input  depicted 
in  Figure  18,  the  time  domain  output  is  cyclically  extended  to  161  samples,  creating  the 
required  cyclic  prefix  and  ensuring  a  total  duration  of  8.0  //sec.  [4] 


The  time  domain  representation  of  the  long  training  symbols  appears  in  Figure  20. 
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Figure  20  Long  training  symbols  (time  domain,  in-phase) 


In  Figure  20,  the  time  axis  for  the  long  training  symbols  reflects  the  position  of  the  long 
training  symbols  within  the  OFDM  packet.  The  last  sample  is  at  16  //sec,  reflecting  the 
last  opportunity  to  use  preamble  data  for  synchronization.  The  802.1  lg  receiver  uses 
these  training  symbols  under  the  same  principle  as  the  short  training  symbols  -  a  priori 
knowledge  of  the  long  training  symbols  can  used  to  achieve  synchronization  [8], 

4.  Pilot  Tones 

Each  OFDM  data  symbol  contains  4  dedicated  sub-carriers  that  also  carry  a  priori 
sample  values.  These  pilot  tones  are  used  to  prevent  loss  of  phase  and  frequency 
synchronization.  The  pilots  are  inserted  into  sub-carriers  -21,  -7,  7,  and  21  prior  to  IFFT 
processing.  To  reserve  those  sub-carrier  positions,  the  sub-carriers  carrying  data,  the 

logical  sub-carriers  0. .  .47,  are  remapped  using  the  mapping  function 
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(4.2) 


i  -  26  0  <  /  <  4 

i-25  5  <  /  <1 7 
i-24  18  <  i  <  23 

i-23  24  <  /  <  29 
i-22  30  <  /  <  42 

i- 21  43  <  /  <  47. 

In  Equation  (4.2)  i  represents  the  sub-carrier  index.  Equation  (4.2)  shows  that 
positions  -21,  -7,  7,  and  21  are  skipped  to  provide  for  pilot  insertion.  The  pilot  a  priori 
values  are  {1,1, 1,-1}  or  {-1,-1, -1,1}  respectively.  [4] 

The  polarity  of  the  pilot  tones  is  dependent  upon  the  first  chip  of  an  internally 
generated  pseudorandom  scrambling  sequence  [4].  The  scrambling  is  performed  to 
avoid  exceeding  the  average  power  spectral  density  of  the  OFDM  spectrum  during 
transmission  [13]. 

5.  Summary 

Due  to  the  nature  of  the  packet-based  WLANs,  synchronization  must  occur 
rapidly  and  includes  the  detection  of  the  actual  packet  as  the  first  synchronization  task. 
Furthermore,  because  of  the  nature  of  OFDM  modulation,  specific  synchronization 
techniques  must  be  employed  that  differ  from  techniques  in  single  frequency 
communication.  The  governing  IEEE  standard  does  not  explicitly  mandate 
synchronization  algorithms  [8],  but  data-aided  algorithms  appear  to  be  best  suited  for 
packet-based  WLANs  and  generally  provide  good  perfonnance  [12].  Furthermore,  the 
IEEE  standard  does  mandate  various  training  symbols  and  pilot  tones  to  aid  data-aided 
synchronization.  Since  this  overhead  is  mandatory,  there  is  no  additional  cost  for  data- 
aided  synchronization  to  in  terms  of  throughput.  The  remainder  of  this  chapter  will  focus 
on  synchronization  algorithms  that  use  training  symbols  in  the  OFDM  packet’s  preamble 
and  embedded  pilot  tones.  Figure  21  depicts  the  time-frequency  structure  of  an  802.1  lg 
OFDM  packet,  where  all  training  symbols  and  pilots  are  shaded  in  red.  It  clearly 
illustrates  the  emphasis  placed  by  the  IEEE  standard  on  data-aided  synchronization. 


MM  =  \ 
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-8.125  MHz 


Figure  21  Time-frequency  structure  of  OFDM  802. llg  packet  -  from  Reference 

[20] 
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C.  IEEE  802.1 1G  PACKET  DETECTION 

1.  Introduction  and  Link  Budget  Analysis 

Packet  detection  means  finding  the  beginning  of  the  OFDM  packet.  Without 
packet  detection,  the  follow-on  synchronization  functions  cannot  be  performed.  [8] 

Furthermore,  as  discussed  previously,  the  packet-based  nature  of  WLANs 
requires  that  synchronization  be  performed  rapidly  once  the  packet  has  been  detected. 
The  minimum  average  received  signal  power  that  the  802. 1 1  g  OFDM  receiver  must 
detect  before  bit  error  rate  becomes  unacceptable  is  known  as  receiver  sensitivity.  The 
sensitivity  is  determined  by  the  IEEE  standard  and  is  depicted  in  Table  2  [4]. 


Table  2  Receiver  Performance  Requirements  -  from  Reference  [4] 


Data  rate 

(Mbps) 

Minimum  Sensitivity 

(dBm) 

6 

-82 

9 

-81 

12 

-79 

18 

-77 

24 

-74 

26 

-70 

48 

-66 

54 

-65 

These  mandated  receiver  sensitivities  will  provide  good  bit-error-rate 
performance  under  typical  indoor  noise  and  channel  conditions  with  the  range  increasing 
coupled  with  decreasing  data  rate  [34].  Table  2’s  receiver  perfonnance  requirement  is 
predicated  on  a  packet  error  rate  of  less  than  10  %  at  a  PDSU  length  of  1000  bytes  [5]. 
The  table  implies  that  higher  data  rates  are  associated  with  greater  signal  strength  at  the 
receiver.  Received  signal  strength  decreases  rapidly  with  increasing  range  and 
obstructions  and  is  known  as  path  loss  [35]. 
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A  common  expression  for  2.4  GHz  radio  frequency  path  loss  is 

/5Z(dB)  =  /,Z(l)  +  101og(Dx)  (4.3) 

where  PL  is  the  path  loss  in  dB,  D  is  distance  in  meters,  and  y  is  the  path  loss 
exponent.  The  value  of  y  depends  on  the  attenuation  characteristics  of  the  wireless 
environment  -  y  =  2  for  open  space,  y  =  3.3  for  open  office,  and/  =  4.5  for  home 
environments.  .PZ(l)  is  the  path  loss  at  1  meter,  which  is  equal  to  41  dB.  [35] 

The  power  received  for  a  typical  802. 1  lg  system  can  be  written  as 


p  _  P„GnGn 

rx  PL 


In  Equation  (4.4),  Prx  is  signal  power  out  of  the  receiver’s  antenna,  Pa  is  transmitted 
power,  Glx  is  the  transmitter’s  antenna  gain,  and  Grx  is  the  receiver’s  antenna  gain. 
Equation  (4.4)  may  be  re-written  as 

/>„(dB)  =  Pt,(dB)  +  G„(dBi)  +  Gt.(dBi)-PI(dB).  (4.5) 

Equation  (4.5)  shows  that  received  power  decreases  with  increasing  path  loss.  Equation 
(4.3)  showed  that  this  path  loss  was  dependent  upon  the  characteristics  of  the  wireless 
channel,  but  nonetheless,  exponential  with  increasing  distance.  The  detection  range  of 
the  802.1  lg  system  could  be  increased  by  improving  receiver  sensitivity.  For  indoor 
environments,  Equation  (4.4)  can  be  re-written  to  solve  for  distance,  D. 

D  =  4.5  P*G«Gr*  (4.6) 

V  WW) 


In  Equation  (4.6),  where D  is  in  meters,  if  Ptx,Gtx,  and  Gn  are  held  constant,  then  to 
double  reception  range  would  require  an  improvement  in  receiver  sensitivity,  or  decrease 
of  acceptable  Prx,  by  a  factor  of  approximately  22.6  (=245),  or  13  dB.  For  outdoor 

environments,  the  perfonnance  requirement  would  not  be  as  onerous,  due  to  the 
decreased  path  loss  exponent.  Figure  22  depicts  the  increased  reception  range  with  lower 
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minimum  receiver  sensitivity,  assuming  the  bit-error-rate  (BER)  can  be  held  to 
acceptable  standards.  The  plot  is  based  on  the  path  loss  model  from  Equation  (4.3)  using 
typical  commercial  sector  values  for  Ptx,  Gtx,  and  Grx  [36]. 


Figure  22  Power  received  as  a  function  of  distance 


While  doubling  range  performance  requires  a  13  dB  decrease  in  sensitivity,  the 
same  analysis  would  show  that  a  3  dB  decrease  in  receiver  sensitivity  yields,  in  an  indoor 
environment,  a  commercially  significant  increase  of  16  percent  in  reception  range.  The 
analysis  depicted  in  Figure  22  only  considers  increasing  receiver  sensitivity  to  achieve 
increases  in  range  performance.  Commercial  vendors  obviously  have  the  option  to 
combine  increased  transmitter  power  or  improved  antenna  gain  to  achieve  the  same 
effects.  This  thesis  will  only  focus  on  synchronization  techniques  to  improve  receiver 
sensitivity  and  thereby  increase  reception  range. 
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2.  Packet  Detection  Theory 

Packet  detection  is  based  on  two  complementary  hypotheses  -  either  a  packet  is 
present  or  it  is  not  [8].  In  order  to  build  up  to  overall  packet  detection,  we  begin  by 
considering  the  simpler  detection  problem  of  packet  detection  using  the  received  signal 
value  at  single  sample  instant,  tm .  This  is  can  be  mathematically  represented  as  [12] 

Ho  ■r(tm)  =  n(tm)  (4.7) 

Hi  :K0  =  s(0  +  n(0  (4-8) 

where  H0 is  hypothesis  zero  (i.e.  no  packet  is  present,  only  noise  present)  and//,  is 
hypothesis  one  (a  packet  is  present  with  noise),  r(tm)is  the  received  signal,  n(tm) is 
noise,  and  s(tm)is  the  incoming  packet.  The  receiver  should  choose  the  correct 
hypothesis  with  high  probability,  i.e.  decide  if  a  packet  is  present  (D,)or  if  a  packet  is 
not  present  (D0). 

The  goal  is  for  the  receiver  to  maximize  the  probability  that  a  packet  is  detected  in 
the  presence  of  noise,  while  minimizing  the  probability  that  a  receiver  mistakenly 
declares  a  packet  present  when  only  noise  is  present.  These  probabilities  are  depicted  in 
Figure  23.  The  horizontal  axis  shows  the  values  for  decision  statistic %.  Ordinarily^  is 
computed  by  the  receiver,  but  for  this  initial  step  in  the  development  of  a  packet  detector, 
X  is  merely  a  sample  of  the  received  signal.  In  other  words,  xm  =  r( tm ) ,  and  is 

dependent  on  whether  r  ( tm  )  contains  only  noise  (Z/0)  or  the  incoming  packet  and  noise 
(//,).  xm  is  compared  to  threshold,  yT,  to  allow  the  receiver  to  make  decision 
D0  or/),.  The  blue  curve  on  the  left  of  Figure  23  represents  the  probability  density 
function  (PDF)  for  xm  given  hypothesis  H0  and  is  called  the  likelihood  of  H0.  The  red 
curve  on  the  right  represents  the  PDF  for  xm  given  hypothesis  //,  and  is  called  the 
likelihood  of  //,.  [10] 
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The  curves  in  Figure  23  depict  classic  Gaussian  probability  density  functions. 
This  is  because  the  values  of  n(tm)  are  Gaussian  distributed.  Because  %m  is  simply  a 

linear  transformation  of  n(tm),  it  is  also  Gaussian  distributed.  The  probability  density 
function  can  be  mathematically  described  as 


P(Zm) 


where  p(%m)  's  the  notation  used  to  describe  the  probability  density  function  of 
Gaussian  random  variable  %m,  mx  is  the  mean  of  and  <J2X  is  the  variance  of  %m. 
Note  that  my  is  zero  under  H0  and  non-zero  under  Hv  [37] 


Figure  23  Probabilities  of  False  Alarm  and  Detection  in  an  AWGN  channel 
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To  maximize  the  probability  that  the  receiver  will  detect  a  packet  in  the  presence 
of  noise  is  to  maximize  Pd,  which  is  known  as  the  probability  of  detection  and  is  defined 
as 

Pd  =  PrjDj  \H^.  (4.9) 

To  minimize  the  probability  that  the  receiver  will  not  mistakenly  declare  a  packet  present 
when  only  noise  is  present  is  to  minimize  the  probability  of  false  alarm,  Pfa,  which  is 

defined  as 

pf.=HD<  N  <4-10> 

[37] 

As  Figure  23  shows,  both  probabilities  cannot  be  optimized  simultaneously.  Pd  is 
the  area  bounded  by  the  likelihood  of  Hl  on  the  right  and  yT  on  the  left.  Pfa  is  the  area 
bounded  by  likelihood  H0  on  the  right  and  yT  on  the  left.  Increasing  Pd  by  lowering  the 
detection  threshold,  yr,  results  in  an  increase  in  the  probability  of  a  false  alann. 
Conversely,  minimizing  Pfa  results  in  lowered  Pd  The  relationship  between  Pd  and  Pfa , 
and  therefore  D0  and  D,,  is  dependent  upon  and  yT.  The  relationship  can  be 
expressed  as 

A 

zm><yT  (4-ii) 

Do 

[10].  The  inequality  relationship  in  Equation  (4.11)  indicates  that  the  receiver  makes 
decision  Dx  if  %m  >  yT  or  D0  if  xm  <yT ■  To  achieve  maximum  performance  from  the 

receiver,  %m  must  be  computed  with  sufficient  accuracy.  Furthennore,  a  value  of  yT 
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must  be  chosen  to  minimize  error  and  optimize  detection.  While  this  thesis  will  explore 
various  synchronization  algorithms  that  are  used  to  compute  %m ,  we  will  first  evaluate  an 
optimal  value  for  yT . 


a.  Minimizing  Error  in  Packet  Detection 

One  popular  criterion  to  optimize  receiver  perfonnance  is  to  minimize  the 
probability  of  receiver  error.  To  minimize  the  probability  of  error,  the  inequality 

A 

Pr{X<yT\H0}  >  PrjZ/J 
Pr{X>yT\H,}  <  Pr{//0} 

A 

is  formed.  Equation  (4.12)  is  known  as  the  likelihood  ratio  test.  A  is  the  ratio  of  the 
conditional  probabilities  Pr[%m  <  yT  |  H(j  j  and  Pr  {  %m  >  yr  |  A }  for  random  variable  % 

at  sample  time  t  =  m  given  hypothesis  H0  or  //, ,  and  threshold  yT.  Pr  {//0  }  and 
Pr{//j}  represent  the  a  priori  probabilities  for  H0  and  //,  respectively.  These 
conditional  probabilities  represent  the  likelihood  of  H0  or  Hv  [10] 

If  the  a  priori  probabilities  of  //n  or  II ]  are  unknown,  often  the  a  priori 

probabilities  are  set  equal  to  each  other.  In  that  case,  Equation  (4.12)  may  be  re-written 
as 

A 

Pr{zm  <  yT  | A,}  ^  Pr{zm  >  yT  |A}  (4-13) 

A 

Receiver  decision-making  based  on  Equation  (4.13)  is  known  as  the  maximum  likelihood 
criterion.  [10] 

The  maximum  likelihood  criterion  can  be  adjusted  to  fit  other  hypothesis 
testing  algorithms,  such  as  Bayesian,  min-max,  and  Neyman-Pearson.  The  Bayesian 
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criterion  assigns  cost  functions  to  incorrect  decisions  by  the  receiver  with  known  a  priori 
probabilities  of  the  hypotheses.  The  min-max  criterion  is  similar  to  the  Bayesian 
criterion  but  is  based  on  unknown  a  priori  probabilities.  The  Neyman-Pearson  criterion  is 
based  on  maximizing  Pd  for  a  fixed  Pfa.  All  three  criterions  result  in  a  likelihood  ratio 

test,  similar  to  Equation  (4.11)  but  with  different  yT.  In  those  cases,  the  value  of  yT 
could  be  dependent  on  cost  functions  associated  with  incorrect  decisions  made  the  by  the 
receiver,  a  priori  probabilities  of  H0  or  Hv  or  on  a  fixed  value  for  Pfa.  [37] 

Because  of  the  desire  to  maximize  Pd  to  improve  performance,  the 
Neyman-Pearson  criterion  will  be  the  focus  of  the  remaining  part  of  this  section. 


b.  Neyman-Pearson  Criterion 

Because  the  Neyman-Pearson  criterion  assigns  a  fixed  probability  of  false 
alarm  and  then  attempts  to  maximize  probability  of  detection  using  a  likelihood  ratio  test, 
the  likelihood  ratio  now  becomes 


D, 


A  = 


Pr{x,n  >  yT\Hl) 

Pr{zm<yT  |#o} 


> 


< 

Dn 


(4.14) 


where  Ar  represents  the  required  likelihood  ratio  threshold  and  A  is  the  likelihood  ratio 
[12].  at  is  determined  by  false  alarm  or  detection  probabilities  and  is  not  equivalent  to 
yT ,  the  signal  threshold.  If  the  receiver  passes  only  RF  in  the  band  of  interest,  the  noise 
tenns  in  Equations  (4.7)  and  (4.8)  are  assumed  to  be  band  limited  AWGN  (BLAWGN), 
with  zero  mean  and  a  variance  a2,  then  the  conditional  probability  density  function  of 
random  variable  %m  given  hypothesis  Hx  can  be  described  as  [37] 

W.U-)  =  ^iU-)  =  T=Te  2CT"  '  (4-15) 

pnal 
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Equation  (4.15)  shows  that  %m  is  Gaussian  distributed  with  its  mean  value  equal  to  the 
signal  value  sm  at  sample  time  t  =  m  [12].  The  probability  density  function,  under 
hypothesis  H0  would  be  [37] 


Zm  <yr  \H0 


{xm)  =  P0(xm)  = 


-(xm  f 

2  al 


(4.16) 


In  Equation  (4.16),  %m  remains  Gaussian  distributed  at  sample  time  t  =  m,  with  the  same 
variance  but  zero  mean  in  accordance  with  H0.  The  likelihood  ratio  of  the  two  density 
functions  is 


Aizm) 


_  PiUJ  _ 

P0(X,n)  1 


W 


=  e 


2<j‘ 


CT.,  L  J 


(4.17) 


Taking  the  natural  log  of  both  sides  and  simplifying  yields 


In ( A ( ))  =  [sm  ( 2 xm  ~  K  )] • ■ 


(4.18) 


Substituting  Equation  (4.18)  into  the  left  hand  side  of  Equation  (4.14)  and  solving  for 
Xm  yields 


A 

<  S,n  2 

Dn 


Equation  (4.19)  shows  that  yT  does  not  equal  Zr  [12].  Setting 


yr  = 


\SmJ 


ln(4)  +  y 


(4.19) 


(4.20) 
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yields  the  hypotheses  depicted  in  Figure  23  and  mathematically  expressed  as  Equation 
(4. 1 1).  As  Figure  23  depicts,  the  probability  of  false  alarm  may  be  defined  as 

Pfa=£iPAzm)dZm  (4-21) 

and  the  probability  of  detection  may  be  defined  as 

Pd=[yPAxm)dxm  (4-22) 

[37].  Equations  (4.21)  and  (4.22)  show  that  a  desired  Pfa  or  Pd  can  be  achieved  by 
specifying  a  particular  value  of  yT.  This  is  illustrated  for  clarity  in  Figure  24,  where  a 

detection  threshold  is  depicted  for  a  hypothetical  decision  statistic.  The  decision  statistic 
is  shown  in  blue  before,  during  and  after  a  packet  is  present  at  the  receiver.  Figure  24 
suggests  that  a  lower  detection  threshold  will  result  in  a  greater  Pd  but  also  result  in 

greater  Pfa.  The  desired  Pd  will  be  driven  by  perfonnance  metrics  for  packet  detection 

for  a  given  signal  strength  in  the  presence  of  noise,  which  is  based  on  the  receiver’s 
sensitivity. 


58 


Figure  24  Signal  threshold  given  hypothesis  H0  andH{ 


Figure  24  also  suggests  that  separating  the  probability  density  functions 
for  hypotheses  H0  and  //,  would  increase  Pd  while  lowering  Pfa.  Changing  the 

probability  density  function  for  Hl  could  be  easily  done  simply  by  increasing  transmitter 

power.  The  IEEE  standard,  however,  limits  maximum  transmitted  power  in  accordance 
with  govermnent  regulations  [5].  Additionally,  commercial  vendors  limit  transmitted 
power  to  meet  power  management  schemes  [31],  and  must  limit  large  peak- average 
power  ratios  within  the  OFDM  transmitter  to  maintain  performance  and  acceptable  cost 
of  the  digital-to-analog  converter  (DAC)  [31]. 

Another  way  to  change  the  PDFs  for  H0  and  //,  would  be  for  the  receiver 
to  use  multiple  signal  samples  to  compute  the  decision  statistic.  This  would  have  the  net 
effect  of  shifting  the  PDF  for  lix  away  from  the  PDF  for  H0  while  the  detection 
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threshold  could  be  maintained  at  the  same  level.  The  variance  for  //,  and  H0  would  be 
increased  but  Pd  would  also  be  increased  without  increasing  Pfa.  [12] 

The  OFDM  802.1  lg  packet  facilitates  the  use  of  multiple  samples  for 
packet  detection.  The  “multiple  samples”  are  the  samples  of  the  repetitive  training 
symbols  [31]  that  are  mandated  at  the  start  of  the  packet  and  shall  be  the  focus  of  the  next 
section. 


c.  Using  Multiple  Samples  to  Facilitate  Packet  Detection 

The  previous  section  highlighted  the  fact  that  multiple  samples  of  the 
incoming  signal  could  be  used  to  increase  Pd  while  holding  Pfa  constant.  For  multiple 

samples,  the  noise  is  uncorrelated  because  it  is  assumed  white  [12]. 

With  multiples  samples,  the  two  hypotheses  from  Equations  (4.7)  and 

(4.8)  become 

H0-rm=n,n  (4-23) 

Hl'rm=Sm+nm  (4-24) 

for  all  sample  times  m  =  |l,2,...,M -1  ,M}.  Equations  (4.23)  and  (4.24)  show  that  rm  is 
the  sum  of  Gaussian  random  variable  nm  and  its  mean  value,  sm ,  when  a  packet  is 
present.  Hence,  rm  is  Gaussian  distributed  with  a  mean  of  zero  under  H0  and  a  mean  of 
sm  under  Hv  The  variance  of  rm  under  //,  is  given  by 


=E\jm -£(OT  =E[sm  +nm -smf  =E{nm)2  =  a2„, 


(4.25) 


where  cr  is  the  noise  variance  under  the  same  band  limited  conditions  discussed 
previously.  [37] 

The  joint  probability  density  functions  are 


A)(o=n 


-G,)2 

.  2  a, : 


(4.26) 
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where  X  is  the  new  decision  statistic  based  on  the  sum  of  individual  values  of  rmsm  taken 
from  sample  times  1,2,..., M  -  l,M.  [37] 

Note  that  Equation  (4.30)  simplifies  to  (4.11).  Because  the  probability 
functions  for  H0  and  //,  are  known,  yT  can  be  determined  by  specifying  the  desired 
probability  of  false  alarm  or  the  probability  of  detection  [12]. 

Due  to  multiple  sampling,  Equation  (4.21)  becomes 
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^=j>.(xyx 


(4.31) 


and  Equation  (4.22)  becomes 


^=(>,(xyx. 


(4.32) 


Considering  X  as  a  Gaussian  random  variable  because  it  is  the  sum  of  Gaussian  random 
variables  [12],  Equation(4.32)  may  be  re-written  as 


p,-r 

J  y 


-(X-mx  )2 

dX, 


(4.33) 


where  mx  is  the  mean  and  crx  is  the  variance  of  X. 


The  area  to  the  right  of  any  threshold  for  a  zero  mean,  normalized, 


Gaussian  random  variable  is  defined  as 


\  oo 

Q^yT^  =  ~rxz\v-e  2 du 


(4.34) 


which  is  sometimes  referred  to  as  the  complementary  cumulative  distribution  function  for 
a  zero  mean,  unit  variance  Gaussian  random  variable  [34].  Substituting  Equation  (4.33) 

into  Equation  (4.34)  by  letting  u  =  ^  m^)/  yields  the  following  relationship: 


Pd  =  Q\ 


(4.35) 


Solving  for  the  threshold,  Equation  (4.35)  becomes 

Xr  =  mx  +  axQ  {Pd')' 


(4.36) 


In  Equation  (4.36),  mx  is  defined  as 


mx  = 


I  M  I  I  M  I  M  I  M 

=  £{X)  =  E  £  j.,  =  E  £  Vm  =£  £(i.  +a.)i.  =  2{i+s.,EM]  =  Afe2(4.37) 
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where  the  last  step  holds  if  sm  =  s  is  a  constant  and  M  is  the  number  of  samples  used  in 
the  detection  algorithm.  The  variance  of  X  is 


<=E- 


=  E< 


M 

Tj(S  +  n>n)S~nh 


“1 2 


m= 1 
f  M 


X-s'2 +'s'Zn»  -Ms'2 


=  E< 


V  m—\  m= 1 


(  M  A 

V  W=1  J 


Equation  (4.38)  can  be  further  simplified  to 


2  2  77 

<JX=  s  t< 


f  M  M  V 

2X 

Vm=l  J\l= 1  / 


m=l  /=1 
M  M  _ 

W=1  /=1 

=  s2[m(<t2) 


=  52Mct; 


where  /n  is  the  Kronecker  delta  function  and  is  defined  as 


1  if  /  =  m 
0  if  /  ^  in 


(4.38) 


(4.39) 


(4.40) 


[12] 

Equation  (4.37)  shows  that  the  mean  value  of  X  under  Hx  can  be  increased 
by  using  multiple  samples.  Equation  (4.39)  shows  that  this  does  increase  the  variance  of 
the  decision  statistic  under  Hv  Using  Equations  (4.36),  (4.37),  and  (4.39),  the 

threshold, yT,  becomes 
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yT  =  Ms2  +  Vm so„Q 1  (Pd  ) 


(4.41) 


[12]. 


For  hypothesis  H0,  the  mean  of  X  is 

mx=E\Yjnm  1  =  0. 

t  m=l  J 


(4.42) 


The  variance  of  X  under  //0is  the  same  as  Equation  (4.39).  Substituting  Equation  (4.31) 
into  Equation  (4.34)  yields 


Pfa=Q 


\axJ 


(4.43) 


which,  after  making  substitutions  from  Equations  (4.39)  and  (4.41)  simplifies  to 

f  _  A 


Pfa=Q 


4m—+q-'(pj) 

V  an  J 


(4.44) 


[12] 

Equation  (4.35)  shows  that  a  desired  Pd  can  be  obtained  by  adjusting  the 
value  of  yT.  Using  multiple  samples  to  accomplish  this,  allows  us  to  lower  the  threshold 
to  obtain  a  desired  Pd  while  minimizing  Pfa.  There  are  packet  detection  schemes  that 
incorporate  this  algorithm.  These  schemes  shall  be  covered  in  detail  in  the  next  sections. 


3.  Delay  and  Correlate  Packet  Detection 

The  Delay  and  Correlate  packet  detector  is  one  example  of  a  robust  packet 
detection  scheme  that  employs  an  algorithm  similar  to  one  discussed  previously  in  this 
chapter.  The  Delay  and  Correlate  detector  is  based  on  a  double  sliding  window  algorithm 
[8]. 

The  double  sliding  window  algorithm  calculates  received  signal  energy  for  two 
consecutive  sliding  windows  of  samples.  A  comparison  of  the  ratio  of  received  energy  in 
both  windows  is  used  to  compute  the  decision  statistic.  The  Delay  and  Correlate 
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detector’s  windows  are  samples  that  comprise  the  signal  and  a  delayed  version  of  the 
signal.  This  detection  scheme  will  exploit  the  repetitive  preamble  data  (short  training 
symbols)  at  the  start  of  the  OFDM  packet.  The  signal  flow  and  structure  of  the  detector 
is  illustrated  below.  [8] 
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Figure  25  Signal  flow  structure  of  the  Delay  and  Correlate  algorithm  -  from 

Reference  [8] 


Figure  25  shows  that  the  decision  statistic  is  equal  to 


X  = 


(4.45) 


where  M  is  the  number  of  samples  used  to  compute  the  decision  statistic.  The  delay, 
z~D  in  Figure  25,  is  used  to  create  the  delayed  version  of  the  received  signal.  To  take 
advantage  of  the  periodicity  of  the  short  training  samples,  the  value  of  the  delay  is  set  to 
16  samples.  The  C  and  P  blocks  are  the  two  sliding  algorithms  used  to  compare 
whether  a  packet  is  present.  The  C  algorithm,  as  shown  in  Figure  25,  represents  the 
cross  correlation  between  the  signal  and  the  delayed  version  of  the  signal  (hence  the 
name  Delay  and  Correlate  detector).  The  P  algorithm  represents  received  signal  energy 
during  the  cross  correlation  window  and  is  used  to  normalize  the  cross  correlation  value, 
so  it  is  not  dependent  upon  the  absolute  received  power  level.  [8] 
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Figure  26  shows  the  response  of  the  decision  statistic  when  the  802.1  lg  short 
training  samples  are  detected  at  SNR  =  10  dB  and  M  =  8.  Included  in  Figure  26  are  the 
results  of  the  cross  correlation  window  and  received  signal  energy  window  for  clarity. 
The  short  training  samples  are  embedded  in  AWGN  and  are  fully  correlated  at  sample 
1000,  with  the  start  of  the  preamble  at  sample  839. 


Moving  sum  Rx  power 


Figure  26  Delay  and  Correlate  algorithm  using  802.1  lg  preamble  and  showing 

pm->  X 


The  low  value  of  X  prior  to  sample  1000  is  explained  by  the  fact  that  the  cross 
correlation  window  consists  only  of  uncorrelated  noise,  whose  mean  value  is  zero.  At  the 
start  of  the  packet,  the  identical  short  training  samples  in  the  cross  correlation  window 
cause  a  sudden  increase  in  the  value  of  the  decision  statistic.  This  will  give  a  good 
estimation  of  the  start  of  the  packet.  [8], 
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Inspection  of  Figure  26  shows  that  the  response  of  Xpennits  the  creation  of  a 
threshold  to  declare  whether  a  packet  is  present  or  not.  As  previously  discussed,  the 
threshold  can  be  determined  from  a  desired  Pd  and  acceptable  Pfa.  In  this  case,  Equation 

(4.30)  can  be  re-written  as 


X  = 


ivi  —i 

I' 

k=0 
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M- 1 

II 

k= 0 


m+k-D 


|2  < 

Dn 


yT- 


(4.46) 


In  Equation  (4.46),  because  the  variance  of  the  numerator  and  denominator  on  the 
left  hand  side  of  the  inequality  are  much  smaller  than  their  means,  the  decision  statistic 
can  be  approximated  as  a  Gaussian  random  variable  [38].  In  that  case,  the  probability 
distribution  of  X  can  be  approximated  as 


(X-Mxf 

/W  =  TTe  24 

yj  2  7T(JX 


(4.47) 


where  mx  is  the  mean  value  and  <jx  is  the  variance  of  the  Gaussian  approximation  of  X  . 
According  to  [38], 


mx 


cr 


(4.48) 


and 


2<Jt  [( 1  +  mx )  o)crn  +  ( 1  +  2mx )  cr4n  _ 

M(°s+°n)2 


(4.49) 


In  Equations  (4.48)  and  (4.49),  a]  is  the  variance  of  short  training  samples  and  cr[  is  the 


variance  of  the  AWGN  in  the  channel.  Recognizing  that 

SNR  =  10*  log 


\anJ 


(4.50) 


Equation  (4.48)  may  be  re-written  as 
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mx 


jq5SVR/10 
^qSTW/IO  _j_  j 


(4.51) 


where  SNR  is  in  units  of  dB.  Equation  (4.49)  may  be  re-written  as 


M 


-nu 


( 1  +  mx )  1 0 SNR !  1 0  +  1  +  2  mx 

(10s™  + 1)2 


(4.52) 


Equation  (4.51)  shows  that  the  mean  can  be  made  to  be  only  dependent  upon 
SNR.  Equation  (4.52)  shows  that  the  variance  of  the  decision  statistic  is  inversely 
proportional  to  the  number  of  samples  used  to  compute  the  decision  statistic.  As  more 
samples  are  used,  this  decreases  cr^.  From  Equation  (4.43),  this  shows  that  Pfa  is 

reduced.  Using  Equations  (4.51),  (4.52),  (4.36),  and  (4.43),  a  designer  need  only  specify 
a  desired  Pd  for  the  Delay  and  Correlate  detector  and  for  a  given  SNR,  yT  and  Pfa  may 

be  computed. 

Figure  27  illustrates  this  below.  The  plots  in  Figure  27  depict  yT  and  Pfa  when 
Pd  is  fixed  and  SNR  is  varied.  The  top  plot  shows  yT  and  Pfa  when  8  samples  ( M  =  8  ) 
are  used  to  compute  The  bottom  plot  depicts  the  same  when  M  =  16.  Figure  27 
clearly  shows  that  using  an  increasing  number  of  samples  to  compute  X  raises  yT  and 
lowers  Pfa.  Close  inspection  of  Figure  27  suggests  that  this  appears  most  evident  at 

single  digit  SNR  when  the  performance  of  the  802.1  lg  system  is  expected  to  be  at  its 
worst.  This  means  that  the  Delay  and  Correlate  detector  is  a  promising  synchronization 
technique  for  packet  detection.  The  repetitive  short  training  symbols,  consisting  of  161 
samples,  easily  lend  themselves  to  this  type  of  algorithm  and  can  be  leveraged  to  increase 
range  and  performance. 
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SNR  (dB) 


Figure  27 


Detection  Threshold  and  Pfa  for  Delay  and  Correlate  Detector  when 


=  .9 


4.  Matched  Filter  Detector 

Another  type  of  detector  that  can  be  used  to  exploit  the  short  training  symbols  is 
the  Matched  Filter  detector.  The  Matched  Filter  detector  algorithm  is  based  on  using  the 
short  training  symbols  sample  values  to  detect  the  incoming  signal.  The  signal  structure 
and  flow  is  depicted  below. 
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Figure  28  Matched  Filter  Detector  -  from  Reference  [13] 


In  Figure  28,  the  receiver  uses  delay  blocks  to  multiply  stored  coefficients  of  short 
training  sample  values  (c0....cM_j)  with  received  values  rw  The  short  training 

samples,  as  specified  by  the  802. 1  lg  standard,  are  periodic  every  sixteen  samples  [4],  so 
in  Figure  28,  M  -  16.  These  multiplied  values  are  then  summed  to  compute  the  decision 
statistic.  There  are  10  short  training  symbols,  so  in  theory  more  samples  could  be  used, 
thereby  improving  perfonnance  but  at  the  cost  of  increased  receiver  complexity.  This 
tradeoff  will  be  addressed  when  synchronization  perfonnance  is  analyzed  in  the  Chapter 
0.  Figure  29  depicts  the  absolute  values  of  the  decision  statistic  when  SNR  =  10  dB  (as 
was  done  previously  with  the  Delay  and  Correlate  detector). 


70 


0.2 


Received  samples  (signal  +  noise) 


T - 1 - 1 - 1 - r 


-0.2 - 1 - I - i - £ 


500 

1000  1500  2000 

2500 

3000 

Sample  count  - 

- ► 

Matched  Filter  Output 

Figure  29  Decision  statistic  output  for  Matched  Filter  Detector 


Figure  29  shows  that  there  is  a  large  increase  in  the  value,  or  “spike”  of  the 
decision  statistic  when  the  portion  of  the  received  signal  containing  a  short  training 
symbol  is  correlated  with  the  Matched  Filter  coefficients. 

The  Matched  Filter  coefficients  are  complex  conjugates  of  the  known  training 
symbol  values.  As  can  be  seen  in  Figure  29,  the  10  short  training  symbols  are  easily 
detected  at  the  tested  SNR.  The  low  values  of  the  decision  statistic  between  the  “spikes”, 
or  sidelobes,  represent  partial  or  no  correlation  between  the  Matched  Filter  coefficients 
and  the  signal.  [13] 

From  Figure  29  it  is  apparent  that  the  Matched  Filter  detector  presents  a  viable 
technique  to  detect  the  802.1  lg  packet.  The  first  spike  will  appear  when  the  first  sixteen 
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short  training  samples  match  the  filter  coefficients,  indicating  the  start  of  the  packet. 
Subsequent  spikes  could  possibly  be  used  to  improve  timing  synchronization. 

Using  the  same  method  of  analysis  from  the  previous  section,  the  mean  value  of 
decision  statistic,  X,  becomes 

{m  'l  I'm  t  I'm  t 

=  \  =  E\Tj{Sm+nm)Sm  \-  (4-53) 

m= 1  J  L  nt=l  J  L  m=\  ) 


Because  the  noise  in  the  channel  is  independent  of  the  values  of  short  training  samples, 
Equation  (4.53)  can  be  simplified  to 


=  E  \  Tj(Sm  +  nm  K  \  =  (j„)  +  E(sm)E  (nm  )}  =  E  \  J 


k  m- 1 


m=\ 


(4.54) 


^  m- 1 


Equation  (4.54)  is  similar  to  Equation  (4.37)  except  for  the  fact  that  the  values  of  the 
short  training  samples  are  not  constant  during  the  window  of  summation.  The  values  of 
the  short  training  samples  are  known,  however,  and  the  expectation  of  their  summation  is 
simply  the  summation  of  the  sample  values.  Hence,  Equation  (4.54)  may  be  re-written  as 

M 

™x=2X-  (4-55) 

m= 1 


Equation  (4.55)  is  confirmed  by  [12].  The  variance  of  the  decision  statistic  can  be 
expressed  as 
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(4.56) 


In  Equation  (4.56),  8mp  is  the  Kronecker  delta  function  defined  in  Equation 

(4.40).  The  analysis  in  Equation  (4.56)  is  similar  to  Equations  (4.38)  and  (4.39),  but  once 
again,  sm  is  not  a  constant  value.  Now  that  the  mean  value  and  variance  of  the  decision 
statistic  for  the  Matched  Filter  detector  have  been  established,  Equations  (4.36)  and 
(4.43)  can  be  used  to  compute  detection  threshold  and  Pfa  for  a  given  Pd. 

In  the  case  of  the  Matched  Filter  detector,  the  coefficients  of  the  Matched  Filter 
are  the  complex  conjugates  of  the  known  periodic  short  training  samples.  Hence,  an 
actual  value  can  be  computed  for  Equation  (4.55).  In  this  case,  mx=.2031. 

Furthermore,  Equation  (4.56)  can  be  simplified  to  crx  =  .203  lcr .  Figure  30  shows  how 
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detection  threshold  and  probability  of  false  alarm  varies  with  SNR  when  Pd  =  .9.  If  the 
mean  of  the  signal  strength  is  expressed  in  volts,  then  the  computed  threshold  can  be 
expressed  in  volts  as  shown  below. 


Figure  30  Matched  Filter  yT  and  Pfa  for  varying  SNR 


Figure  30  suggests  that  at  low  SNR,  the  Matched  Filter  detector  has  a  relatively 
high  false  alarm  rate  and  a  low  threshold  of  detection  is  required  to  ensure  desired  Pd  is 

achieved.  As  Figure  30  shows,  with  increasing  SNR,  Pfa  decreases,  as  expected,  while 
the  detection  threshold  can  be  raised. 

5.  Summary  and  Comparison  of  Packet  Detection  Techniques 

Link  budget  analysis  suggested  that  improving  packet  detection  would  improve 
range  and  performance  of  the  802. 1 1  g  system.  This  led  to  a  review  of  basic  detection 
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theory,  the  role  that  packet  detection  plays  in  synchronization,  and  a  review  of  various 
packet  detection  algorithms.  The  Neyman-Pearson  algorithm  seems  to  be  promising 
because  it  focuses  on  maximizing  Pd  for  a  particular  acceptable  Pfa.  Costs  associated 

with  incorrect  decision  making  need  not  be  computed.  Finally  the  use  of  multiple 
samples  to  compute  the  decision  statistic  was  examined.  High  Pd  could  be  specified 
without  increasing  Pfa  by  using  multiple  samples.  The  Delay  and  Correlate  detector,  as 

described  in  [8],  was  examined  in  detail.  It  employs  multiple  samples  to  optimize 
performance.  Subsequent  analysis  depicted  in  Figure  27  showed  that  the  Delay  and 
Correlate  detector  could  meet  a  desired  Pd  at  low  SNR.  While  threshold  values  remained 

low,  Pfa  was  held  to  an  acceptable  value  by  using  multiple  samples.  The  Matched  Filter 

detector,  as  described  by  [13],  also  uses  multiple  samples  as  it  seeks  to  leverage  the  16 
sample  periodicity  of  the  short  training  samples.  Unfortunately,  inspection  of  Figure  30 
revealed  much  higher  values  of  Pfa  at  low  SNR  compared  to  the  Delay  and  Correlate 

detector  values  shown  in  Figure  27.  This  suggests  that  the  Matched  Filter  detector’s 
performance  lags  that  of  the  Delay  and  Correlate  detector.  This  shall  be  confirmed  in 
Chapter  V,  Synchronization  Perfonnance.  In  Chapter  V,  multiple  simulations  for  both 
detectors  shall  be  run  at  low  SNR  and  detection  and  false  alarm  shall  be  verified. 
Additionally,  the  Delay  and  Correlate  detector  and  the  Matched  Filter  detector  will  be 
evaluated  using  more  than  one  16  sample  short  training  symbol. 

D.  FREQUENCY  SYNCHRONIZATION 
1.  Introduction 

As  shown  in  Chapter  III,  Equation  (3.1),  OFDM  receivers  suffer  from  degradation 
in  SNR  due  to  lack  of  frequency  synchronization.  This  lack  of  frequency 
synchronization,  or  frequency  offset,  is  a  critical  factor  in  OFDM  system  design.  Failure 
to  synchronize  in  frequency  results  in  ICI  and  violates  the  orthogonality  of  the  sub¬ 
carriers.  [37] 
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a. 


Single  and  Multi-Carrier  Frequency  Synchronization 
Differences 


At  the  heart  of  nearly  all  synchronization  circuits  is  some  version  of  a 
phase-locked  loop  or  its  functional  equivalent  [10].  Many  single  carrier  systems  use  a 
narrow  band  phase-locked  loop  to  acquire  and  track  the  carrier,  thereby  maintaining 
frequency  synchronization  [12].  The  phase-locked  loop  is  used  to  generate  an  error 
signal  proportional  to  the  difference  between  the  received  frequency  and  a  reference 
frequency.  This  error  signal  is  then  used  to  drive  a  voltage-controlled  oscillator  to  match 
the  received  frequency  and  drive  the  error  signal  to  zero.  Unfortunately,  in  OFDM 
systems,  the  multiple  sub-carriers  present  in  the  received  signal  make  it  extremely 
difficult  to  create  a  single  error  signal  that  is  proportional  to  the  difference  between  the 
received  signal  and  a  reference  signal  [12].  Consequently,  various  frequency 
synchronization  techniques  have  been  proposed. 

b.  General  Frequency  Synchronization  Technique  in  OFDM 

To  achieve  frequency  synchronization  with  minimal  computational  effort 
by  the  OFDM  receiver,  the  synchronization  process  is  divided  into  two  phases  - 
acquisition  and  tracking.  This  also  serves  to  minimize  the  amount  of  redundant 
information  that  must  be  added  to  the  transmitted  packet.  In  the  acquisition  (or  coarse) 
phase,  an  initial  estimate  of  the  frequency  offset  must  be  made.  Consequently,  the 
acquisition  phase  is  characterized  by  more  complex  and  robust  synchronization 
algorithms  which  focus  on  those  portions  of  the  packet  that  contain  data  that  can  be 
leveraged  for  synchronization.  Once  acquisition  has  been  achieved,  tracking  (or  fine 
synchronization)  can  be  maintained  with  a  less  complex  algorithm  that  only  needs  to 
correct  small  frequency  deviations.  [39] 

Furthermore,  frequency  synchronization  techniques  can  be  characterized 
as  data-aided  (relying  on  special  training  symbols  embedded  within  the  packet)  or 
dependent  on  redundancy  (e.g.  cyclic  prefix  based)  within  the  packet  [33].  Data-aided 
techniques  will  achieve  synchronization  using  the  802.1  lg  mandated  training  symbols,  as 
discussed  in  the  previous  chapter.  Redundancy  techniques  exploit  the  cyclic  prefix  that  is 
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part  of  every  802. 1  lg  OFDM  symbol.  The  802. 1  lg  standard  does  not  specify  a  particular 
synchronization  technique;  however  latter  short  training  symbols  are  intended  by  the 
802.1  lg  standards  group  for  coarse  frequency  synchronization  while  the  two  long 
training  symbols  are  intended  for  fine  frequency  synchronization  [40].  As  discussed  in 
earlier  in  this  chapter,  data-aided  synchronization  offers  a  strong  candidate  for  OFDM 
WLAN  applications  because  of  the  repetitive  nature  of  the  training  symbols  within  the 
preamble.  The  cost  in  data  efficiency  must  be  borne  anyway  to  remain  in  compliance 
with  the  IEEE  standard.  Cyclic-prefix  based  synchronization  offers  the  advantage  of  not 
having  to  wait  for  the  arrival  of  the  OFDM  preamble,  but  the  802. 1  lg  short  training 
symbols  lack  a  cyclic  prefix.  .  Synchronization  schemes  exploiting  the  short  training 
symbols  are  categorized  as  data-aided  and  will  be  the  focus  of  this  thesis. 

c.  Summary  of  Various  Synchronization  Algorithms 

Per  the  IEEE  802.1  lg  standard,  any  coarse  algorithm  must  be  sufficiently 
robust  to  correct  for  a  maximum  frequency  offset  of  120  kHz,  as  shown  in  Equation  (3.2). 
Ideally,  the  frequency  offset  should  be  computed  prior  to  FFT  processing  by  the  receiver 
[8].  Several  frequency  synchronization  algorithms  have  been  suggested  for  OFDM; 
Daffara  and  Chouly  [41]  and  Moose  [16]  have  suggested  a  maximum  likelihood  (ML) 
estimator  to  compute  frequency  offset.  Li,  Liu,  and  Giannakis  [40]  have  recommended  a 
non-linear  least  squares  approach.  Schmidl  and  Cox  [38]  and  Van  de  Beek,  Sandell,  and 
Borjesson  [42]  have  suggested  a  CP  synchronization  scheme  that  exploits  the  redundancy 
between  samples  in  the  data  and  CP  portions  of  the  OFDM  frame.  The  first  two 
algorithms  shall  be  addressed  in  this  chapter  as  the  focus  remains  on  data-aided 
synchronization  schemes. 

2.  Maximum  Likelihood  Synchronization 

The  maximum  likelihood  estimator  is  one  of  several  promising  techniques  that 
have  been  presented  by  a  variety  of  researchers,  in  slightly  different  forms  [8].  The 
underlying  premise  is  to  exploit  the  repetitive  training  symbols  (data-aided)  or  repeated 
data  symbols  (redundancy  -  CP  based).  Any  phase  difference  between  the  sub-carriers  of 
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repeated  OFDM  symbols  can  be  used  to  estimate  the  frequency  difference.  This  is  a 
reasonable  assumption  because  from  Proakis  [11],  the  received  signal  may  be  expressed 
as 

r(t)  =  s(t  -  x)  +  n(t)  (4.57) 

where  r  is  the  propagation  delay,  n(t)  is  the  AWGN  of  the  channel  and 

s(t)  =  Re[s,(t)ej2*f“‘].  (4.58) 

In  Equation  (4.58),  s(t )  is  the  transmitted  packet  and  s,(t)  is  the  equivalent  low-pass 
signal  of  s(t),  and  ftx  is  the  carrier  frequency  of  the  transmitted  signal.  Consequently, 
Equation  (4.57)  may  be  re-written  as 

r(t )  =  Re  j[s,  (t  -  x)  +  n,  | 

=  Re  |  [(.s',  (t  -  x)  +  n,(t))e~j2^T\ej1^  (4.59) 

=  Re  [  [s,  (t  -  x  )e~7(i  +  n,  (t)e [j  ei2rcfa'  j 

where  (f>  is  the  carrier  phase  error  and  is  equal  to  2^/’tr  [11].  If  there  is  a  difference 
between  fa  and  frx,  the  receiver’s  center  frequency,  an  additional  error  will  result  as  the 

received  signal  is  down  converted  to  baseband  at  the  receiver.  This  is  caused  by  the 
oscillator  mismatch,  between  transmitter  and  receiver  [8].  The  preceding  statement 
ignores  any  type  of  Doppler  effect  between  receiver  and  transmitter.  As  shown  in 
Chapter  0,  these  effects  were  considered  negligible  at  typical  walking  speeds  and  are  not 
considered  for  the  remainder  of  this  thesis. 

Disregarding  noise  (letting  n,  (/  )  =  0)  and  letting  x  =  0  at  sample  time  t  =  mts  for 

simplicity,  once  the  signal  in  Equation  (4.59)  has  been  received  and  down  converted  to 
baseband,  it  could  be  multiplied  by  the  complex  conjugate  of  a  repeated  symbol.  If  there 
were  no  oscillator  mismatch  between  repeated  symbols,  i.e.  ftx  -  frx  =  0  =  A/,  then  the 

result  of  the  correlation  would  be  strictly  real.  If  there  were  an  oscillator  mismatch,  i.e. 
A/  ^  0,  then  any  remaining  imaginary  portion  is  dependent  upon  the  phase  error,  which 
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would  be  proportional  to  the  oscillator  mismatch.  This  assumes  that  any  fading  is  flat 
fading.  Consider  the  sampled  complex  baseband  signal  with  an  oscillator  mismatch 
between  transmitter  and  receiver  (disregarding  noise  again).  In  this  case,  the  sampled 
signal,  at  t  =  mts,  would  be 


r  —  s  gjZxfixint, e~j2)[f„mts 


r  =  se 


J2x(Af)mts 


(4.60) 


where  A f  is  the  frequency  offset  between  transmitter  and  receiver  caused  by  oscillator 
mismatch  and  ts  is  the  sample  duration  [8].  Doppler  effects  are  ignored. 


Using  intennediate  variable  z,  an  estimation  of  the  frequency  offset  can  be 
formed  as 
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(4.61) 


Equation  (4.61)  shows  that  intennediate  variable  z  consists  of  a  real  and  imaginary 
component.  If  there  were  no  frequency  offset  (A/  =  0)  z  would  be  purely  real.  Instead, 
the  imaginary  portion  of  z  is  dependent  upon  the  frequency  offset.  If  the  sampling  rate  is 
high  (ts  «:  l),  then  the  imaginary  portion  of  z  can  be  simplified  because  if  2tzA fDts  «:  1, 

then  sin ( 2n/\fDts )  ~  2n!\fDts  and  cos (2 7rAfDts)~  1.  Consequently,  we  can  estimate  the 
frequency  offset  from  the  imaginary  portion  of  z  as 
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(4.62) 


A/  = 


2  nDt, 


-Zz, 


where  Z  is  defined  as  the  angle  argument  of  complex  number  z.  [8] 

This  estimator  is  referred  to  as  the  maximum  likelihood  estimator  because  the 
received  values  can  be  modeled  using  the  conditional  probability  distribution  of  the 
received  values  given  A/,  or  PjrjA/}.  In  this  case,  A/ is  considered  unknown,  but 

deterministic.  A f  is  the  value  of  A f  that  maximizes  the  P\rn  | A/'} .  [11] 


a.  Limits  on  the  Operating  Range  of  the  ML  Estimator 

ML  estimators,  such  as  the  one  described  above,  can  only  correct  up  to  a 
certain  maximum  frequency  offset.  This  is  because 

Zz  =  27rDtsAf  (4.63) 


which  is  only  unambiguously  defined  in  the  range  [-71, 7r] .  If  the  frequency  offset 
between  repeated  samples  exceeds 

I  A/ 
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n 


9  rr  79/ 


9  Fit 


(4.64) 


then  the  frequency  offset  estimate  will  be  inaccurate.  Typically,  the  maximum  allowable 
frequency  offset  is  normalized  to  sub-carrier  spacing.  If  the  delay,  D,  is  equal  to  symbol 
length,  which  for  immediately  repeated  samples  means  D  =  1,  Equation  (4.64)  may  be  re¬ 
written  as 


A/ 


<-/. 

2Js 


(4.65) 


Hence,  this  type  of  estimator  is  usually  limited  to  a  maximum  frequency  offset  that  is  one 
half  of  the  sample  rate  if  the  frequency  correction  is  performed  prior  to  the  DFT 
demodulation  by  the  receiver.  [12] 

From  Equation  (4.64),  when  D  becomes  larger,  the  range  of  the  estimator 
is  reduced.  For  data-aided  synchronization,  if  the  short  training  samples  are  used  for  the 

80 


ML  estimator  ( D  =  16,  ts  =50  //sec)  [8],  then  Equation  (4.64)  shows  that  the  maximum 

frequency  offset  becomes  625  kHz.  Redundancy  synchronization  algorithms  that  use  the 
cyclic  prefix  have  a  delay  that  is  longer  than  the  symbol  period  and  consequently  have 
less  range  than  the  short  training  samples.  This  can  be  alleviated  by  increasing  the 
sampling  rate.  The  long  training  samples,  as  discussed  previously,  are  periodic  every  64 
samples  (D  =  64),  giving  the  estimator  a  range  of  up  to  156.25  kHz.  From  the  numbers 

discussed  above,  it  is  apparent  that  either  data-aided  synchronization  technique  meets  the 
802.1  lg  standard.  The  performance  such  a  time-domain  ML  estimator  will  be  discussed 
in  the  next  section. 

b.  Implementation  of  the  ML  Estimator 

As  discussed  on  page  75,  frequency  synchronization  is  usually 
accomplished  by  the  OFDM  receiver  in  two  phases  -  acquisition  and  tracking.  The  short 
training  symbols,  with  their  16  sample  periodicity  at  the  start  of  the  frame,  can  be 
exploited  to  accomplish  ML  frequency  synchronization  [8].  Because  the  short  training 
symbols  are  at  the  beginning  of  the  frame,  they  lend  themselves  well  to  the  acquisition 
task.  Figure  3 1  depicts  the  output  of  MATLAB  code,  in  which  the  ML  algorithm  has 
been  used  to  generate  a  frequency  offset  correction  that  is  applied  to  subsequent 
incoming  short  training  samples.  The  top  graph  in  Figure  31  shows  the  received  short 
training  symbols  (161  samples)  after  passing  through  an  AWGN  channel  (SNR  is  10  dB). 
The  middle  graph  shows  those  same  received  samples  when  the  OFDM  receiver’s 
oscillator  is  experiencing  a  frequency  offset  of  40  percent  of  sub-carrier  spacing.  Recall 
from  Equation  (3.2)  that  the  worst  case  tolerable  frequency  offset  per  the  802. 1 1  g 
standard  is  120  kHz,  which  represents  only  38.4  percent  relative  of  the  sub-carrier 
spacing.  So  the  middle  graph  represents  a  worst  case  scenario  for  frequency  offset, 
rounded  up  for  simplicity.  To  facilitate  comparison  between  the  ML  algorithm  and  other 
frequency  synchronization  techniques,  a  standardized  frequency  offset  of  -40%  of  sub¬ 
carrier  spacing,  or  125  KHz  is  employed.  Such  a  frequency  offset  will  cause  a 
constellation  rotation  of  2 n  during  the  160  short  training  symbol  preamble,  resulting  in 
a  severe  degradation  in  performance.  The  bottom  graph  of  Figure  31  shows  the  short 
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training  samples  with  the  ML  frequency  offset  correction  applied.  To  show  the  effects  of 
the  ML  algorithm  on  short  training  symbol  reception  for  clarity,  the  ML  algorithm  uses 
an  “on-the  fly”  correction  that  applies  the  computed  frequency  offset  correction  to 
subsequent  incoming  short  training  symbols.  This  serves  to  significantly  mitigate  the 
effects  of  the  original  frequency  offset  despite  the  noisy  channel. 


Figure  31  Frequency  offset  synchronization  using  the  ML  algorithm 


This  is  evident  because  the  samples  depicted  in  the  bottom  third  of  Figure 
31  more  closely  resemble  the  original  noisy  samples  without  the  frequency  offset  (top 
third  of  the  figure).  The  ML  derived  frequency  offset  correction  is  applied  to  drive  the 
frequency  offset  ideally  to  exactly  zero.  This  is  never  accomplished  due  to  the  effects  of 
noise.  The  noise  interference  introduces  inaccuracies  in  the  initial  frequency  offset 
estimation  and  a  continual,  or  “on  the  fly”  frequency  correction  algorithm  must  be  used 

for  follow  on  short  training  symbols.  Figure  32  shows  the  resulting  frequency  offset 
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correction  relative  to  sub-carrier  spacing  for  the  short  training  symbols  received  in  Figure 
31  above.  Figure  32  shows  that  the  initial  -40%  frequency  offset  is  detected  and 
corrected  within  the  arrival  of  the  first  two  short  training  symbols.  The  residual 
frequency  offset  correction  is  applied  to  maintain  the  correction  factor  until  completion  of 
the  short  training  symbols. 


Figure  32  ML  Frequency  correction  for  125  kHz  offset  (10  dB  SNR) 


Additionally,  Figure  32  shows,  the  final  frequency  offset  is  less  than  1% 

of  the  sub-carrier  spacing,  or  3125  KHz,  when  SNR  =  10dB.  The  horizontal  axis  in 

Figure  32  depicts  the  frequency  offset  correction  required  based  on  the  current  and 

previous  OFDM  short  training  symbols.  Data  points  in  Figure  32  coincide  with  the 

arrival  of  short  training  symbols  at  which  point  the  frequency  offset  estimate  is  refined 

and  a  new  frequency  correction  is  generated.  The  fidelity  of  the  synchronizer  depends  on 

the  number  of  short  training  symbols.  If  4  or  5  symbols  are  only  detected,  vice  the  full 
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complement,  performance  of  the  synchronizer  suffers.  The  signal  flow  and  structure  of 
this  ML  frequency  synchronizer  is  depicted  below. 


Figure  33  ML  Frequency  Synchronizer 


The  synchronizer  depicted  in  Figure  33  will  gradually  improve  or  maintain 
the  frequency  offset  estimate  with  each  short  training  symbol  received.  Due  to  noise,  a 
final  residual  frequency  offset  may  remain  but  can  be  further  mitigated  using  frequency 
tracking,  or  fine  correction,  by  exploiting  the  long  training  symbols  which  follow  later  in 
the  802.1  lg  preamble.  In  Chapter  V,  the  performance  of  this  synchronizer  at  low  SNR 
and  various  frequency  offsets  will  be  examined. 

3.  Non-linear  Least  Squares  Synchronization 

Another  approach  to  frequency  synchronization  involves  a  Non-linear  Least 
Squares  (NLS)  approximation  of  the  frequency  offset.  This  method  depends  on  the 
known  signal  structure  of  the  preamble.  In  the  case  of  802. 1  lg,  the  preamble  consists  of 
the  10  identical  short  training  symbols  used  for  packet  detection  and  2  identical  long 
training  symbols.  [40] 
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Recall  from  Equation  (4.60)  that 


rm  =  sm  exp(j2xAfmts ),  (4.66) 

where  rm  is  the  received  sample  at  discrete  time  t  =  m,  sm  is  the  equivalent  low  pass 
transmitted  sample,  and  Af  is  the  frequency  offset  caused  by  a  mismatch  between  the 
oscillators  of  the  receiver  and  transmitter  with  all  other  causes  of  frequency  offset 
ignored  for  simplicity.  Because  the  10  short  training  symbols  are  identical  with  16 
sample  periodicity  and  if  we  assume  flat  fading  conditions,  the  received  short  training 
samples  may  be  written  as 

^m+nMs  ^  m+nM  s  [exp  ( j2ntsAfnMs  )] ,  (4.67) 

where  n  is  the  symbol  number  with  the  range  n  =  (0,..., Ns  - 1),  and  in  is  now  the 
sample  index  with  range  m  =  (l,...,Ms).  For  the  short  training  symbols,  Ms  is  the 
duration  of  a  training  symbol  in  samples  and  Ns  is  the  number  of  training  symbols.  Per 
the  802.1  lg  standard,  Ms  =  16  and Ns  =  10.  Equation  (4.67)  can  also  be  applied  to  the 
long  training  symbols,  in  which  case,  NL  =  2  and  ML  =  64.  [40] 

The  transmitted  values  of  the  short  and  long  training  samples  are  known  a  priori 
to  the  receiver  and  are  assumed  to  be  only  affected  by  the  frequency  offset  between 
transmitter  and  receiver.  Therefore,  Equation  (4.67)  can  be  re-written  as 

rm+nM  =  sm  exp  {j2xAfts  [m  +  nMs  ]}  if  1  <  m  <  Ms  and  0  <n<Ns 

(4.68) 

rm+nMs  =  exp \  j2jiAfts  [m  +  nM L  ]}  if  1  <  m  <  Ms  and  0  <  n  <  NL 

for  the  short  and  long  training  samples  respectively.  From  Equation  (4.68),  it  can  be  seen 
that  the  frequency  offset  acts  as  a  complex  exponential  multiplied  by  the  first  training 
symbol.  If  this  can  be  considered  as  a  multiplicative  transfer  function  in  the  frequency 
domain,  then  in  the  time  domain,  this  can  be  treated  as  an  impulse  response  for  the  first 
training  symbol,  i.e.  a  convolution  between  the  frequency  offset  and  the  first  training 
symbol  [10]. 
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Time  domain  convolution  can  be  easily  expressed  as  matrix  multiplication  [8].  In 
this  case,  the  matrix  multiplication  expression  for  Equation  (4.68)  becomes 

4sxi=4sXm/msxi  (4.69) 

where  r  is  the  column  vector  of  received  training  samples  that  constitute  one  symbol,  c  is 
the  column  vector  comprised  of  complex  exponentials  that  represent  the  frequency  offset 
and  s  is  the  circular  convolution  matrix  fonned  from  the  known  training  data.  For  the 
short  training  samples, 
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•'  *3 
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*16 

*15 

*14  " 

•'  *1 

(4.70) 


A  similar  circular  matrix  can  be  constructed  for  the  long  training  samples.  [8]. 

From  Equation  (4.69),  the  column  vector  of  complex  exponentials  caused  by  the 
frequency  offset  can  be  estimated  as 

c=?fr  (4.71) 

where  s 1  is  the  Moore-Penrose  generalized  inverse  of  s  and  represents  the  shortest  least 
squares  solution  to  Equation  (4.69)  [43].  The  circular  convolution  matrix  can  be  stored 
in  the  receiver’s  memory  thereby  reducing  the  quantity  of  computations  done  when  a 
signal  is  detected.  Follow-on,  identical  symbols  can  then  be  leveraged  to  refine  the 
estimate  of  c.  Equation  (4.72)  depicts  such  an  algorithm. 

4=-?t(,I+ . +  fi)  (4-72) 

n 

Equation  (4.72)’s  algorithm  can  be  applied  to  both  short  and  long  training  symbols, 
which  then  may  be  used  to  refine  the  estimate  of  the  vector  of  complex  exponentials.  In 
this  case,  n  =  (l,2,...Ns )  or  n  =  (l,2,...NL)  for  the  short  and  long  training  symbols 

respectively.  Once  c  has  been  accurately  estimated,  the  frequency  offset  can  be 
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computed  by  taking  the  Fast  Fourier  Transfonn  (FFT)  of  c  and  then  the  magnitude  of  the 
frequency  offset  is  found  using  Equation  (4.62). 

a.  Implementation  of  the  NLS  Estimator 

The  NLS  synchronizer  described  above  was  constructed  in  MATLAB. 
The  algorithm  exploits  the  periodicity  of  both  the  short  and  long  training  symbols.  The 
frequency  offset  is  -40  percent  of  sub-carrier  spacing  and  SNR  is  10  dB,  which  are  the 
same  conditions  for  the  ML  frequency  synchronizer.  Once  again,  “on-the-fly”  frequency 
correction  is  implemented  during  the  reception  of  the  short  and  long  training  symbols  to 
evaluate  the  performance  of  the  algorithm.  Figure  34  depicts  the  frequency  offset 
correction  value  that  is  applied  to  subsequent  incoming  training  symbols.  As  expected, 
due  to  noise  in  the  channel,  the  offset  is  never  completely  eliminated,  but  the  algorithm 
routinely  computes  corrections  that  result  in  a  final  frequency  offset  that  is  typically  less 
than  1  percent  of  sub-carrier  spacing  using  both  the  short  and  long  training  symbols. 
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Figure  34  Frequency  offset  correction  using  NLS  algorithm 


In  the  Figure  34  above,  “LTS”  indicates  the  long  training  symbols  per  the 
802.1  lg  standard.  These  two  long  symbols,  with  the  CP  discarded,  constitute  the  final 
correction  to  the  frequency  offset  using  the  NLS  algorithm.  The  CP  is  discarded  to 
simulate  what  an  actual  802.1  lg  receiver  may  have  to  do  when  correcting  for  frequency 
offset  with  a  real  signal.  The  signal  flow  and  structure  of  the  NLS  synchronizer  is 
depicted  in  Figure  35.  It  depicts  the  implementation  of  the  NLS  algorithm  coupled  with 
the  FFT  of  c.  It  is  then  that  the  actual  frequency  offset  is  estimated  using  Equation 
(4.62)  from  section  D2.  This  does  not  limit  the  algorithm  to  a  maximum  frequency  offset 
of  one-half  of  sub-carrier  spacing  apparently.  Using  the  NLS  algorithm,  offsets  of  up  to 
95%  sub-carrier  spacing  were  detected  and  reasonably  corrected.  The  NLS  algorithm 
discussed  not  only  depends  on  a  priori  knowledge  of  the  training  sample  values  but  also 
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assumes  perfect  symbol  timing.  If  symbol  timing  is  unknown,  or  sample  values  are  not 
known  a  priori,  the  NLS  algorithm  may  require  an  iterative  solution,  where  the  receiver 
must  make  an  initial  guess  of  the  frequency  offset.  The  perfonnance  of  the  NLS 
algorithm  will  be  explored  more  thoroughly  in  Section  V  under  various  SNR  coupled 
with  partial  reception  of  the  training  symbols  and  imperfect  symbol  timing. 


A 


Figure  35  NLS  Frequency  Synchronizer 


4.  Summary  of  Frequency  Synchronization  Algorithms 

Analysis  in  Chapter  III  established  the  necessity  of  frequency  synchronization  in 
OFDM  systems.  Many  synchronization  algorithms  exist,  but,  as  discussed  in  this  section 
A3  of  this  chapter,  data-aided  synchronization  seems  to  be  the  most  applicable  for  OFDM 
WLANs  because  synchronization  can  occur  rapidly  with  the  widest  acquisition  range  and 
best  tracking  [12].  The  IEEE  standard  mandates  the  use  of  training  symbols  in  the 
OFDM  preamble  anyway,  discounting  any  argument  against  their  use  in  terms  of 
throughput.  Two  specific  synchronization  algorithms  were  examined  in  this  sub-chapter 
-  maximum  likelihood  (ML)  and  non-linear  least  squares  (NLS).  ML  provided  the 
advantage  of  a  low  complexity  algorithm,  with  rapid  correction  to  a  predetermined 
frequency  offset.  The  ML  did  suffer  from  a  limited  acquisition  range,  but  the  algorithm’s 
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range,  using  either  the  long  or  short  training  symbols,  is  large  enough  to  synchronize  in 
frequency  provided  the  oscillators  used  in  the  transmitter  and  receiver  meet  the  IEEE 
standard.  The  NLS  algorithm,  seemed  to  provide  less  of  a  residual  frequency  offset  than 
the  ML  algorithm  with  much  greater  acquisition  range.  This  shall  be  verified  in  Chapter 
V.  The  NLS  algorithm,  however,  was  more  complex  and  took  a  longer  time  to  achieve 
frequency  synchronization  as  a  comparison  between  Figure  32  (ML)  and  Figure  34 
(NLS)  shows. 
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V.  PERFORMANCE  OF  THE  SYNCHRONIZATION 
ALGORITHMS 


A.  INTRODUCTION 

This  chapter  focuses  on  evaluating  the  perfonnance  of  the  packet  detection  and 
frequency  synchronization  algorithms  discussed  in  Chapter  IV.  Ultimately,  the 
performance  of  any  algorithm  is  dependent  upon  the  nature  of  the  channel  and  strength  of 
the  received  signal  [12].  The  channel  for  all  of  the  Monte  Carlo  simulations  in  this 
chapter  will  be  an  AWGN  channel.  The  discussion  in  Chapter  IV  emphasized  data-aided 
synchronization  and  the  necessity  of  fast  synchronization  in  a  WLAN  environment; 
hence  the  following  perfonnance  evaluations  will  focus  on  the  802.1  lg  16 //sec 
preamble. 

B.  PACKET  DETECTION  PERFORMANCE 

The  two  packet  detection  algorithms  analyzed  in  Chapter  IV  were  called  the 
Matched  Filter  detector  and  the  Delay  and  Correlate  detector.  Both  packet  detection 
algorithms  are  evaluated  in  this  Chapter  by  comparing  their  respective  probability  of 
detection  against  probability  of  false  alann  while  varying  a  single  digit  SNR  (it  is 
assumed  that  both  algorithms  perform  adequately  at  double-digit  SNR).  Such  a 
comparison,  when  depicted  graphically  is  known  as  a  receiver  operating  characteristic 
(ROC)  curve  [37].  To  compute  probability  of  detection,  continuous  short  training 
symbols  are  transmitted  in  an  AWGN  channel.  When  a  packet  detection  algorithm 
computes  a  decision  statistic  that  is  greater  than  the  threshold  value,  a  successful 
detection  is  judged  to  have  occurred.  To  compute  probability  of  false  alarm,  a  noise-only 
signal  is  sent  to  the  receiver,  and  once  again,  the  decision  statistic  is  computed  and  then 
compared  to  the  appropriate  threshold.  When  the  decision  statistic  exceeds  the  threshold 
value,  a  false  detection  is  judged  to  have  occurred.  The  number  of  detections  and  false 
detections  are  tabulated  and  compared  to  the  number  of  decision  statistics  to  determine 
Pd  and  Pfa.  The  threshold  range  is  detennined  by  the  mean  value  of  the  decision  statistic 

which  is  either  fixed  or  dependent  upon  the  evaluated  SNR  values  depending  on  the  type 
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of  detector  being  evaluated.  For  this  perfonnance  evaluation,  because  any  algorithm  is 
expected  to  provide  acceptable  performance  at  medium  to  high  levels  of  SNR,  arbitrary 
values  of  SNR  between  -6  dB  and  +6  dB  have  been  chosen.  Once  ROC  data  has  been 
generated  between  the  Delay  and  Correlate  algorithm  and  the  Matched  Filter  algorithm, 
they  may  be  compared. 

1.  Delay  and  Correlate  Packet  Detection  Performance 

The  threshold  values  used  to  detennine  Pd  and  Pfa  for  various  low  levels  of  SNR 

were  based  on  the  mean  value  of  the  decision  statistic  in  the  test  signal.  Equation  (4.51) 
gave  the  mean  value  of  the  short  training  samples  as  a  function  of  SNR.  The  threshold  is 
then  allowed  to  vary  from  mx  -0.5  to  mx  +0.5,  where  mx  is  the  mean  value  of  the 

decision  statistic  when  a  signal  is  present.  This  threshold  range  is  employed  in  order  to 
detect  the  packet  at  low  SNR.  The  use  of  multiple  samples  within  the  Delay  and 
Correlate  detector  is  used  to  keep  Pfa  at  a  relatively  low  value.  The  table  below  shows 

the  evaluated  SNR  values,  the  mean  value  of  the  decision  statistic  as  computed  from 
Equation  (4.5 1),  and  the  resulting  threshold  range. 


Table  3  Threshold  range  for  evaluated  SNR 


SNR  (dB) 

Mean  value  of  decision  statistic 

K) 

Threshold  range  (yr) 

-6 

.0403 

-.4597  -+  .5403 

-3 

.1115 

-.3885  -+  .6115 

0 

.25 

-.25  -+  .75 

3 

.4437 

-.0563  -+  .9437 

6 

.6388 

-.1388  —+  1.1388 

Once  the  threshold  range  has  been  detennined  for  each  evaluated  SNR,  the 
performance  analysis  is  begun  using  the  test  signal  discussed  earlier.  The  test  signal 
consists  of  16,384  short  training  samples  which,  due  to  the  16  sample  periodicity  and 
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structure  of  the  Delay  and  Correlate  algorithm,  resolves  to  1023  decision  statistics  for 
each  evaluated  SNR.  This  is  because,  in  the  Delay  and  Correlate  algorithm,  the  first 
symbol  is  compared  to  the  second  symbol;  the  second  symbol  is  compared  to  the  third, 
and  so  on.  The  test  signal  has  a  total  of  1024  short  training  symbols,  which  resolves  to 
1023  decision  statistics.  The  test  signal  is  repeated  to  compute  Pfa  but  in  this  case,  the 

test  signal  consists  only  of  noise.  The  resulting  ROC  curve  is  shown  below  for  the 
various  evaluated  SNRs  (shown  as  asterisk  data  points  in  Figure  36).  Schmidl  and  Cox’s 
[38]  Gaussian  approximation  (GA)  of  the  theoretical  performance  of  the  Delay  and 
Correlate  packet  detector  is  depicted  as  solid  lines.  The  curves  are  derived  using  the  GA 
as  it  was  discussed  in  Chapter  0  coupled  with  Equations  (4.36),  (4.43),  (4.51),  and  (4.52). 
From  Figure  36,  the  Gaussian  approximation  (GA)  serves  as  a  lower  bound  at  higher 
tested  SNRs  while  it  closely  approximates  the  computed  values  at  low  SNR.  For  clarity, 
the  portion  of  the  ROC  curves  for  which  Pfa  <  .  1 5  is  depicted. 


Figure  36  ROC  curves  for  Delay  and  Correlate  Detector 
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While  the  IEEE  standard  makes  no  mention  of  probability  of  false  alarm 
requirements,  it  is  certainly  possible  that  designers  would  want  to  minimize  the 
probability  of  false  alarm  because  the  receiver  may  be  processing  a  false  alann  and  miss 
an  actual  packet.  Furthermore,  I  speculate  that  this  increased  signal  processing  might 
result  in  a  noticeable  degradation  in  the  time  between  subscriber  unit  battery  recharges, 
particularly  for  small  handheld  devices.  The  IEEE  standard  does,  however,  mandate  a 
maximum  packet  error  rate  of  10%  for  a  PDSU  length  of  1000  bytes  given  the  mandated 
receiver  sensitivity  for  the  various  data  rates  [5].  If  an  arbitrary  Pfa  of  .05  is  acceptable, 

then  Figure  36  shows  that  the  Delay  and  Correlate  detector,  using  16  sample  periodicity, 
in  an  AWGN  channel,  meets  the  IEEE  packet  error  rate  requirement  when  SNR  is 
evaluated  at  3  and  6  dB  SNR.  This  assumes  that  the  packet  cannot  be  demodulated 
properly,  and  hence  in  error,  if  it  is  never  detected.  Therefore,  for  the  purposes  of  this 
performance  analysis,  an  acceptable  baseline  Pd  is  .9. 

a.  Improved  Performance  from  the  Delay  and  Correlate  Detector 

The  performance  of  the  Delay  and  Correlate  detector  can  be  further 
improved  if  additional  samples  are  used  to  detect  the  packet.  This  should,  as  discussed  in 
Chapter  IV,  allow  for  a  low  threshold  while  minimizing  false  alarms.  The  short  training 
symbols,  at  the  start  of  the  802.1  lg  packet  consist  of  10  symbols,  each  16  samples  long. 
If  the  802.1  lg  receiver  contained  a  buffer  that  could  store  more  of  the  short  training 
samples,  then  these  additional  samples  could  be  used  in  the  Delay  and  Correlate  detector. 
The  buffered  sample  values  could  also  be  used  for  other  synchronization  tasks  as  well. 
Such  a  detector,  however,  would  cost  more.  Along  with  greater  receiver  complexity  and 
additional  signal  processing,  additional  hardware  may  be  required  (e.g.  a  high  gain 
antenna)  to  successfully  demodulate  the  packet  at  very  low  SNR,  thereby  increasing  unit 
cost.  Such  a  specialized  detector’s  performance,  using  the  Delay  and  Correlate 
algorithm,  is  well  illustrated  in  the  table  below.  The  Delay  and  Correlate  algorithm  and 
the  simulation  parameters  remain  the  same  as  in  the  previous  perfonnance  analysis, 
except  that  the  size  of  the  delay  block  has  been  changed  to  80  samples  from  16  samples. 
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Although  there  are  160  samples  in  the  short  training  symbols,  because  the  first  half  of  the 
preamble  is  compared  to  the  second  half,  the  algorithm  is  limited  to  80  repetitive 
samples. 


Table  4:  Pd  and  Pfa  for  80  sample  Delay  and  Correlate  Detector 


SNR  (dB) 

Pd 

/  Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

/Pfa 

Pd 

Pfa 

-6 

0.98 

0.77 

0.96 

0.65 

0.95 

0.52 

0.95 

0.47 

0.92 

0.39 

0.91 

0.35 

0.91 

0.30 

0.90 

0.28 

-3 

1.00 

0.40 

1.00 

0.25 

0.99 

0.17 

0.99 

0.13 

0.99 

0.06 

0.98 

0.05 

0.96 

0.02 

0.94 

0.01 

0.92 

0.01 

0 

1.00 

0.15 

1.00 

0.06 

1.00 

0.02 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

0.99 

0.00 

0.98 

0.00 

3 

1.00 

0.03 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

6 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

yT 


In  Table  4,  Pd  and  Pfa  are  shown  as  a  function  of  SNR  and  threshold.  For 

this  performance  analysis  each  run  is  conducted  at  a  tested  SNR,  which  is  adjusted  by 
increasing  the  noise  variance.  The  threshold,  yT,  is  gradually  increased  for  the  SNR 
values  listed  in  the  left  most  column  of  the  table.  As  discussed  in  Chapter  0,  the 
threshold  can  be  varied  to  achieve  a  desired  Pd.  yT  is  set  at  an  initial  value  of  my  / 1 0 

and  allowed  to  gradually  increase  to  a  value  of  my  j2.  In  order  to  present  the  table  data 

clearly,  the  probabilities  have  been  rounded  to  2-decimal  accuracy.  Table  values  that 
have  a  red  background  indicated  a  failure  to  meet  the  previously  discussed  standard  of 
Pd=. 9.  Cells  with  a  yellow  background  show  where  P,  >  .9  but  Pfa  >  .05.  A  blue 

background  indicates  the  desired  criteria  of  Pd  >  .9  and  Pfa  <  .05  has  been  met.  From 

Table  4,  the  performance  of  the  80  sample  Delay  and  Correlate  detector  clearly  exceeds 
the  performance  of  the  same  algorithm  using  only  16  samples.  As  discussed  earlier,  the 
16  sample  detector  met  the  desired  criteria  for  +3  dB  and  +6  dB  SNR.  The  80  sample 
detector  also  met  that  criteria  and  meets  the  criteria  at  0  dB  and  -3  dB  for  higher 
threshold  levels.  Additionally  at  0  dB,  the  80  sample  detector  missed  the  arbitrarily  set 
Pfa  by  only  a  small  amount  at  low  thresholds. 
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2. 


Matched  Filter  Packet  Detection  Performance 


The  other  packet  detection  algorithm  discussed  in  Chapter  IV  was  the  Matched 
Filter  detector.  Because  the  values  of  the  short  training  samples  are  known  a  priori  to  the 
receiver,  this  could  be  used  to  facilitate  Matched  Filter  detection.  For  the  performance 
analysis,  the  same  test  signal  used  in  section  B1  was  also  used  to  evaluate  the  Matched 
Filter  detector’s  performance.  The  detector,  which  has  the  same  signal  structure  as 
depicted  in  Figure  28,  generates  16368  decision  statistics  from  the  test  signal.  The  value 
of  the  decision  statistic,  is  dependent  upon  the  tested  SNR  but  closely  matches  the 

bottom  half  of  Figure  29.  As  Figure  29  shows,  the  values  of  xM  “spike”  upwards  when 
there  is  a  perfect  correlation  between  an  incoming  short  training  symbol  and  the  Matched 
Filter  coefficients.  The  “spikes”  occur  every  16  samples  (the  other  values  of 
XM  represent  partial  or  no  correlation  with  the  filter  coefficients).  These  particular  values 

of  Xm  are  the  focus  of  the  performance  analysis  for  the  Matched  Filter  detector.  Hence, 
there  are  actually  1023  values  of  jMthat  can  be  used  for  detecting  a  short  training 

symbol  in  noise.  The  performance  analysis  was  run  using  the  same  tested  values  of  SNR 
that  were  used  in  the  evaluation  of  the  Delay  and  Correlate  detector.  For  the  analysis,  the 
threshold  was  allowed  to  vary  in  accordance  with  the  mean  value  of  the  decision  statistic 
for  each  tested  SNR.  For  the  16  sample  Matched  Filter  detector,  the  ROC  data  is  best 
illustrated  in  the  table  below.  To  analyze  Pfa,  the  test  signal  consists  of  only  noise. 

Table  5:  Pd  and  Pfa  for  Matched  Filter  Detector 


yT - ► 
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In  Table  5,  the  left  most  columns  represent  low  threshold  and  the  right  most 
columns  represent  high  threshold.  For  the  analysis,  yT  =  mx/%  to  mx/\.5.  The 

performance  of  the  16  sample  Matched  Filter  detector  shows  good  performance  at  higher 
tested  SNRs  but  fails  to  meet  the  arbitrary  standard  of  Pfa  <  .05  at  low  SNR  for  the 

various  thresholds.  The  Matched  Filter  Detector  also  fails  to  meet  the  arbitrary  standard 
of  Pd  >  .9  and  Pfa  <  .05  only  at  the  lowest  tested  SNR  and  the  highest  values  of  yT.  The 

performance  of  this  detector,  like  the  Delay  and  Correlate  detector  can  be  improved  by 
using  more  samples.  The  drawback  will  be  greater  complexity  of  the  detector. 

a.  Improved  Performance  from  the  Matched  Filter  Detector 

An  improved  Matched  Filter  detector  can  be  built  by  using  more  samples. 
The  improved  detector  would  require  an  additional  buffer  or  delay  blocks  to  store  the 
samples  so  they  can  be  compared  against  an  increased  number  of  filter  coefficients.  A 
modified  structure  is  shown  in  Figure  37.  This  structure  is  specific  for  exploiting  the 
repetitive  short  training  symbols.  The  blocks  annotated  as  z  D  represent  single  sample 
delay  blocks  while  the  block  annotated  as  z  M  represents  a  16  sample  delay  block.  This 
detector  represents  a  significant  increase  in  complexity  relative  to  the  Matched  Filter 
detector  analyzed  previously.  For  comparison,  the  original  detector’s  block  diagram  is 
depicted  in  Figure  28. 
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Figure  37  32  sample  Matched  Filter  Detector 


The  resulting  performance  analysis  of  the  new  or  “double”  Matched  Filter 
detector  shows  that  the  use  of  more  samples  allows  for  a  higher  threshold,  thereby 
reducing  the  false  alarm  rate.  For  the  performance  analysis,  the  new  mean  of  the  decision 
statistic  becomes  twice  the  old  mean,  or  mx  =  A.  The  threshold,  yT,  is  allowed  to  vary 

from  mxj 8  to  mxj  1.5.  The  results  of  the  performance  analysis  are  shown  in  Table  6. 


Table  6:  Pd  and  Pfa  for  Improved  Matched  Filter  Detector 


yT  - ► 
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3.  Summary  of  Packet  Detection  Performance 

Sub-Sections  B1  and  B2  discussed  the  perfonnance  analysis  of  the  Delay  and 
Correlate  detector  and  the  Matched  Filter  detector  respectively.  It  is  shown  that  the 
Delay  and  Correlate  Detector  had  difficulty  meeting  the  arbitrary  standard  of 
Pd  >  .9  and  Pfa  <  .05  at  low  SNR.  Although  some  improvement  was  experienced  by  using 

additional  samples  (with  the  drawback  of  likely  greater  cost  and  complexity), 
performance  still  lagged  at  very  low  SNR.  The  “improved”  Delay  and  Correlate  detector 
was  dependent  upon  receiving  all  ten  short  training  samples  because  it  incorporated  a 
delay  buffer  equivalent  to  80  samples.  Although  it  is  beyond  the  scope  of  this  thesis,  the 
performance  of  the  80  sample  detector  could  be  evaluated  with  only  partial  receipt  of  the 
ten  short  training  samples.  It  is  expected  that  the  perfonnance  of  the  detector  would 
further  deteriorate. 

The  Matched  Filter  detector’s  performance  is  initially  analyzed  using  16  filter 
coefficients  due  to  the  16  sample  periodicity  of  the  short  training  symbols.  The 
performance  of  the  detector  is  comparable  to  the  Delay  and  Correlate  Detector  using  80 
samples.  This  is  graphically  evident  in  the  similar  color  schemes  of  Table  4  and  Table  5. 
To  determine  if  the  desired  values  of  Pd  and  Pfa  could  be  obtained,  additional  samples 

were  used  to  create  more  filter  coefficients  and  then  placed  in  series  with  the  single 
Matched  Filter.  In  this  case,  the  receipt  of  at  least  two  consecutive  short  training  symbols 
would  be  required  for  detection.  The  resulting  Double  Matched  Filter  detector,  as 
depicted  in  Figure  37,  would  be  more  complex  and  likely  more  costly  than  a  single 
Matched  Filter  detector,  but  met  the  performance  requirements  across  the  tested  range  of 
SNR  for  middle  to  high  threshold  values. 

In  conclusion,  the  Double  Matched  Filter  detector  provided  the  best  perfonnance 
at  the  tested  SNR  range  but  meeting  the  performance  requirements  meant  using 
additional  filter  coefficients,  leading  to  greater  detector  complexity.  This  result  is  not 
completely  surprising.  The  Matched  Filter  is  the  best  filter  for  maximizing  SNR  [10]. 
Intuitively,  it  makes  sense  that  the  same  filter  would  provide  the  best  performance  at  a 
given  SNR.  For  a  commercial  vendor,  such  a  detector  may  not  be  economically  viable 


99 


due  to  the  requirement  to  minimize  subscriber  unit  production  costs.  On  the  other  hand, 
such  high  performance  detectors  may  find  specialized  uses,  particularly  when  coupled 
with  hardware  solutions,  such  as  high  gain  antennas. 

C.  FREQUENCY  SYNCHRONIZATION  PERFORMANCE 

In  Chapter  IV,  the  Maximum  Likelihood  (ML)  and  Non-linear  Least  Squares 
(NLS)  frequency  synchronization  techniques  were  introduced.  This  section  focuses  on 
the  performance  analysis  and  comparison  of  those  techniques.  As  discussed  previously, 
the  necessity  to  detect  and  correct  any  frequency  offset  error  suggests  leveraging  the 
short  and  long  training  symbols  that  are  placed  at  the  beginning  of  the  802.1  lg  preamble 
in  accordance  with  the  IEEE  standard.  For  the  performance  evaluation,  the  test  signal 
consists  of  the  ten  short  training  symbols  sent  to  the  receiver  in  an  AWGN  channel.  As 
in  the  detection  analysis,  multi-path  effects  are  ignored  for  the  purpose  of  this  thesis.  The 
test  signal  is  sent  multiple  times  to  the  receiver.  In  each  case,  the  tested  algorithm  must 
correct  the  frequency  offset.  Although  in  practice,  there  will  always  be  some  residual 
frequency  offset  due  to  phase  noise,  it  is  not  considered  in  this  evaluation.  Hence,  the 
final  frequency  offset  should  be  zero  in  theory  after  the  application  of  the  frequency 
correction  algorithm.  The  results  of  each  algorithm’s  multiple  frequency  correction 
attempts  are  displayed  as  histograms  in  the  following  sub-sections. 

1.  Maximum  Likelihood  Frequency  Synchronization  Performance 

The  perfonnance  analysis  of  the  ML  algorithm  consists  of  the  test  signal,  depicted 
in  Figure  31,  corrupted  by  noise  and  a  significant  frequency  offset  (40  percent  of  sub¬ 
carrier  spacing).  In  each  case,  the  ML  is  able  to  detect  and  correct  the  frequency  offset, 
using  only  the  short  training  symbols.  For  the  initial  analysis,  all  ten  short  training 
symbols  are  used  to  refine  the  frequency  estimate.  Additionally,  the  ML  algorithm  is 
examined  for  robustness  by  using  only  4  short  training  symbols  and  2  long  training 
symbols.  These  results  will  be  compared  to  the  performance  data  when  the  algorithm 
uses  all  ten  symbols.  The  expectation  is  that  the  ML  algorithm  will  be  relatively  immune 
to  the  reduction  in  the  number  of  short  training  symbols  available  for  synchronization. 
This  is  suggested  from  the  implementation  results  displayed  in  Figure  32.  If  less  training 
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symbols  are  required  for  frequency  synchronization,  then  the  initial  training  symbols  may 
be  discarded  after  being  used  for  packet  detection,  perhaps  reducing  complexity  and  cost 
of  the  receiver.  The  results,  for  500  trials  at  10  dB  SNR  and  6  dB  SNR,  are  shown  in  the 
histograms  of  Figure  38. 


Figure  38  Remaining  Frequency  offset  using  ML  algorithm  (10  and  6  dB  SNR) 

using  10  short  training  symbols 

The  performance  analysis  results  depicted  in  the  histograms  suggest  that  the  ML 
algorithm  is  an  accurate  technique  to  synchronize  in  frequency.  The  average  final 
frequency  offset,  after  10  short  training  symbols,  was  less  than  1  percent  frequency  offset 
at  10  dB  SNR.  For  all  500  trials,  the  final  frequency  offset  never  exceeds  two  percent  of 
sub-carrier  spacing  for  the  histogram  on  the  left  despite  an  initial  offset  of  40  percent  of 
sub-carrier  spacing  (or  about  125  kHz).  As  noted  previously,  the  maximum  allowable 
frequency  offset  for  the  ML  algorithm  is  50  percent  relative  to  sub-carrier  spacing.  As 
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expected,  the  performance  of  the  ML  algorithm  degrades  as  SNR  is  reduced.  This  is 
evident  from  the  greater  variance  in  the  final  frequency  offset  percentage  visibly  present 
in  the  histogram  on  the  right  side  of  Figure  38.  Additionally,  there  are  far  fewer  trials 
that  resulted  in  a  zero  or  near  zero  offset  condition. 

Furthermore,  as  SNR  is  reduced,  perfonnance  continues  to  degrade.  The 
degradation  is  graphically  illustrated  by  the  histogram  depicting  final  frequency  offset  at 
3  dB  SNR  in  Figure  39.  As  shown  in  Figure  39,  the  ML  algorithm  may  also  employ  the 
long  training  symbols  in  addition  to  the  short  training  symbols  for  frequency 
synchronization.  The  purpose  of  showing  the  use  of  the  long  training  symbols  in 
frequency  synchronization  on  the  left  hand  side  of  Figure  39  is  to  highlight  the 
performance  advantage  inherent  in  the  use  of  the  long  training  symbols  to  “fine”  tune  the 
frequency  offset  as  recommended  by  the  IEEE  standard  [4],  This  is  clearly  evident  when 
the  two  histograms  in  Figure  39  are  compared.  The  histogram  on  the  left  hand  side  of 
Figure  39  depicts  ML  algorithm  performance  without  the  long  training  symbols.  Another 
reason  is  because  the  next  section  will  examine  the  “robustness”  of  the  ML  algorithm. 
The  robustness  of  the  algorithm  is  analyzed  by  limiting  the  number  of  training  symbols 
available  to  the  synchronization  algorithm.  This  performance  analysis  involves 
simulating  the  reception  of  only  a  few  of  the  short  training  symbols  and  the  long  training 
symbols  vice  the  full  complement. 
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Figure  39  Remaining  frequency  offset  percentage  using  ML  algorithm  at  3  dB 

SNR 


a.  Robustness  of  the  ML  Frequency  Correction  Algorithm 

To  test  the  robustness  of  the  ML  algorithm,  the  analysis  was  re-run  using 
only  four  short  training  symbols  and  then  again  with  four  short  training  symbols  and  the 
two  long  training  symbols  to  highlight  the  critical  role  the  long  training  symbols  play  in 
synchronization.  In  each  case,  the  new  test  signal  was  re-sent  500  times  to  the  receiver  at 
10  dB  SNR.  The  results,  displayed  in  Figure  40,  can  be  contrasted  with  the  left  hand  side 
of  Figure  38.  While  both  Figures  are  histograms  that  depict  the  performance  of  the  ML 
algorithm  at  10  dB  SNR,  the  results  in  Figure  38  are  based  on  the  reception  of  10  short 
training  symbols  and  the  two  long  training  symbols.  The  results  in  Figure  40  show  only  a 
slight  degradation  in  performance.  This  robustness  of  the  ML  algorithm  is  not  surprising 
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because  of  the  ability  of  the  ML  algorithm  to  quickly  detect  and  correct  the  frequency 
offset  as  suggested  from  the  implementation  results  depicted  in  Figure  32. 
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Figure  40  Remaining  Frequency  offset  percentage  using  ML  algorithm  and 

reduced  training  symbols 


2.  Non-linear  Least  Squares  Frequency  Synchronization  Performance 

The  second  frequency  technique  that  was  discussed  in  Section  IV  was  Non-linear 
Least  Squares  (NLS).  Although  this  technique  could  be  used  in  a  non  data-aided 
synchronization  scheme  as  discussed  in  Chapter  IV,  the  necessity  in  the  WLAN 
environment  to  synchronize  in  frequency  required  focusing  on  a  data-aided  scheme,  using 
both  the  short  and  long  training  symbols.  This  will  also  allow  for  a  straight  forward 
comparative  analysis  with  the  ML  technique.  The  implementation  results  depicted  in 
Figure  34  show  that  the  NLS  technique  could  synchronize  in  frequency,  but  the  long 
training  symbols  are  needed  to  approach  the  same  performance  standard  as  that  of  the 
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ML  algorithm.  I  judge  this  to  be  acceptable  because  the  IEEE  standard  specifically 
recommends  using  the  two  long  training  symbols  for  fine  frequency  correction  [4].  The 
performance  of  the  NLS  algorithm  is  evaluated  using  the  same  conditions  that  governed 
the  analysis  of  the  ML  technique  -  large  initial  frequency  offset,  low  SNR,  multi-path 
effects  disregarded,  and  eventually  examining  the  robustness  of  the  NLS  algorithm  by 
using  only  4  short  training  symbols.  The  initial  frequency  offset  is  -40%  of  sub-carrier 
spacing  (or  125  kHz  less  than  the  sub-carrier’s  center  frequency).  In  the  last  condition 
(robustness),  it  is  expected  that  the  NLS  algorithm  will  lag  in  performance  relative  to  the 
ML  algorithm.  To  graphically  depict  the  perfonnance  of  the  NLS  technique,  histograms 
are  produced  (as  in  the  ML  algorithm  analysis)  by  using  a  test  signal  re-sent  to  the 
receiver  500  times.  In  this  case,  the  test  signal  consists  of  the  10  short  training  symbols, 
the  cyclic  prefix  of  the  long  training  symbols  (which  is  discarded  by  the  receiver)  and  the 
two  long  training  symbols  themselves.  The  first  of  these  histograms,  in  Ligure  41,  shows 
the  remaining  frequency  offset  percentage  relative  to  sub-carrier  spacing  at  10  dB  SNR 
(left  hand  side)  and  6  dB  SNR  (right  hand  side). 
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Figure  41  Remaining  Frequency  offset  percentage  using  NLS  algorithm  using  10 

short  and  2  long  training  symbols 

From  Figure  41,  it  is  evident  that  the  NLS  algorithm  effectively  performs  the 
frequency  synchronization  task.  All  of  the  trials  resulted  in  a  frequency  offset  of  less 
than  two  percent  for  both  10  dB  and  6  dB  SNR.  As  expected,  as  SNR  is  reduced,  the 
performance  of  the  NLS  is  degraded.  Results  for  6  dB  SNR  are  depicted  as  the  right 
histogram  in  Figure  4 1 .  The  histogram  shows  a  similar  performance  degradation  that  was 
apparent  for  the  ML  algorithm  -  greater  variance  in  the  final  frequency  offset  percentage 
(relative  to  sub-carrier  spacing)  and  less  trials  that  results  in  zero  or  near  zero  remaining 
frequency  offset.  Similar  decreasing  performance  can  be  observed  when  the  SNR  is 
reduced  to  3  dB.  That  histogram  will  be  highlighted  in  the  final  section  of  this  chapter 
when  a  direct  comparison  with  the  ML  algorithm  is  made. 
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a.  Robustness  of  the  NLS  Frequency  Correction  Algorithm 

The  performance  results  shown  in  the  NLS  histograms  suggest  similar 
performance  relative  to  the  ML  algorithm.  The  NLS  histograms  presented  so  far  depict 
performance  of  the  algorithm  when  all  ten  short  training  symbols  and  two  long  training 
symbols  are  used  for  frequency  synchronization.  To  examine  the  robustness  of  the 
algorithm,  additional  analysis  is  performed  under  the  same  conditions  as  the  ML 
algorithm.  Once  again,  the  test  signal  consists  of  four  short  training  symbols  and  the  two 
long  training  symbols.  To  further  test  the  robustness  of  the  algorithm,  a  second  test 
signal  consisting  only  4  short  training  symbols  is  used.  The  results  for  500  trials  are 
shown  in  Figure  42. 
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The  histograms  in  Figure  42  reinforce  the  power  of  the  long  training 
symbols  in  frequency  synchronization.  This  is  especially  true  for  the  NLS  algorithm,  as  a 
64  x  64  matrix  is  used  to  for  “fine”  frequency  adjustment.  In  this  case,  fine 
synchronization  using  the  long  training  symbols  quickly  eliminates  any  residual 
frequency  offset.  The  Figure  42  histograms  can  also  be  viewed,  conversely,  as  showing 
that  the  NLS  algorithm  is  highly  dependent  upon  the  long  training  symbols  for  frequency 
synchronization.  Unlike  the  ML  algorithm,  the  NLS  algorithm  slowly  synchronizes  in 
frequency  with  each  received  short  training  symbol.  This  result  is  not  unexpected  and 
was  evident  in  the  implementation  results  shown  in  Figure  34.  Hence,  the  NLS 
algorithm  is  robust  only  if  the  long  training  symbols  are  received  and  used  for  frequency 
synchronization.  The  algorithm  will  not  sufficiently  synchronize  the  receiver  using  only 
4  short  training  symbols. 

3.  Summary  of  Frequency  Synchronization  Performance 

Two  types  of  frequency  synchronization  algorithms  were  evaluated  at  low  SNR  - 
the  Maximum  Likelihood  algorithm  and  the  Non-linear  Least  Squares  algorithm. 
Histograms  were  produced  from  multiple  test  signals  sent  to  receivers  using  the  two 
algorithms.  Both  algorithms  performed  effectively.  The  robustness  of  each  algorithm 
was  also  examined.  Here,  the  performance  of  the  algorithms  contrasted  sharply.  The  ML 
algorithm  was  able  to  quickly  eliminate  most  of  the  frequency  offset  using  only  the  short 
training  symbols  (only  4  symbols  in  the  tested  case).  The  long  training  symbols  were  an 
important  facet  of  the  synchronization  scheme  in  the  ML  algorithm  but  the  algorithm 
could  correct  a  frequency  offset  without  them  if  required.  The  NLS  algorithm  was  highly 
dependent  upon  the  use  of  the  long  training  symbols  to  eliminate  the  frequency  offset. 
Without  the  long  training  symbols,  the  offset  could  not  be  eliminated  with  only  4  short 
training  symbols.  As  discussed  previously,  this  was  because  the  NLS  algorithm  only 
slowly  corrects  the  frequency  offset  as  opposed  to  the  ML  algorithm,  which  quickly 
eliminates  the  offset.  Although  the  NLS  scheme  was  less  robust,  due  to  the 
dependency  of  the  long  training  symbols,  the  performance  of  the  NLS  algorithm  was 
nevertheless  superior  to  the  ML  algorithm.  This  is  evident  in  a  final  comparison  of 
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both  synchronization  algorithms  at  very  low  SNR  (3  dB),  using  only  4  short  training 
symbols  and  the  long  training  symbols  as  shown  in  Figure  43. 


Frequency  offset  percentage  relative  to  sub-carrier  spacing 


Figure  43  Comparison  of  remaining  Frequency  offset  percentage  at  3  dB  SNR 
with  4  short  and  2  long  training  symbols 

In  summary,  this  chapter  analyzed  the  synchronization  algrorithms  discussed  in 
the  Chapter  IV.  The  various  algorithms  for  packet  detection  and  frequency 
synchronization  were  comparatively  analyzed.  The  following  conclusion  chapter  will 
summarize  the  results  of  this  performance  analysis. 
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VI.  CONCLUSION 


This  thesis  addresses  promising  synchronization  algorithms  that  may  be  used  to 
extend  the  range  of  802.1  lg  WLANs.  The  scope  of  the  thesis  centers  on  the  analysis  of 
packet  detection  (frame  synchronization)  and  frequency  synchronization  algorithms. 
Furthermore,  to  focus  on  the  implementation  of  those  algorithms,  only  performance  in  an 
AWGN  channel  is  considered. 

Frame  synchronization  algorithms  that  seem  to  hold  promise  are  the  Delay  and 
Correlate  packet  detector  and  the  Matched  Filter  Detector.  Promising  frequency 
synchronization  algorithms  highlighted  in  the  thesis  were  the  Maximum  Likelihood 
algorithm  and  the  Non-Linear  Least  Squares  algorithm.  The  implementation  of  all  of 
these  algorithms  were  successfully  completed  in  MATLAB  and  discussed  in  Chapter  IV. 
The  MATLAB  code  is  included  as  Appendix  A.  The  performance  analysis  of  these 
algorithms,  under  test  conditions  of  low  SNR,  is  covered  in  detail  in  Chapter  V.  The 
results  are  discussed  in  the  following  section. 

A.  RESULTS  OF  THE  PERFORMANCE  ANALYSIS 

The  performance  analysis  of  Chapter  V  focuses  on  two  specific  packet  detection 
synchronization  algorithms  using  data-aided  synchronization  at  low  SNR.  Specifically, 
the  packet  detection  algorithms  are  evaluated  in  Chapter  V  with  an  SNR  range  of 
+6  dB  to  -6dB,  using  the  short  training  symbols  at  the  beginning  of  the  802. 1 1  g 
preamble.  The  long  training  symbols  are  not  used,  recognizing  their  necessity  for  the 
follow-on  task  of  frequency  synchronization.  While  Receiver  Operating  Characteristic 
curves  of  Pd  vs  Pfa  are  generated,  the  results  for  both  algorithms  are  best  displayed  in 

color-coded  table  format.  The  results  for  the  80-sample  Delay  and  Correlate  detector, 
displayed  in  Table  4,  showed  that  the  detector  could  not  meet  the  arbitrary  standard  of 
Pd  >  .9  and  Pfa  <  .05  at  lowest  tested  SNR  (-6  dB)  regardless  of  the  threshold  setting.  Of 

note,  the  Delay  and  Correlate  detector  suffered  from  a  very  high  probability  of  false 

alarm  at  the  lowest  tested  SNR.  This  was  the  most  robust  of  the  Delay  and  Correlate 

detectors  implemented  in  this  thesis.  The  80  sample  Delay  and  Correlate  detector  uses 
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the  most  samples  available,  assuming  full  reception  of  all  ten  short  training  symbols. 
Although  beyond  the  scope  of  this  thesis,  it  is  likely  that  fewer  received  symbols  would 
result  in  a  degradation  of  performance  for  this  detector. 

The  Matched  Filter  detector’s  performance  was  similar  to  that  of  the  Delay  and 
Correlate  detector  across  the  tested  SNR  range.  There  was  one  difference:  Pfa  was 

consistently  lower  than  that  of  the  80-sample  Delay  and  Correlate  detector  at  the  lowest 
tested  SNR.  Despite  that,  the  Matched  Filter’s  Pfa  was  still  above  the  arbitrary  standard 

at  the  lowest  tested  SNR.  As  discussed  in  Chapter  IV,  the  use  of  more  samples  within  a 
detector  algorithm  would  allow  for  a  lower  threshold  while  maintaining  the  same  false 
alarm  rate  or  would  lower  the  false  alarm  rate  for  a  given  threshold.  This  led  to  an 
implementation  and  evaluation  of  a  32  sample  (or  double)  Matched  Filter  detector.  The 
results,  shown  in  Table  6,  highlight  the  superior  performance  of  this  detector  throughout 
the  range  of  tested  SNR  values.  The  32  sample  Matched  Filter  detector  outshines  the 
maximum  sample  (80)  Delay  and  Correlate  detector.  Furthermore,  its  superior 
performance  is  enhanced  because  it  requires  the  reception  of  only  two  vice  10  short 
training  symbols,  which  the  80-sample  Delay  and  Correlate  detector  is  dependent  upon. 

The  performance  results  of  the  frequency  synchronization  algorithms  are  best 
displayed  as  histograms.  Once  again,  the  emphasis  is  on  data-aided  synchronization.  In 
this  case,  both  short  and  long  training  symbols  are  employed.  Each  frequency 
synchronization  algorithm  was  initially  tested  using  10  short  training  symbols  and  the 
two  long  training  symbols.  Although  the  effect  of  multi-path  was  beyond  the  scope  of 
this  thesis,  the  cyclic  prefix  of  the  first  long  training  symbol  was  discarded  by  the 
receiver  to  more  realistically  simulate  the  performance  of  these  algorithms  in  a  WLAN 
environment.  After  initial  testing,  the  algorithms  were  tested  for  “robustness”  -  in  which 
the  number  of  available  short  training  symbols  was  reduced  from  10  to  4.  This  was  done 
to  simulate  that  short  training  symbols  used  for  packet  detection  are  not  buffered  by  the 
receiver  and  therefore  not  available  for  frequency  synchronization.  Although  it  is 
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certainly  technically  possible  for  an  OFDM  receiver  to  buffer  the  short  training  symbols, 
the  receiver  may  not  buffer  the  samples  because,  as  discussed  in  Chapter  III,  it  may  not 
be  economically  viable  to  do  so. 

The  histograms  in  Chapter  V  show  that  both  the  Maximum  Likelihood  and  Non¬ 
linear  Least  Squares  algorithms  adequately  correct  a  large  frequency  offset  at  10  dB 
SNR.  The  histograms  for  both  algorithms  show  very  similar  perfonnance,  with  both 
algorithms  suffering  performance  degradation  as  SNR  is  reduced  to  the  minimum  tested 
SNR  (3  dB).  The  most  telling  results  are  contained  in  Figure  43.  The  histograms 

depicted  there  provide  a  direct  comparison  between  the  ML  and  NLS  frequency 
synchronization  algorithms  at  the  lowest  tested  SNR  coupled  with  a  reduced  number  of 
training  symbols  available  to  correct  a  frequency  offset.  The  NLS  algorithm  offers 
superior  performance  relative  to  the  ML  algorithm.  The  disadvantage  of  the  NLS 
algorithm  is  its  complexity.  Use  of  the  NLS  algorithm  requires  that  the  receiver  store  a 
16  x  16  circular  convolution  matrix  for  short  training  symbols.  For  the  long  training 
symbols,  the  receiver  must  store  a  64  x  64  matrix!  While  the  NLS  algorithm’s 
performance  was  superior,  such  complexity  may  not  be  economically  viable  for  a 
commercial-off-the-shelf  (COTS)  receiver. 

B.  RECOMMENDATIONS  FOR  FOLLOW  ON  RESEARCH 

The  use  of  OFDM  as  a  viable  modulation  technique  to  achieve  high  data  rates, 
while  mitigating  ISI,  suggests  its  continued  proliferation  within  WLANs.  This  thesis 
highlighted  several  promising  synchronization  techniques  that  successfully  detected  and 
frequency  synchronized  the  802.1  lg  packet  using  the  IEEE  mandated  training  symbols  at 
low  SNR.  No  direct  comparison  was  made  to  any  existing  commercial  synchronization 
techniques  however.  The  software  analysis  of  the  promising  techniques  discussed  in 
Chapter  V  naturally  suggests  a  follow  on  comparison  between  those  techniques  and  the 
synchronization  algorithms  currently  in  use  in  commercial  products  to  verify  the 
improved  performance  and  to  accurately  assess  tradeoffs  between  performance  and  cost. 
Due  to  the  proprietary  nature  of  commercial  synchronization  algorithms,  this  could  be 
done  using  COTS  equipment. 
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Also,  because  multiple  synchronization  tasks  are  required  to  successfully 
demodulate  OFDM,  a  wide  variety  of  synchronization  algorithms  have  been  developed. 
This  wide  variety  of  algorithms  and  tasks  cannot  be  adequately  addressed  within  the 
limited  scope  of  a  single  thesis.  This  offers  three  additional  areas  for  follow  on  research. 

Although  phase  noise  correction  and  symbol  timing  were  reviewed  in  Chapter  III, 
analysis  of  synchronization  algorithms  that  address  those  issues  offer  opportunity  for 
follow-on  research.  Next,  because  the  IEEE  standard  does  not  specify  the  types  of 
synchronization  algorithms  to  be  used  for  the  various  synchronization  tasks,  new 
algorithms  are  being  developed  that  require  comparative  analysis.  Finally,  the 
performance  analysis  of  the  evaluated  synchronization  algorithms  was  conducted  in  a 
simulated  AWGN  channel.  While  incorporating  multi-path  into  the  performance  analysis 
would  be  technically  more  difficult,  a  more  realistic  simulation  could  be  performed 
because  WLANs  often  operate  in  a  multi-path  environment. 
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APPENDIX 


The  following  MATLAB  code  implements  and  provides  output  for  the  various 
detectors  and  synchronizers  used  throughout  this  thesis.  Use  of  ‘functions’  is  employed, 
allowing  the  recycling  of  usable  code. 


delay_correlate_detector.m  -  This  m-filc  implements  the  Delay  and  Correlate  detector 
that  was  extensively  analyzed  in  this  thesis.  The  output  of  this  file  is  a  graph  of  the 
detector’s  decision  statistic  along  with  the  ‘C’  and  ‘P’  windows. 

clear  all 

%  this  file  implements  the  Delay  and  Correlate  detector.  It  uses  the  functions  st  sym  and 
%  short_tmg_samples. 

%  author:  Vik  Sardana,  MSEE 
%  22  DEC  2006 

%  output  of  this  file  is  a  graph  of  the  decision  statistic,  P  and  C  windows. 


%  get  short  training  symbols 

ml=st_sym; 

m2=ifft(ml); 

% 

%  add  noise  to  the  delay  buffer  (or  else  divide  by  zero)  and  training  samples 
noise_variance=6.3477e-4;  %average  signal  power  is  6.3477e-3 
nsd=sqrt(noise_variance);  %standard  deviation  of  noise,  used  in  randn  function 

% 

%cyclically  extend  short  trng  symbols  to  161  samples  and  add  noise 
rn3=st_trng_samples(rn2)+nsd*randn(  1 , 1 6 1  )+nsd*randn(  1 , 1 6  l)*j ; 


buffi  =nsd*randn(  1 , 1 000)+nsd*randn(  1 , 1 000)  *j ; 
buff2=nsd*randn(  1 , 1 600)+nsd*randn(  1 , 1 600)  *j ; 

% 

%  %create  incoming  signal  with  noise  **short  trng  samples  start  at  sample  1000  ** 
rx_signal=  [buffi  m3  buff2]; 

% 

% 

% 

%  search_window=1600;  %scan  from  samples  1  to  1600 
%  D=16;  %size  of  delay  window  =  periodicity  of  short  trng  samples 
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D=16;  %size  of  delay  window,  should  equal  periodicity  of  training  samples 
L=16;  %number  of  multiple  samples  to  use  to  compute  decision  statistic 
nn=l+D; 

rx_len=length(rx_signal) ; 

rx_cross=rx_signal(nn:rx_len).*conj(rx_signal(nn-D:rx_len-D)); 
rx_pwr=(abs(rx_signal(nn-D:rx  len-D))).A2; 


%  compute  moving  sums  of  rx  power  and  cross  corr  using  delayed  signal  and  L  samples 

n=l; 

% 

while  n<=length(rx_pwr)-L-l 

ms_rx_p  wr(n)=sum(rx_p  wr(n :  n+L- 1 )) ; 
ms_rx_cross(n)=abs(sum(rx_cross(n:n+L- 1 ))); 

ds(n+L)=(ms_rx_cross(n)A2)/(ms_rx_pwr(n))A2;  %compute  decision  statistic 
n=n+l; 
end 


% 

tt=[  1 : 1  :length(ms_rx_pwr)] ; 

% 

figure(2) 

subplot(3,l,l) 

plot(tt,ms_rx_pwr) 

grid 

axis([900  1200  min(ms_rx_pwr)  max(ms_rx_pwr)]) 
title(’Moving  sum  Rx  power') 

% 

subplot(3,l,2) 

plot(tt,ms_rx_cross) 

grid 

axis([900  1200  min(ms_rx_cross)  max(ms_rx_cross)]) 
title(’Moving  sum  Cross  Correlation') 

% 

vv=[l:l:length(ds)]; 

subplot(3,l,3) 

plot(vv,ds) 

grid 

%  axis([900  1200  min(ds)  max(ds)]) 
title('Decision  statistic') 
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st_sym.m  -  This  function  is  used  to  generate  the  short  training  symbols  used  by  the 
Delay  and  Correlate  detector  and  the  Matched  Filter  detector.  This  function  also 
performs  the  necessary  coefficient  mapping  as  shown  in  Figure  18. 
function[mapped_st_trng_symbols]=st_sym 

%  this  function  returns  the  row  vector  of  mapped  short  tmg  symbols  per  the 
%  802.1  la  standard 


%short  training  symbols:(-26:26)  ref  IEEE  802.1  la  std  17.3.3 

st_tmg_symbols_initial=sqrt(13/6)*[0  0  1+j  0  0  0  -1-j  0  0  0  1+j  0  0  0  -1-j  0  0  0  -1-j  0  0  0 
1+j  0  0  0  0  0  0  0  -1-j  0  0  0  -1-j  0  0  0  1+j  0  0  0  1+j  0  0  0  1+j  0  0  0  1+j  0  0]; 

%  mapped  short  training  symbols  per  figure  109,  802.1 1  standard  (inputs  of  the  idft) 
mapped_st_trng_symbols=[0  st_trng_symbols_initial(28:53)  00000000000 
st_tmg_symbols_initial(  1 :26)]; 

short_trng_samples.m  -  This  function  creates  the  161  values  of  the  short  training 
samples  and  applies  the  windowing  function  per  the  IEEE  standard  [4]. 

function[mapped_short_trng_symbols]=st_trng_samples(short_trng_symbols) 

%  this  function  cyclically  extends  the  64  short  training  symbols  to  161 
%  samples  per  the  802. 1 1  standard  (REFitable  G.3) 

mapped_short_tmg_symbols=[short_trng_symbols  short_tmg_symbols 
short_tmg_symbols(  1 :33)]; 

%  802.1 1  windowing  function  is  also  applied  here 
mapped_short_tmg_symbols(  1  )=mapped_short_trng_symbols(  1 )  *  .5 ; 
mapped_short_trng_symbols(  1 6 1  )=mapped_short_tmg_symbols(  1 6 1 )  *  .5 ; 

match_filter_detector.m  -  This  m-file  implements  the  matched  filter  detector  that  was 
extensively  analyzed  in  this  thesis.  The  output  of  this  file  is  a  graph  of  the  detector’s 
decision  statistic  and  a  depiction  of  the  embedded  short  training  samples  in  an  AWGN 
channel. 

clear  all 

%  this  file  is  the  match  filter  detector.  It  uses  the  functions  st  sym  and 
%  short_tmg_samples. 

%  author:  Vik  Sardana,  MSEE  student 
%  22  DEC  2006 


%  get  short  training  symbols 
ml=st_sym; 
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rn2=ifft(rnl); 

% 

%  add  noise  to  the  delay  buffer  (or  else  divide  by  zero)  and  training  samples 
noise_variance=6.3477e-4;  %average  signal  power  is  6.3477e-3 
nsd=sqrt(noise_variance);  %standard  deviation  of  noise,  used  in  randn  function 


%cyclically  extend  short  trng  symbols  to  161  samples  and  add  noise 
sig_plus_noise=st_trng_samples(m2)+nsd*randn(  1 , 1 6 1  )+nsd*randn(  1,161  )*j; 

%  sig_plus_noise=st_tmg_samples(m2);  %test  only  NO  NOISE  for  now 

%create  1600  samples  of  noise 

buffi  =nsd*randn(  1 , 1 000)+nsd*randn(  1 , 1 000)*j ; 

buff2=nsd*randn(  1 , 1 600)+nsd*randn(  1 , 1 600)*j ; 

% 

%create  incoming  signal  with  noise  **short  trng  samples  start  at  sample  1000  ** 
rx_signal=[buffl  sig_plus_noise  buff2]; 


% 


D=16;  %size  of  delay  window,  should  equal  periodicity  of  training  samples 


short_trng_sample_values=st_tmg_samples(rn2); 
stored_values=conj(short_tmg_sample_values(  17:32)); 


%compute  matched  filter  decision  stat 

nn=[l :  1  :length(rx_signal)-(D+l)]; 
ds=zeros(  1  ,length(nn)); 

nn=l; 

while  nn<=length(rx_signal)-D- 1 

ds(nn)=abs(sum(rx_signal(nn:nn+(D-l)).*stored_values)); 

nn=nn+l; 

end 


%  length(rx_signal)-D- 1 
tt=[l:l:length(ds)]; 
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gg=[l :  1 :  length(rx_signal)] ; 

figure(l) 

subplot(2,l,2) 

plot(tt,ds) 

grid 

title('Matched  Filter  Output') 

ylabel('Decision  Statistic') 

axis([900  1200  0  .25]) 

subplot(2,l,l) 

plot(gg,rx_signal) 

grid 

ylabel('Received  Signal  Amplitude') 

ml_cfo_corrector_withLTS.m-  This  m-filc  implements  the  maximum  likelihood  (ML) 
frequency  synchronizer  that  was  extensively  analyzed  in  this  thesis.  This  file  uses  only  4 
short  training  symbols  and  the  long  training  symbols  to  evaluate  the  algorithm’s 
“robustness”.  The  number  of  the  training  symbols  and  the  SNR  can  be  modified  at  the 
user’s  discretion.  The  output  of  this  file  is  a  series  of  histograms  that  depict  the  final 
frequency  offset  percentage  (relative  to  sub-carrier  spacing)  for  each  trial  run  of  the  test 
signal. 

%  This  m-file  represents  the  data-aided  maximum  likelihood  frequency  offset  corrector 
%  for  802.1  lg  OFDM  systems 

%  author  Vik  Sardana 
%  used  for  author's  MSEE  thesis  2006 
clear  all 

%  step  1  :  get  short  training  symbols  to  be  used  in  data-aided  cfo  correction 
short_training_samples_64=st_sym; 

%  step  2:  take  ifft  in  the  transmitter 

post_ifft_short_training_samples=ifft(short_training_samples_64); 

%  step  3:  cyclically  extend  64  short  tmg  samples  to  161  samples 
short_training_samples_  161  =st_tmg_samples(post_ifft_short_training_samples); 

%  step  4:  step  up  short  tmg  samples  to  transmit  carrier  frequency 
tx_packet=transmitted_signal(short_training_samples_  161); 

%  tx_packet=[tx_packet(l:64)  zeros(l,97)]; 
tx_packet=tx_packet(l  :64); 
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%  step  5:  corrupt  transmitted  signal  with  noise  in  channel 

%  noise_variance=6.3477e-4;  %  for  10  dB  -  average  signal  power  is  .10 156/16=6. 3477e- 
3 

%  noise_varianee=.00 159447;  %  for  6  dB 
noise_variance=. 003 17385;  %. 00317385  for  3  dB 

%  noise_variance=6.3477e-3;  %  for  0  dB 

%  noise_variance=. 02527065;  %. 02525871  for  -6  dB 

nsd=sqrt(noise_variance);  %standard  deviation  of  noise,  used  in  randn  function 


%  step  5a:  loop  to  collect  data  on  frequency  offset  with  correction 

loop_count=l; 

while  loop_count<=500 

tx_packet_plus_channel=tx_packet+nsd*randn(  1  ,length(tx_packet))+nsd*randn(  1  ,length( 
tx_packet))*j; 

CFO_correction_short_symbols(loop_count)=received_signal_with_correction(tx_packet 

pluschannel); 

frequency_offset_with_correction(loop_count)=long_symbol_corrector... 

(CFO_correction_short_symbols(loop_count),nsd); 

rx_packet=received_signal(tx_packet_plus_channel); 

frequency_offset_without_correction(loop_count)=cfo_calculator_short_symbols(rx_pac 

ket); 

loop_count=loop_count+ 1 ; 
end 


bins=[30:.5:50]; 

binsl=[-10:.5:10]; 

figure(l) 

hist(CFO_correction_short_symbols,binsl) 

grid 

xlabel('Final  frequency  offset  percent  relative  to  sub-carrier  spacing') 
title(’WITHOUT  long  training  samples') 

figure(3) 

hist(frequency_offset_with_correction,bins  1 ) 
grid 

xlabel('Final  frequency  offset  percent  relative  to  sub-carrier  spacing') 
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title(’With  long  training  samples') 


mean_percent_offset=mean(frequency_offset_with_correction) 

var_percent_offset=var(frequency_offset_with_correction) 

mean_percent_offset_uncorrected=mean(frequency_offset_without_correction) 

var_percent_offset_uncorrected=var(frequency_offset_without_correction) 

figure(4) 
subplot(  1,2,1) 

hist(CFO_correction_short_symbols,binsl) 

grid 

xlabel('Final  frequency  offset  percent  relative  to  sub-carrier  spacing') 
title(’WITHOUT  long  training  samples') 
ylabel('Trials') 
subplot)  1,2,2) 

hist(frequency_offset_with_correction,bins  1 ) 
grid 

xlabel('Final  frequency  offset  percent  relative  to  sub-carrier  spacing') 

title(’With  long  training  samples') 

ylabel(’Trials') 

nls_receiver_with_correction6.m  -  This  m-file  implements  the  non-linear  least  squares 
(NLS)  frequency  synchronizer  that  was  extensively  analyzed  in  this  thesis.  This  file  uses 
only  4  short  training  symbols  and  the  long  training  symbols  to  evaluate  the  algorithm’s 
“robustness”.  The  number  of  the  training  symbols  and  the  SNR  can  be  modified  at  the 
user’s  discretion.  The  output  of  this  file  is  a  histogram  that  depicts  the  final  frequency 
offset  percentage  (relative  to  sub-carrier  spacing)  for  each  trial  run  of  the  test  signal. 

%  this  is  an  m-file  which  implements  the  NLS 
%  frequency  synchronization  algorithm  as  described  in  my  thesis 
%  author:  Vik  Sardana,  DEC  2006 

%  THIS  PARTICULAR  VERSION  USES  ONLY  4  SHORT  TRAINING  SYMBOLS 
clear  all 

%  step  1  :  get  short  training  symbols  to  be  used  in  data-aided  cfo  correction 
short_training_samples_64=st_sym; 

%  step  2:  take  ifft  in  the  transmitter 

post_ifft_short_training_samples=ifft(short_training_samples_64); 

%  step  3:  cyclically  extend  64  short  tmg  samples  to  161  samples 

short_training_samples=st_tnig_samples(post_ifft_short_training_samples); 

%now  get  long  training  samples 
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long_symbols=lg_sym; 

long_training_samples=lt_sym_time_domain; 

%discard  first  32  samples  of  long  training  samples  (it  is  the  CP) 
long_training_samples_no_cp=long_training_samples(33 : 1 60); 

%  add  noise  to  simulate  channel 
noise_variance=. 003 17385  %3  dB  SNR 
%  noise_varianee=.00 159372  %6  dB  SNR 

%  noise_variance=6.3447e-4;  %10  dB  SNR  %average  signal  power  is  6.3447e-3 
nsd=sqrt(noise_variance);  %standard  deviation  of  noise,  used  in  randn  function 

%  add  noise  to  the  short  training  samples 

received_noisy_samples=short_training_samples+nsd*randn(  1 , 1 6  l)+nsd*randn(  1 , 1 6  l)*j 


%add  noise  to  the  long  training  samples 

noisy_long_samples=long_training_samples_no_cp+nsd*randn(  1 ,  128)+nsd*randn(  1 , 128 

)*i; 

%create  circular  convulation  matrix  from  long  training  samples 
stored_long_training_samples=[long_training_samples_no_cp 
longtrainingsamplesnocp  longtrainingsamplesnocp . . . 

longtrainingsamplesnocp  long  training  samples  no  cp 
long  training  samples  no  cp  long  training  samples  no  cp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 
long  training  samples  no  cp  long  training  samples  no  cp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 
long  training  samples  no  cp  long  training  samples  no  cp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 

long  training  samples  no  cp  long  training  samples  no  cp 

longtrainingsamplesnocp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 

long  training  samples  no  cp  long  training  samples  no  cp 

longtrainingsamplesnocp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 

long  training  samples  no  cp  long  training  samples  no  cp 

longtrainingsamplesnocp . . . 

long  training  samples  no  cp  long  training  samples  no  cp 

long  training  samples  no  cp  long  training  samples  no  cp] ; 

nn=l; 

Dlong=64; 
for  col=[l :  1  :Dlong] 
for  row=[  1 : 1  :Dlong] 
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XL(row,col)=stored  long_training_samplcs(nn);  %circular  convultion  matrix  for 
long  tmg  samples 
nn=nn+l; 
end 

nn=nn-l; 

end 

XL_inv=pinv(XL);  %creates  psuedo  inverse  matrix  of  long  training  symbol  convultion 
matrix 


%now  simulate  frequency  offset 

%  step  1 :  initialize  variables 

ts=50e-9;  %sample  time 

f_tx=2.9e9;  %  transmitted  carrier  frequency 

freq_spacing=3 1 2 . 5  e3 ; 

D=16;  %periodicity  of  the  short  training  samples 

carrier_freq_offset=freq_spacing*-.40;  %freq  offset  of  -40%  of  sub-carrier  spacing 
f_rx=f_tx+carrier_freq_offset;  %  received  carrier  frequency=tx_freq  +  some  fireq  offset. 
f_delta=f_rx-f_tx; 


n=l; 

while  n<=161 

yn(n)=received_noisy_samples(n).*(exp(j*2*pi*ts*n*f_delta));  %creating  fireq 
offset  for  noisy  samples 
n=n+l; 
end 


%  step  3a:  now  form  circular  convolution  matrix  (a  16x16  matrix) 

%  based  on  short  trng  sample  values.  Convultion  matrix  called  ’X’ 
stored_values=short_training_samples(  17:48); 

rxx=[stored_values  storedjvalues  stored_values  stored_values  storedjvalues 
stored_values...  %buffering  short  training  samples  to  create  matrix 
stored_values  stored_values  stored_values  stored_values]; 

nn=l; 

for  col=[l:l:D] 
for  row=[l:l:D] 

X(row,col)=rxx(nn); 

nn=nn+l; 

end 

nn=nn-l; 

end 
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X_inv=pinv(X);  %creates  psuedo  inverse  matrix  of  short  training  symbols 

%  implement  loop  here  to  compute  distribution  of  final  percentage  offset 

no_of_trials=l; 

while  no_of_trials<=500 

received_noisy_samples=short_training_samples+nsd*randn(l,length(received_noisy_sa 

mples))... 

+nsd*randn(  1  ,length(received_noisy_samples))*j ; 
n=l; 

while  n<=length(received_noisy_samples) 

m(n)=received_noisy_samples(n).*(exp(j*2*pi*ts*n*f_delta));  %creating  freq  offset 
for  noisy  samples 
n=n+l; 
end 

symboll=rn(l:16);  %create  offset  symbols  1  and  2 
received_sig=conj  (symbol  1 )'; 
n=  1 7; 

symbol_number=  1 ; 

while  symbol_number<=4  %CHANGE  #1  FROM  10  TO  4  TO  TEST  FOR 
ROBUSTNESS  (4  SHORT  SYMBOLS  ONLY) 
xx=l/symbol_number; 
h_hat=X_inv*receivcd_sig*xx; 

H_HAT=fft(h_hat); 
mean_H_HAT=mean(H_HAT); 
angle=atan(imag(mean_H_HAT)/real(mean_H_HAT)); 
fo=angle/(2*pi*ts*D); 
f  corfisymbol  number)=- 1  *fo; 

percent_corrected(symbol_number)=100*(f_delta+f_corr(symbol_number))/freq_spacing 

nn=l; 

initial_sample=n; 

if  n<=64  %  CHANGE  #2  FROM  1 60  TO  64 

while  n<initial_sample+D 

next_symbol_corrected(nn)=received_noisy_samples(n).*(exp(j*2*pis|sts*n*(f_delta+f_c 

orr(symbol_number)))); 

rn_corrected(n)=received_noisy_samples(n).*(exp(j*2*pi*tsH5n*(f_delta+f_corr(symbol_ 

number)))); 

n=n+l; 

nn=nn+l; 

end 

end 
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received_sig=received_sig+conj(next_symbol_corrected)'; 
symbol_number=symbol_number+ 1 ; 
end 

final_percent_corrected_st(no_of_trials)=percent_corrected(length(percent_corrected)); 
n=[l :  1 : 128];  %now  apply  short  samples  freq  correction  to  long  trng  samples 

noisy_long_samples=long_training_samples_no_cp+nsd*randn(  1 ,  128)+nsd*randn(  1 , 128 

)*j; 

m_long=noisy_long_samples(n).*(exp(j*2*pi*ts*n*(f_delta+f_corr(symbol_number- 

i)))); 

received_long_samples=conj(rn_long(l  :64)+rn_long(65: 128))'; 
h_est= .  5  *XL_inv  *  rec  eived_long_samples ; 

H_est=fft(h_est); 

H_est_mean=mean(H_est); 

angle_long=atan(imag(H_est_mean)/real(H_est_mean)); 
folo  n  g=a  n  g  1  c  1  o  n  g/(  2  *pi*ts  *Dlong); 

final_freq_offset_correction=f_corr(symbol_number- 1  )-fo_long; 
f_corr(symbol_number)=final_freq_offset_correction; 

fmal_offset_percentage(no_of_trials)=100*(f_delta+final_freq_offset_correction)/freq_s 
pacing;  %CHANGE  4  COMMENT  LINE  145,  ADD  LINE  146 
no_of_trials=no_of_trials+ 1 ; 
end 

figure(2) 

bins=[-20:.5:20]; 

hist(final_offset_percentage,bins) 

grid 

xlabel('Frequency  offset  percentage  relative  to  sub-carrier  spacing') 
ylabel(’Trials') 

lg_sym.ni  -  this  is  a  function  used  to  create  the  long  training  symbols. 
function[mapped_lg_trng_symbols]=lg_sym 

%  this  function  returns  the  row  vector  of  mapped  long  trng  symbols  per  the 
%  802.1  la  standard 


%long  training  symbols:(-26:26)  ref  IEEE  802.11a  appendix  G.3.2  pg  59 
lg_trng_symbols_initial=[  1  1-1-111-11-1111111-1-111-11-111110  1-1- 
111-11-11-1-1-1-1-111-1-11-11-1111  1]; 

%  mapped  short  traning  symbols  per  figure  109,  802.1 1  standard  (inputs  of  the  idft) 
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mapped_lg_trng_symbols=[0  lg_trng_symbols_initial(28:53)  00000000000 
lg_trng_symbols_initial(  1 :26)]; 

lg_sym_time_domain.m  -  this  function  creates  the  time  domain  representation  of  the 
long  training  samples.  It  also  calls  the  function  cyclic_prefix_creator.m 

function[long_training_samples]=lg_sym_time_domain 
%  Vik  Sardana  MSEE  thesis 

%  this  function  used  to  generate  the  long  training  samples  in  the  time  domain 

%generate  long  training  samples: 

long_training_symbols=lg_sym; 

post_ifft_long_tmg_symbols=ifft(long_training_symbols); 
long_training_samples=cyclic_prefix_creator(post_ifft_long_trng_symbols); 
cyclic_prefix_creator.m  -  this  function  creates  the  cyclic  prefix  for  the  long  training 
symbols  and  also  windows  the  first  and  last  sample  in  accordance  with  the  standard. 

function[output_samples]=cyclic_prefix_creator( input) 

%  this  function  cyclically  extends  the  64  post  ifft  inputs  to  161 
%  samples  per  the  802. 1 1  standard  (REFitable  G.3) 

output_samples=[input(33:64)  input  input  input(l)]; 

%  802.1 1  windowing  function  is  also  applied  here 
output_samples(  1  )=output_samples(  1  )*  .5; 
output_samples(  1 6 1  )=output_samples(  1 6 1 )  * .  5 ; 
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